rat.growth {SMPracticals} | R Documentation |
Data on the weights of 30 rats each week for 5 weeks.
data(rat.growth)
A data frame with 150 observations on the following 3 variables.
a factor with levels 1-30
takes values 0-4
rat weight (units unspecified)
Gelfand, A. E., Hills, S. E., Racine-Poon, A. and Smith, A. F. M. (1990) Illustration of Bayesian inference in normal data models using Gibbs sampling. Journal of the American Statistical Association, 85, 972–985.
Davison, A. C. (2003) Statistical Models. Cambridge University Press. Page 460.
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))