https://github.com/cran/emplik
Raw File
Tip revision: edc37b332ad887da14df3d32cf4bbea803e0843b authored by Mai Zhou on 26 April 2015, 00:00:00 UTC
version 1.0-1
Tip revision: edc37b3
LTRC.R
LTRC <- function(x,d,w=rep(1, length(d)),y=rep(-Inf, length(x)) ) {

#### This works for Left Truncated and Right Censored data.
#### this is actually Lynden-Bell or WJT estimator.
#### Or the Kaplan-Meier/Nelson Aalen est. for right censor only data.

temp <- Wdataclean2(x,d,w)
dd <- temp$dd
ww <- temp$weight
dd[length(dd)] <- 1
xx <- temp$value

######why not use DnR?

temp <- DnR(xx,dd,ww,y=y)

NelAal <- temp$n.event/temp$n.risk
survP <- cumprod( 1 - NelAal )
NelAal <- cumsum(NelAal)
jumps <- -diff( c(1, survP) )

list(times=xx[dd==1], survjump=jumps, surv=survP, CumHaz=NelAal)
}
back to top