R/sample_inc_rate.R
sample_inc_rate.RdThe number of cases are back calculated using the incidence rate and the size of the population. Samples are then taken of notifications, assuming a poisson distribution. Incidence rates are then calculated based on these samples. Multiple samples of a matrix input may be taken with the output formated as a list.
sample_inc_rate(inc_rate, pop, inc_per = 1e+05, samples = 100)
| inc_rate | Incidence rates formated as |
|---|---|
| pop | Population estimates for each cohort as a vector. |
| inc_per | The scaling factor for the incidence rate, defaults to 100,000. |
| samples | Numeric, the number of samples to take. Defaults to 100. |
Incidence rate samples as a list, formated as a 3 by 3 matrix.
sample_inc_rate(Unvaccinated, sutherland_data$BCG.Unvaccinated[1:3], samples = 5)#> Error in sample_inc_rate(Unvaccinated, sutherland_data$BCG.Unvaccinated[1:3], samples = 5): object 'Unvaccinated' not found