AvgDailyGain {SASmixed} | R Documentation |
The AvgDailyGain
data frame has 32 rows and 6 columns.
This data frame contains the following columns:
0
< 10
< 20
< 30
indicating the amount of medicated feed additive added to
the base ration.
Treatment
as a numeric variable
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 5.3).
library(SASmixed) data(AvgDailyGain) options( contrasts = c(unordered = "contr.SAS", ordered = "contr.poly")) ## plot of adg versus Treatment by Block coplot(adg ~ as.numeric(Treatment) | Block, data = AvgDailyGain, show = FALSE) fm1Adg <- lme( adg ~ InitWt * Treatment - 1, data = AvgDailyGain, random = ~ 1 | Block) summary( fm1Adg ) # compare with output 5.1, p. 178 VarCorr( fm1Adg ) # show the variance component estimates anova( fm1Adg ) # checking significance of terms fm2Adg <- update( fm1Adg, adg ~ InitWt + Treatment ) # common slope model summary( fm2Adg ) VarCorr( fm2Adg ) anova( fm2Adg ) summary( update( fm1Adg, adg ~ InitWt + Treatment - 1 ) )