Constructs confidence intervals on directly standardized rates via several parameterizations.
The DESCRIPTION file: This package was not yet installed at build time.
Index: This package was not yet installed at build time.
Maintainer: Michael Nelson <[email protected]>
Fay & Feuer (1997). 'Confidence intervals for directly standardized rates: a method based on the gamma distribution.' Statistics in Medicine. 16: 791–801. https://doi.org/10.1002/(SICI)1097-0258(19970415)16:7<791::AID-SIM500>3.0.CO;2-%23
Fay & Kim (2017). 'Confidence intervals for directly standardized rates using mid-p gamma intervals.' Biometrical Journal. doi: 10.1002/bimj.201600111.
Ng, Filardo, & Zheng (2008). 'Confidence interval estimating procedures for standardized incidence rates.' Computational Statistics and Data Analysis 52: 3501–3516. doi: 10.1016/j.csda.2007.11.004
Tiwari, Clegg, & Zou (2006). 'Efficient interval estimation for age-adjusted cancer rates.' Statistical Methods in Medical Research 15: 547–569. doi: 10.1177/0962280206070621
Swift, MB (1995). 'Simple confidence intervals for standardized rates based on the approximate bootstrap method', Statistics in Medicine, 14, 1875–1888. doi: 10.1002/sim.4780141704.
Dobson, AJ, Kuulasmaa, K, Eberle, E and Scherer, J (1991) 'Confidence intervals for weighted sums of Poisson parameters', Statistics in Medicine, 10: 457–462. doi: 10.1002/sim.4780100317
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
## Not run: ## using the example from asht::wspoissonTest ## birth data on Down's syndrome from Michigan, 1950-1964 ## see Table II of Fay and Feuer (1997) ## xfive = counts for mothers who have had 5 or more children ## nfive and ntotal are number of live births xfive<-c(0,8,63,112,262,295) nfive<-c(327,30666,123419,149919,104088,34392) ntotal<-c(319933,931318,786511,488235,237863,61313) ## use mult =10^5 to give rates per 100,000 ## gamma method of Fay and Feuer (1997) is default dsr(xfive, nfive, ntotal, ci.method = "gamma", type = "max", mult = 1e5) ## Dobson et al (1991) dsr(xfive, nfive, ntotal, ci.method = "moments", type = "dobson", mult = 1e5) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.