Displaying 1 result from an estimated 1 matches for "grpdat".
Did you mean:
grpdata
2011 Jul 07
2
How do I overlay two trellis plots of lme fitted lines produced by plot.augPred?
..., 28L, 29L, 30L, 31L,
32L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 77L, 78L,
79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L, 91L,
92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L,
104L, 105L, 106L, 107L, 108L, 109L), class = "data.frame")
library(nlme)
grpdat = groupedData(y ~ age | ptno,
data = tmpdat,
FUN = mean,
labels = list(x="Age", y="Y"))
fit.lme = lme(grpdat)
plot(augPred( fit.lme ), ylim=c(20,65))
plot(augPred(update(fit.lme, y ~...