swh:1:snp:7d062e61dfe39b25a891852629fe3416f36554de
Tip revision: 6d2aaee84bc7d13a97cf480fbaf9e89663bcb3cd authored by Mai Zhou on 09 January 2009, 00:00:00 UTC
version 0.9-5
version 0.9-5
Tip revision: 6d2aaee
MMRtime.R
MMRtime <- function(x, d, age) {
#### estimate Mean/Median Residual lifetime over age.
temp <- WKM( x=x , d=d )
tivec <- temp$times
pivec <- temp$jump
if( age >= tivec[length(tivec)] ) stop("age too large")
if( age < tivec[1] ) warning("age smaller than first event time")
pivec[ tivec < age ] <- 0
Sage <- sum( pivec )
fenzi <- sum( (tivec - age)*pivec )
MRtime <- fenzi/Sage
Ptheta <- Sage/2
Cprob <- cumsum(pivec)
posi <- sum(Cprob < Ptheta)
theta <- tivec[posi+1]
list(MeanResidual = MRtime, MedianResidual = theta - age)
}