data <- read.delim("clipboard") data$Treat <- as.factor(data$Treat) nom <- bartlett.test(data$value~data$Treat,data = data) nom oneway<-aov(data$value~data$Treat,data = data) anova(oneway) library(agricolae) out <- LSD.test(oneway,"data$Treat",p.adj="none") out