if (FALSE) {
library(GerminaR)
library(dplyr)
smr <- ger_summary(SeedN = "seeds"
, evalName = "D"
, data = prosopis) %>%
mutate(across(rep:temp, as.factor))
av <- aov(grp ~ nacl*temp, smr)
anova(av)
mc <- ger_testcomp(aov = av
, comp = c("nacl", "temp"))
plotdt <- mc$table
plot <- fplot(data = plotdt
, type = "bar"
, x = "temp"
, y = "grp"
, group = "nacl"
, sig = "sig"
#, error = "ste"
, color = T
, ylimits = c(0, 120, 20)
)
plot
}
Run the code above in your browser using DataLab