Description Usage Arguments Value Author(s) References See Also Examples
View source: R/pred.jplm.cumhaz.R
This function calculates a predicted baseline cumulative hazard function, evaluated at given time points.
1 | pred.jplm.cumhaz(object, at=NULL, CI=FALSE)
|
object |
a Joint Model fit object, i.e., the result of |
at |
a vector of fixed time points to be evaluated. |
CI |
logical value; if |
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.
Sehee Kim
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.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | # 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
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.