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)
library(dplyr)
data_long <- data %>%
  gather(Treatment, Fish_Length, -ID)

print(data_long)
library(car)
anova_result <- aov(Fish_Length ~ Treatment, data = data_long)
summary(anova_result)