joineRML (version 0.2.1)

plot.ranef.mjoint: Plot a

Description

Displays a plot of the BLUPs and approximate 95% prediction interval for each subject.

Usage

# S3 method for ranef.mjoint
plot(x, ...)

Arguments

x

an object inheriting from class ranef.mjoint, representing the estimated random effects for the mjoint object from which it was produced.

...

additional arguments; currently none are used.

Value

an object inheriting from class ggplot, which displays a trellis plot with a separate panel for each effect, showing a dotplot (with optional error bars indicating approximate 95% prediction intervals if the argument postVar = TRUE is set in the call to ranef) for each subject (by row).

References

Pinheiro JC, Bates DM. Mixed-Effects Models in S and S-PLUS. New York: Springer Verlag; 2000.

See Also

ranef.mjoint.

Examples

## Not run: ------------------------------------
# require(ggplot2)
# data(heart.valve)
# hvd <- heart.valve[!is.na(heart.valve$log.grad) & !is.na(heart.valve$log.lvmi), ]
# set.seed(1)
# 
# fit1 <- mjoint(formLongFixed = log.lvmi ~ time,
#     formLongRandom = ~ time | num,
#     formSurv = Surv(fuyrs, status) ~ 1,
#     data = hvd,
#     timeVar = "time")
# 
# plot(ranef(fit1, postVar = TRUE))
## ---------------------------------------------