agridat (version 1.8.1)

gomez.stripsplitplot: Rice strip-split-plot experiment

Description

A strip-split-plot experiment with three reps, genotype as the horizontal strip, nitrogen fertilizer as the vertical strip, and planting method as the subplot factor.

source

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

Examples

dat <- gomez.stripsplitplot

# Layout
desplot(gen~x+y, dat, out1=rep, out2=gen, col=nitro, text=planting,
cex=1)

# Gomez table 4.19, ANOVA of strip-split-plot design
dat <- transform(dat, nf=factor(nitro))
m1 <- aov(yield ~ nf * gen * planting +
          Error(rep + rep:nf + rep:gen + rep:nf:gen), data=dat)
summary(m1)

# There is a noticeable linear trend along the y coordinate which may be
# an artifact that blocking will remove, or may need to be modeled.
# Note the outside values in the high-nitro boxplot.
require("HH")
interaction2wt(yield ~ nitro + gen + planting + y, dat, x.relation="free")