powered by
Pixel
data(Pixel) Pixel$DS = getGroups(Pixel, ~ 1 | Dog/Side, level = 2) fm1 <- lme(pixel ~ day + I(day^2), data = Pixel, random = list(DS = ~ 1, Dog = ~ 1)) summary(fm1) VarCorr(fm1)
Run the code above in your browser using DataLab