rat.growth | 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))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.