JointModel (version 1.0)

pred.jplm.cumhaz: Predict the baseline cumulative hazard function at any given time point

Description

This function calculates a predicted baseline cumulative hazard function, evaluated at given time points.

Usage

pred.jplm.cumhaz(object, at=NULL, CI=FALSE)

Arguments

object

a Joint Model fit object, i.e., the result of jplm.

at

a vector of fixed time points to be evaluated.

CI

logical value; if TRUE, a 95% pointwise confidence interval is returned.

Value

If CI=FALSE, it returns a numeric vector of predicted cumulative hazard values at at= . If CI=TRUE, it returns a numeric vector of predicted cumulative hazard values, the standard error estimate of the predicted value, and its lower and upper 95% pointwise confidence interval.

References

Kim, S., Zeng, D., Taylor, J.M.G. (2016) Joint partially linear model for longitudinal data with informative drop-outs. Under revision 0, 000-000.

See Also

jplm

Examples

#	a simulated data set of longitudinal responses
attach(prostate)
#	a simulated data set of drop-out process (or, time-to-event)
attach(dropout)

#	joint fit of a partially linear model and a proportional odds model
#	with a subject-specific random intercept and random slope 
fit1 <- jplm(logPSA.postRT ~ logPSA.base + (1 + VisitTime|ID), 
                 nlm.par=prostate$VisitTime, data.y=prostate, 
                 Surv(DropTime, Status) ~ logPSA.base2, 
                 formula.frailty= ~ 1 + DropTime, 
                 id.vec=dropout$ID2, transf.par=1, data.surv=dropout)

# Evaluate at 20,...,80 percent of the maximum observed survival time 	
pts <- c(0.2, 0.4, 0.6, 0.8)*max(dropout$DropTime)  
pred.jplm.cumhaz(fit1, at=pts)
out <- pred.jplm.cumhaz(fit1, at=pts, CI=TRUE)
out$Value