library(readxl)
library(dplyr)
library(tidyr)
library(ScottKnott)

df <- read_excel("File_Name")

long_df <- df %>%
  pivot_longer(cols = starts_with("Rep"),
               names_to = "Rep",
               values_to = "Value")

anova_model <- aov(Value ~ Genotype, data = long_df)
summary(anova_model)

sk <- SK(anova_model)
print(sk)

groups <- data.frame(
  Genotype = sk$groups$Genotype,
  SK_Group = sk$groups$groups
)
print(groups)