data <- expand.grid(env = c("E1","E2"), rep = 1:3, A = c("A1","A2"),
B = c("B1","B2"), C = c("C1","C2"))
data$yield <- rnorm(nrow(data), 1200, 150)
anova_sspd_pooled(data, "yield", "A", "B", "C", "env", "rep")
Run the code above in your browser using DataLab