getsurvprob: Get survprob from pec packages

In whcsu/sfse: Random Survival Forest with Space Extensions (short line)

Description

Get survprob from pec packages

Usage

getsurvprob(obj, times)

Arguments

obj	survfit object
times

Details

Get survprob from pec packages

Author(s)

pec

References

from pec packages

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (obj, times) 
{
    if (!inherits(obj, "survfit")) 
        stop("obj is not of class survfit")
    class(obj) <- NULL
    lt <- length(times)
    nsurv <- times
    if (length(times) == length(obj$time)) {
        if (all(times == obj$time)) 
            return(obj$surv)
    }
    inside <- times %in% obj$time
    for (i in (1:lt)) {
        if (inside[i]) 
            nsurv[i] <- obj$surv[obj$time == times[i]]
        else {
            less <- obj$time[obj$time < times[i]]
            if (length(less) == 0) 
                nsurv[i] <- 1
            else nsurv[i] <- obj$surv[obj$time == max(less)]
        }
    }
    nsurv
  }

whcsu/sfse documentation built on Dec. 4, 2019, 2:11 p.m.