Fay, Michael (NIH/NIAID) [E]
2013-Jun-05 15:39 UTC
[R-pkgs] bpcp package for pointwise confidence intervals for a survival distribution
Hi all, I just uploaded a new version of the bpcp package. It calculates confidence intervals for a survival distribution for right-censored data using the newly developed beta product confidence procedure. Previously developed methods can have substantial error rate inflation for the lower limit, especially at the right end of the curves when there are small numbers of events. The bpcp method guarantees coverage when there is no censoring (when it reduces to the exact Clopper-Pearson intervals for binomial data at each time) or when there is progressive type II censoring. Simulations show at least nominal coverage for independent censoring. Also, it is fast so it can be used routinely. For a complete description of the new method see Fay, Brittain, and Proschan (2013). Pointwise confidence intervals for a survival distribution with small samples or heavy censoring. Biostatistics doi:10.1093/biostatistics/kxt016 (available at http://www.niaid.nih.gov/about/organization/dcr/brb/staff/Pages/michael.aspx ). Mike ****************************************************************** Michael P. Fay, PhD Mathematical Statistician National Institute of Allergy and Infectious Diseases Tel: 301-451-5124 Fax:301-480-0912 (U.S. postal mail address) 6700B Rockledge Drive MSC 7609 Bethesda, MD 20892-7609 (Overnight mail address) 6700-A Rockledge Drive, Room 5133 Bethesda, MD 20817 http://www3.niaid.nih.gov/about/organization/dcr/BRB/staff/michael.htm ********************************************************************** Disclaimer: The information in this e-mail and any of its attachment...{{dropped:7}}
Reasonably Related Threads
- Update to exact2x2 package: Exact McNemar's test
- Update to exact2x2 package: Exact McNemar's test
- R package: Mchtest - Monte Carlo hypothesis testing allowing Sequential Stopping
- R package: Mchtest - Monte Carlo hypothesis testing allowing Sequential Stopping
- poisson.test from stats package does not pass the conf.level (PR#14195)