data <- data.frame(
rep = rep(1:4, each = 5),
treatment = rep(c("T1", "T2", "T3", "T4", "T5"), 4),
yield = c(rnorm(5, 1200, 50), rnorm(5, 1250, 55),
rnorm(5, 1180, 45), rnorm(5, 1270, 60))
)
anova_rbd(data, response = "yield", treatment = "treatment",
block = "rep", posthoc = "duncan")
Run the code above in your browser using DataLab