https://github.com/cran/emplik
Raw File
Tip revision: 1a00a3dce7a5ac3854ad4978598dbc8ab597c49f authored by Mai Zhou on 23 October 2010, 00:00:00 UTC
version 0.9-6
Tip revision: 1a00a3d
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