Displaying 1 result from an estimated 1 matches for "team9".
Did you mean:
team
2007 Apr 27
1
Example of mcmcsamp() failing with lmer() output
...<<1>>=
gunload <- read.table(hh('datasets/gunload.dat'), header = T)
gunload$method <- factor(gunload$method, labels = c('new', 'old'))
gunload$physique <- factor(gunload$group, labels = c('slight',
'average', 'heavy'))
gunload$team9 <- factor(rep(1:9, each = 2))
@
This is a split plot design with randomly selected teams for each
physique, which was crossed with two methods.
Following the Oats example in chap. 1 of MEMSS, and adapting to lmer:
<<2>>=
require(lme4)
gl.lmer <- lmer(rounds ~ method * physiqu...