data(rat.growth)
library(nlme)
rat.fit <- groupedData( y~poly(week,2) | rat,
data = rat.growth,
labels = list( x = "Week",
y = "Weight" ),
units = list( x = "", y = "(?)") )
summary(lme(rat.fit))
Run the code above in your browser using DataLab