
Weighted post-stratification estimation of coverage over several service delivery units
Source:R/05-post_strat_estimator.R
estimate_coverage.Rd
Weighted post-stratification estimation of coverage over several service delivery units
Usage
estimate_coverage_overall(cov_df, pop_df, strata, u5, p, k = 3)
estimate_coverage(cov_df, cov_type = c("cf", "tc"), k = 3)
Arguments
- cov_df
A
data.frame()
of stratified coverage survey data to get overall coverage estimates of.cov_df
should have a variable namedcases_in
for number of SAM or MAM cases in the programme found during the survey,cases_out
for number SAM or MAM cases not in the programme found during the survey, andrec_in
for children recovering from SAM or MAM who are in the programme found during the survey. A final required variable should be one that contains identifying geographical information corresponding to the location from which each row of the survey data was collected from.- pop_df
A
data.frame()
with at least two variables:strata
for the stratification/grouping information that matches the grouping information incov_df
andpop
for information on population for the given grouping information.- strata
A character value of the variable name in
cov_df
that corresponds to thestrata
values to match with values inpop_df
.- u5
A numeric value for the proportion of the population that is under years old.
- p
Prevalence of SAM or MAM in the given population.
- k
Correction factor. Ratio of the mean length of an untreated episode to the mean length of a CMAM treatment episode
- cov_type
Coverage estimator to report. Either "cf" for case-finding effectiveness or "tc" for treatment coverage. Default is "cf".
Value
A list of overall coverage estimates with corresponding 95% confidence intervals for case-finding effectiveness and treatment coverage.
Examples
cov_df <- survey_data
pop_df <- pop_data |>
setNames(nm = c("strata", "pop"))
estimate_coverage_overall(
cov_df, pop_df, strata = "district", u5 = 0.177, p = 0.01
)
#> $cf
#> $cf$estimate
#> [1] 0.1257481
#>
#> $cf$ci
#> [1] 0.09247579 0.15902045
#>
#>
#> $tc
#> $tc$estimate
#> [1] 0.1706466
#>
#> $tc$ci
#> [1] 0.1371647 0.2041284
#>
#>