Diseño en parcelas divididas (split plot desing)

library(MASS)
library(collapsibleTree)
library(ggplot2)
library(lattice)
data("oats")
avena=oats

Arbol colapsable

collapsibleTree(avena,
                hierarchy = c("B","V",
                              "N","Y"))
bwplot(avena$Y~avena$N|avena$V)

book.url <- "http://stat.ethz.ch/~meier/teaching/book-anova"
john <- readRDS(url(file.path(book.url, "data/john.rds")))

bwplot(john$mass~john$fertilizer|john$variety)

View(john)