sampleSurv: Samples fitted survival function

In icenReg: Regression Models for Interval Censored Data

Description

Samples fitted survival function

Usage

sampleSurv(fit, newdata = NULL, p = NULL, q = NULL, samples = 100)

Arguments

fit	Either an ic_bayes or ic_par fit
newdata	A data.frame with a single row of covariates
p	A set of survival probabilities to sample corresponding time for
q	A set of times to sample corresponding cumulative probability for
samples	Number of samples to draw

Details

For Bayesian models, draws samples from the survival distribution with a given set of covariates. Does this by first drawing a set of parameters (both regression and baseline) from fit$samples and then computing the quantiles of the distribution (if p is provided) or the CDF at q.

If a ic_par model is provided, the procedure is the same, but the sampled parameters are drawn using the normal approximation.

Not compatible with ic_np or ic_sp objects.

Author(s)

Clifford Anderson-Bergman

Examples

data("IR_diabetes")
fit <- ic_par(cbind(left, right) ~ gender, data = IR_diabetes)

newdata <- data.frame(gender = "male")
time_samps <- sampleSurv(fit, newdata, 
                         p = c(0.5, .9), 
                         samples = 100)
# 100 samples of the median and 90th percentile for males                        

prob_samps <- sampleSurv(fit, newdata, 
                         q = c(10, 20),
                         samples = 100)
# 100 samples of the cumulative probability at t = 10 and 20 for males                        

Example output

Loading required package: survival
Loading required package: Rcpp
Loading required package: coda

icenReg documentation built on Oct. 23, 2020, 8:11 p.m.