# Data from 3 years
data <- data.frame(
year = rep(c("Y1", "Y2", "Y3"), each = 20),
rep = rep(rep(1:4, each = 5), 3),
treatment = rep(c("T1", "T2", "T3", "T4", "T5"), 12),
yield = rnorm(60, 1200, 150)
)
anova_rbd_pooled(data, response = "yield", treatment = "treatment",
environment = "year", block = "rep")
Run the code above in your browser using DataLab