agridat (version 1.8.1)

gotway.hessianfly: Hessian fly damage to wheat varities

Description

Hessian fly damage to wheat varities

source

C. A. Gotway and W. W. Stroup. A Generalized Linear Model Approach to Spatial Data Analysis and Prediction Journal of Agricultural, Biological, and Environmental Statistics, 2, 157-178.

Details

The response is binomial.

References

The GLIMMIX procedure. http://www.ats.ucla.edu/stat/SAS/glimmix.pdf.

Examples

dat <- gotway.hessianfly
dat$prop <- dat$y / dat$n
desplot(prop~long*lat, dat, out1=block, text=gen, cex=1, shorten='no',
        main="gotway.hessianfly")


# Block random.  See Glimmix manual, output 1.18.
# Note: (Different parameterization)
require(lme4)
l2 <- glmer(cbind(y, n-y) ~ gen + (1|block), data=dat, family=binomial,
   control=glmerControl(check.nlev.gtr.1="ignore"))
coef(l2)

require('asreml')
a2 <- asreml(prop ~ gen, data=dat, random=~block,
             family=asreml.binomial(),
             weights=n)
coef(a2)