data(Pixel)
fm1 <- lme(pixel ~ day + I(day^2), data = Pixel, random = ~1 | Dog/Side)
#fm1 <- lme(pixel ~ day + I(day^2), data = Pixel,
# random = list(Dog = ~ day, Side = ~ 1))
summary(fm1)
VarCorr(fm1)
Run the code above in your browser using DataCamp Workspace