https://github.com/cran/emplik
Tip revision: 23907dd7c7efb4547f0300a60c2b026df465087d authored by Mai Zhou on 10 September 2004, 00:00:00 UTC
version 0.8-2
version 0.8-2
Tip revision: 23907dd
WKM.R
WKM <- function(x,d,zc=rep(1,length(d)),w=rep(1,length(d))) {
temp <- Wdataclean3(x,d,zc,w)
dd <- temp$dd
ww <- temp$weight
dd[length(dd)] <- 1
######why not use DnR?
allrisk <- rev(cumsum(rev(ww)))
survP <- cumprod( 1 - (dd*ww)/allrisk )
jumps <- -diff( c(1, survP) )
logel <- sum(ww[dd==1]*log(jumps[dd==1])) + sum(ww[dd==0]*log(survP[dd==0]))
list(times=temp$value, jump=jumps, surv=survP, logel=logel)
}
################################################################
##This function compute a weighted Kaplan-Meier estimator
## x = times, d = censoring status, w = weights
################################################################