if (FALSE) {
library(agridat)
data(jansen.apple)
dat <- jansen.apple
libs(lattice)
xyplot(inoculum ~ y/n|gen, data=dat, group=block,
layout=c(3,1),
main="jansen.apple",
xlab="Proportion infected per block/inoculum",
ylab="Inoculum level")
## libs(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