# MinutePlot validation -- One-way ANOVA + Tukey letters (PlantGrowth)
# Install once (skip on later runs):
#   install.packages(c("emmeans", "multcomp", "multcompView"))   # exactly what THIS script uses

# Load packages (every session):
library(emmeans); library(multcomp); library(multcompView)

# Versions (compare against those stated on this page):
cat(R.version.string, "\n")
for (p in c("emmeans", "multcomp", "multcompView")) cat(p, as.character(packageVersion(p)), "\n")

# Load data:
d <- PlantGrowth                 # built-in; or: d <- read.csv("plantgrowth_r.csv")  (same folder as script)

# Analysis (one print per table shown on the page, in page order):
fit <- aov(weight ~ group, data = d)
print(summary(fit))
em <- emmeans(fit, ~ group)
print(summary(em))
print(multcomp::cld(em, adjust = "tukey", Letters = letters, alpha = 0.05))