data <- expand.grid(
rep = 1:4,
irrigation = c("I1", "I2", "I3"),
tillage = c("CT", "MT", "ZT")
)
data$yield <- rnorm(nrow(data), 1200, 150)
anova_strip(data, response = "yield", horizontal_factor = "irrigation",
vertical_factor = "tillage", replication = "rep")
Run the code above in your browser using DataLab