data <- data.frame(
  ID = c("F01", "F02", "F03", "F04", "F05", "F06"),
  Control = c(47, 65, 44, 37, 62, 51),
  Treat1 = c(68, 55, 62, 70, 59, 63),
  Treat2 = c(78, 76, 72, 81, 76, 83),
  Treat3 = c(85, 65, 81, 98, 92, 89)
)
library(tidyr)
data_long <- gather(data, key = "Treatment", value = "Fish_Length", -ID)
library(car)
anova_result <- aov(Fish_Length ~ Treatment, data = data_long)
summary(anova_result)