data(egsingle)
str(egsingle)
fm1 <- lme(math~year*size+female, egsingle,
random = ~1|schoolid/childid,
control = list(EMv = TRUE, msV = TRUE, opt = 'optim'))
summary(fm1)
Run the code above in your browser using DataCamp Workspace