gomez.splitsplit {agridat} | R Documentation

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

### Format

A data frame with 135 observations on the following 7 variables.

`rep`

  block, 3 levels

`nitro`

  nitrogen fertilizer, in kilograms/hectare

`management`

  plot management

`gen`

  genotype/variety of rice

`yield`

  yield

`col`

  column position in the field

`row`

  row position in the field

### Source

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

Used with permission of Kwanchai Gomez.

### References

H. P. Piepho, R. N. Edmondson. (2018). A tutorial on the statistical analysis of factorial experiments with qualitative and quantitative treatment factor levels. Jour Agronomy and Crop Science, 8, 1-27. https://doi.org/10.1111/jac.12267

### Examples

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

    if(require(desplot)){
      desplot(nf ~ col*row, dat,
              # aspect unknown
              out1=rep, col=management, num=gen, cex=1,
              main="gomez.splitsplit")
      desplot(yield ~ col*row, dat,
              # aspect unknown
              out1=rep, main="gomez.splitsplit")
    }

    if(require(HH)){
    position(dat$nf) <- c(0,50,80,110,140)
    interaction2wt(yield~rep+nf+management+gen, data=dat,
                   main="gomez.splitsplit",
                   x.between=0, y.between=0,
                   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) # Similar to Gomez, p. 153.

[Package agridat version 1.16 Index]