agridat (version 1.16)

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.

Format

yield

grain yield in kg/ha

planting

planting factor, P1=broadcast, P2=transplanted

rep

rep, 3 levels

nitro

nitrogen fertilizer, kg/ha

gen

genotype, G1 to G6

col

column

row

row

Details

Note, this is a superset of the the 'gomez.stripplot' data.

Examples

# NOT RUN {
data(gomez.stripsplitplot)
dat <- gomez.stripsplitplot

# Layout
if(require(desplot)){
  desplot(gen ~ col*row, dat,
          out1=rep, col=nitro, text=planting, cex=1,
          main="gomez.stripsplitplot")
}


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

# }
# NOT RUN {
# 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 + row, dat,
                 x.between=0, y.between=0,
                 x.relation="free")
# }