Displaying 1 result from an estimated 1 matches for "mix1".
Did you mean:
mix
2010 Jul 02
1
xyplot: key inside the plot region / lme: confidence bands for predicted
....)
2: With lme(), how can I generate confidence bands or prediction
intervals around the fitted values? Once
I get them, I'd like to add them to the plot.
Example:
library(nlme)
library(lattice)
Ortho <- Orthodont
Ortho$year <- Ortho$year - 8 # make intercept = initial status
Ortho.mix1 <- lme(distance ~ year * Sex, data=Ortho,
random = ~ 1 + year | Subject, method="ML",
# correlation = corAR1 (form = ~ 1 | Subject)
)
Ortho.mix1
# predicted values
grid <- expand.grid(year=0:6, Sex=c("Male", "Female"))
grid$age <- grid$...