agridat: gomez.stripsplitplot – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

gomez.stripsplitplot

Strip-split-plot experiment of rice

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.

Used with permission of Kwanchai Gomez.

Source

Gomez, K.A. and Gomez, A.A.. 1984, Statistical Procedures for Agricultural Research. Wiley-Interscience. Page 155.

Examples

## Not run: 

library(agridat)
data(gomez.stripsplitplot)
dat <- gomez.stripsplitplot

# Layout
libs(desplot)
desplot(dat, gen ~ col*row,
        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)

# 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.
libs("HH")
interaction2wt(yield ~ nitro + gen + planting + row, dat,
               x.between=0, y.between=0,
               x.relation="free")


## End(Not run)

agridat

Agricultural Datasets

v1.18

CC BY-SA 4.0

Authors

Kevin Wright [aut, cre] (<https://orcid.org/0000-0002-0617-8673>)

Initial release