Displaying 1 result from an estimated 1 matches for "mccg".
Did you mean:
mcc
2006 Nov 10
0
Can plot(augPred()) be used with objects that are not groupedData?
..., 117, 107, 85, 125, 107, 128, 84, 104,
100, 114, 117, 110, 127, 106, 92, 125, 96, 123, 101,
100, 103, 105, 132, 122, 133, 107, 99, 116, 91, 128,
113, 88, 105, 112, 130)
mcc <- data.frame(
subject = factor(rep(1:12, 4)),
age =rep(seq(30, 48, by = 6), each = 12),
score = score
)
mccG <- groupedData(score ~ age | subject, data = mcc)
class(mcc)
class(mccG)
mcc.lme <- lme(score ~ age, random = ~ 1 | subject, mcc)
mccG.lme <- lme(mccG)
summary(augPred(mcc.lme, primary = ~ mcc$age)) # works
summary(augPred(mccG.lme)) # works
plot(augPred(mcc.lme, primary = ~ mcc$age)) # do...