KMW | R Documentation |
This function returns a vector with the Kaplan-Meier weights.
KMW(time, status)
time |
Survival time of the process. |
status |
Censoring indicator of the survival time of the process; 0 if the survival time is censored and 1 otherwise. |
Vector with Kaplan-Meier weights.
Luis Meira-Machado and Marta Sestelo
E. Kaplan and P. Meier. Nonparametric estimation from incomplete observations. Journal of the American Statistical Association, 53:457-481, 1958.
PKMW
obj <- with(colonCS, survCS(time1, event1, Stime, event)) kmw <- KMW(time = obj$Stime, status = obj$event) require(survival) colon.surv <- survfit(Surv(Stime, event) ~ 1, obj) times <- summary(colon.surv)$time surv <- summary(colon.surv)$surv nevent <- summary(colon.surv)$n.event p <- match(obj$Stime, times) kmw2 <- -diff(c(1, surv))/nevent kmw2 <- kmw2[p]*obj$event kmw2[is.na(kmw2)] <- 0 all.equal(kmw, kmw2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.