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)