[R] survey package: proportion estimate confidence intervals using svymean

From: Peter Holck <holck_at_hawaii.edu>
Date: Tue, 12 Feb 2008 23:04:07 -1000


Using the survey package I find it is convenient and easy to get estimated proportions using svymean, and their corresponding estimated standard errors. But is there any elegant/simple way to calculate corresponding confidence intervals for those proportions? Of course +/- 1.96 s.e. is a reasonable approximation for a 95% CI, but (incorrectly) assumes symmetrical distribution for a proportion.

About all I can think of is to transform onto some scale like the logit scale, use the delta method to estimate standard errors, and transform back for appropriate confidence intervals. Or perhaps use the machinery of the survey package svyglm function with logit link (family=quasibinomial) to do the calculations for me. Still, seems rather cumbersome. Is there a better solution?

Thanks for any suggestions or pointing out the errors in my logic! Peter Holck



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Wed 13 Feb 2008 - 09:11:47 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 13 Feb 2008 - 16:30:13 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive