agridat (version 1.5)

gomez.splitsplit: Rice split-split plot experiment

Description

Grain yield of three varieties of rice grown in a split-split plot arrangement with 3 reps, nitrogen level as the main plot, management practice as the sub-plot, and rice variety as the sub-sub plot.

Usage

data(gomez.splitsplit)

source

Gomez, K.A. and Gomez, A.A.. 1984, Statistical Procedures for Agricultural Research. Wiley-Interscience. Page 143. Used with permission of Kwanchai Gomez.

Examples

dat <- gomez.splitsplit
dat$nf <- factor(dat$nitro)

desplot(yield~x*y, dat)
desplot(nf~x*y, dat, out1=rep, col=management, num=gen, cex=1)

if(require(HH)){
position(dat$nf) <- c(0,50,80,110,140)
interaction2wt(yield~rep+nf+management+gen, data=dat,
               relation=list(x="free", y="same"),
               rot=c(90,0), xlab="",
               par.strip.text.input=list(cex=.7))
}

# AOV.  Gomez page 144-153
m0 <- aov(yield~ nf * management * gen + Error(rep/nf/management),
         data=dat)
summary(m0) # Compare to Gomez, p. 153.