estimate_asympto {asymptor}

Estimate the proportion of asymptomatic cases by capture/recapture

Description

Estimate the proportion of asymptomatic cases by capture/recapture

Usage

estimate_asympto(date, cases, deaths, bounds = c("lower", "upper"))

Arguments

date

A vector containing the dates

cases

A numeric vector containing the number of new cases at each date (not the cumulative number of cases).

deaths

A numeric vector containing the number of new deaths at each date (not the cumulative number of cases).

bounds

"lower", "upper", or both ⁠c("lower, "upper"⁠) (the default), telling which bounds of the number of asymptomatic cases are computed. associated with the lower and upper bounds. estimation. FALSE is never the recommended setting excepted for testing or debugging purpose.

Value

A data.frame with two or three columns (depending on the value of the bounds argument):

• date: the original date column

• lower: the lower bound of asymptomatic cases

• upper: the upper bound of asymptomatic cases

Note

Please note there is currently some discussion about the validity of the methods implemented in this package. You should read carefully the original articles, alongside this answer from Li et al. (2022) before using this package in your project.

References

Böhning D., Rocchetti I., Maruotti A., Holling H. (2020), Estimating the undetected infections in the Covid-19 outbreak by harnessing capture–recapture methods, International Journal of Infectious Diseases, 97, p197-201, doi:10.1016/j.ijid.2020.06.009.

Rocchetti I., Böhning D., Holling H., Maruotti A., (2020), Estimating the size of undetected cases of the SARS-CoV-2 outbreak in Europe: An upperbound estimator, medRxiv, doi:10.1101/2020.07.14.20153445.

Li M., Dushoff J., Earn D.J.D., Bolker B.M. (2022), Evaluating undercounts in epidemics: responsi to Maruotti et al. 2022, arXiv, doi:10.48550/ARXIV.2209.11334.

Examples

d <- readRDS(system.file("extdata", "covid19_italy.rds", package = "asymptor"))
head(d)

estimate_asympto(d$date, d$new_cases, d$new_deaths)