# Example: Irrigation Tillage (main), Variety Nitrogen (sub)
data <- expand.grid(
rep = 1:3,
irrigation = c("I1", "I2"),
tillage = c("CT", "ZT"),
variety = c("V1", "V2"),
nitrogen = c("N1", "N2", "N3")
)
data$yield <- rnorm(nrow(data), 1200, 150)
anova_spd_ab_cd(data, response = "yield",
main_factor1 = "irrigation", main_factor2 = "tillage",
sub_factor1 = "variety", sub_factor2 = "nitrogen",
replication = "rep")
Run the code above in your browser using DataLab