This function estimates the survival probability at tx when a piecewise
exponential distribution is fitted to (times,cens) cens = 0 for censored, cens = 1 for uncensored.
the change point is tchange and lamest is the estimated parameters
Usage
pexeest(times, cens, tchange, tx)
Arguments
times
All the event/censoring times used to fit the model
cens
censoring status used to fit the model
tchange
Change-points
tx
Time points to estimate the survival probability
Value
quan survival probability
lamest Lambda estimates for time periods divided by the change-points
# NOT RUN {data(pexeest_times_censoring)
data(t100)
times = pexeest_times_censoring[,1]
cens = pexeest_times_censoring[,2]
pexeest(times, cens, 28.03014, t100)
# }