Description Usage Arguments Value Author(s) References See Also Examples
View source: R/pred.jplm.nonlinear.R
This function calculates a predicted nonlinear effect function evaluated at given points.
1 | pred.jplm.nonlinear(object, nlm.par, at=NULL, CI=FALSE)
|
object |
a Joint Model fit object, i.e., the result of |
nlm.par |
a vector of nonlinear effect covariate, as specified in |
at |
a vector of fixed points to be evaluated. |
CI |
logical value; if |
If CI=FALSE
, it returns a numeric vector of predicted nonlinear effect at at=
, the standard error estimate of the predicted value, and test result based on the asymptotic normality. If CI=TRUE
, it returns a numeric vector of predicted nonlinear effect, 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 time-to-event (e.g., drop-out process)
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 measurement time
pts <- c(0.2, 0.4, 0.6, 0.8)*max(prostate$VisitTime)
pred.jplm.nonlinear(fit1, prostate$VisitTime, at=pts)
out <- pred.jplm.nonlinear(fit1, prostate$VisitTime, at=pts, CI=TRUE)
out$Value
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.