De typische vertraging van testafname tot rapportage is 2 dagen: 1 voor het lab en 1 voor het dagrapport. Dan wordt iemand bijvoorbeeld op een maandagmiddag getest, komt er ergens op dinsdag uitslag uit het lab en zit hij woensdag in het dagrapport (dat van di 10:00 - wo 10:00 loopt).
De grootste groep in het cijfer van maandag is dus op zaterdag (dag -2) getest. Maar: ook de tests van zondagochtend kunnen voor nog in het maandagrapport terechtkomen.
Onderzoeksvraag: "wat is de invloed van de zondagse tests op het cijfer van maandag?"
Ik neem de daggegevens van de afgelopen 7 weken (sinds de piek van de golf in december).
Per dag ga ik het aantal gerapporteerde besmettingen op de dag X vergelijken met het aantal in de rollende week (dag X-3 ... dag X+3).
Door een rolling window te gebruiken dat symmetrisch naar links en rechts kijkt, corrigeer je voor de skew in de data. Die skew is het effect dat de eerste dagen in de week een te hoge factor en de latere dagen in de week een te lage factor zouden krijgen vanuit de neergaande trend die over de weken heen loopt.
Dit zijn de getallen:
sun 0.143823
mon 0.122910
tue 0.122125
wed 0.146157
thu 0.162670
fri 0.152967
sat 0.149348
Je ziet dat maandag en dinsdag de laagste factor hebben. Dat komt door een lager aantal testen op zaterdag en zondag.
Als ik de factoren van de "andere dagen" (wo t/m zo) middel, dan kom op op een gemiddelde niet-weekend dagfactor van 0.151
De "delta" kolom is de dagfactor minus die gemiddelde-dag-zonder-weekendeffect-factor:
sun 0.143823 -0.007170
mon 0.122910 -0.028083
tue 0.122125 -0.028868
wed 0.146157 -0.004836
thu 0.162670 0.011677
fri 0.152967 0.001974
sat 0.149348 -0.001645
Je zit op zondag ongeveer een kwart van het weekend-effect (-0.007), en op ma en het volle effect (-0.028).
Uitkomst: Daaruit zou je kunnen schatten dat tests die al met 1 dag vertraging in het rapport zitten i.p.v. 2 dagen, ongeveer een kwart van het totaal betreffen.
Toepassing:
Leg ik deze wijsheid op een willekeurige maandag, dan zou je kunnen stellen dat 75% van een maandagse rapportage van de zaterdag afkomstig is en 25% van de zondag. Vandaag is het gecorrigeerde cijfer 2221.
Als ik omhoog reken om de ontbrekende tests van zondag te verdisconteren (geschat op 25% van het normale volume op een maandag), dan kom ik op omstreeks 2960 nieuwe infecties als er gisteren wel gewoon getest was. Dat cijfer ligt 7% lager dan 1 week geleden.
De marges zijn breed vanwege mijn lompe extrapolatie, maar een structurele omekeer naar een R>1 wordt in elk geval niet bevestigd door dit getal.
[
Voor 5% gewijzigd door
t_captain op 08-02-2021 20:44
]