This function computes the mean residual life for censored data
using the Kaplan-Meier estimate of the survival function. If
When there are ties between censored and uncensored observations, for definiteness our ordering places the censored observations before uncensored.
This function is used by locfit.censor
to compute
censored regression estimates.
km.mrl(times, cens)
A vector of the estimated mean residual life. For uncensored observations, the corresponding estimate is 0.
Obsereved survival times.
Logical variable indicating censoring. The coding is 1
or TRUE
for censored; 0
or FALSE
for uncensored.
Buckley, J. and James, I. (1979). Linear Regression with censored data. Biometrika 66, 429-436.
Loader, C. (1999). Local Regression and Likelihood. Springer, NY (Section 7.2).
locfit.censor
# censored regression using the Kaplan-Meier estimate.
data(heart, package="locfit")
fit <- locfit.censor(log10(surv+0.5)~age, cens=cens, data=heart, km=TRUE)
plotbyfactor(heart$age, 0.5+heart$surv, heart$cens, ylim=c(0.5,16000), log="y")
lines(fit, tr=function(x)10^x)
Run the code above in your browser using DataLab