# NOT RUN {
data(jansen.apple)
dat <- jansen.apple
require(lattice)
dotplot(gen ~ y/n|block, data=dat, group=inoculum, auto.key=TRUE,
main="jansen.apple",
xlab="Proportion infected per block/inoculum", ylab="Genotype")
if(require(lme4)){
# Tentative model. Needs improvement.
m1 <- glmer(cbind(y,n-y) ~ gen + factor(inoculum) + (1|block),
data=dat, family=binomial)
summary(m1)
}
# }
Run the code above in your browser using DataLab