R做Tukey檢驗

mydata <- read.table(a.txt, header=FALSE, sep="\t")
a <- c(a[,1]);


#a <- c(1.092859111,0.916497508,0.526556309,0.402578703,0.372798667,0.52515351,0.413846575,0.564489973)
b=factor(c(rep(1,2),rep(2,2),rep(3,2),rep(4,2)))
c <- data.frame(a,b)
fit1 <- aov(a ~ b)
summary(fit1)
TukeyHSD(fit1)