data <- data.frame(
treatment = rep(c("T1", "T2", "T3", "T4"), each = 5),
yield = c(rnorm(5, 1200, 50), rnorm(5, 1350, 60),
rnorm(5, 1100, 55), rnorm(5, 1450, 65))
)
anova_crd(data, response = "yield", treatment = "treatment", posthoc = "all")
Run the code above in your browser using DataLab