agridat (version 1.5)

gotway.hessianfly: Hessian fly damage to wheat varities

Description

Hessian fly damage to wheat varities

Usage

data(gotway.hessianfly)

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')

# Block random.  See Glimmix manual, output 1.18.
# Note: (Different parameterization)
require('lme4')
l2 <- lmer(cbind(y, n-y) ~ gen + (1|block), data=dat, family=binomial)
coef(l2)

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

# To do: Add spatial.  Throws an error in asreml