1、
multcomp包-aov（）
agricolae包-SNK.test（）
library(multcomp)
RACE <- as.character(RACE)
RACE <- as.factor(RACE)
Height <- as.numeric(Height)
fit1 <- aov(Height ~ RACE,data=data4_1)
summary(fit1)
qqPlot(lm(Height~RACE,data=data4_2),simulate=TRUE,main="Q-Q Plot",labels=FALSE)
bartlett.test(Height~RACE,data=data4_2)
library(multcomp)
RACE <- as.character(RACE)
RACE <- as.factor(RACE)
Height <- as.numeric(Height)
fit2 <- aov(Height ~ RACE,data=data4_2)
summary(fit2)
install.packages("agricolae")
library(agricolae)
out <- SNK.test(fit2,"RACE",console = TRUE)
print(SNK.test(fit2,"RACE",group = FALSE))
