agridat (version 1.12)

gotway.hessianfly: Hessian fly damage to wheat varities

Description

Hessian fly damage to wheat varities

Format

block

block factor, 4 levels

genotype factor, 16 wheat varieties

lat

latitude, numeric

long

longitude, numeric

y

number of damaged plants

n

number of total plants

Details

The response is binomial.

References

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

Examples

# NOT RUN {
## Could the spaMM package be useful here...????

data(gotway.hessianfly)
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)
if(require(lme4)){
## l2 <- glmer(cbind(y, n-y) ~ gen + (1|block), data=dat, family=binomial,
##    control=glmerControl(check.nlev.gtr.1="ignore"))
## coef(l2)
}

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

# }
# NOT RUN {
# }