km.ci | R Documentation |
Computes pointwise and simultaneous confidence intervals for the Kaplan-Meier estimator.
km.ci(survi, conf.level = 0.95, tl = NA, tu = NA, method = "rothman")
survi |
A survival object for which the new confidence limits should be computed. This can be built using the "Surv" and the "survfit" function in the R package "survival". "km.ci" modifies the confidence limits in this object. |
conf.level |
The level for a two-sided confidence interval on the survival curve. Default is 0.95. |
tl |
The lower time boundary for the simultaneous confidence limits. If it is missing the smallest event time is used. |
tu |
The upper time boundary for the simultaneous confidence limits. If it is missing the largest event time is used. |
method |
One of '"peto"', '"linear"', '"log"', "loglog"', '"rothman"', "grunkemeier"', '"hall-wellner"', '"loghall"', "epband"', "logep" |
A simulation study showed, that three confidence intervals produce satisfying confidence limits. One is the "loglog" confidence interval, an interval which is based on the log of the hazard. The other competitive confidence concept was introduced by Rothman (1978) and is using the assumption that the survival estimator follows a binomial distribution. Another good confidence concept was invented by Thomas and Grunkemeier (1975) and is derived by minimizing the likelihood function under certain constraints. Special thanks goes to Robert Gentleman for providing code for the confidence interval by Thomas and Grunkemeier.
The confidence interval using Peto's variance can not be recommended since it yields confidence limits outside the admissible range [0;1] as well as the "linear" and the "log" (which is based on the logarithm of S(t)).
The function can produce simultaneous confidence bands, too. The Hall-Wellner band (1980) and the Equal Precision band by Nair (1984) together with their log-transformed counterpart. From all simultaneous confidence intervals only the log-transformed Equal Precision "logep" band can be recommended. The limits are computed according to the statistical tables in Klein and Moeschberger (2002).
a 'survfit' object;
see the help on 'survfit.object' for details.
Strobl, R.
Strobl, R., Dirschedl, P. and Mansmann, U.. Comparison of simultaneous and pointwise confidence intervals for survival functions. (2005, submitted to Biom. J.).
survfit
,
print.survfit
, plot.survfit
,
lines.survfit
,
summary.survfit
,
survfit.object
, coxph
,
Surv
, strata
.
require(survival) data(rectum.dat) # fit a Kaplan-Meier and plot it fit <- survfit(Surv(time, status) ~ 1, data=rectum.dat) plot(fit) fit2 <- km.ci(fit) plot(fit2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.