if (FALSE) {
library(GerminaR)
library(dplyr)
gdata <- prosopis %>% mutate(across(c(nacl, temp, rep), as.factor))
smr <- ger_summary(SeedN = "seeds", evalName = "D", data = gdata)
av <- aov(grp ~ rep + nacl*temp, smr)
mc <- ger_testcomp(aov = av
, comp = c("nacl", "temp")
)
}
Run the code above in your browser using DataLab