bac<-read.table("c:\\users\\86167\\Desktop\\练习六表.txt",header=T)

boxplot(day~type,data=bac,col="red") 

ba.an<-aov(lm(day~type,data=bac)) 
summary(ba.an)

##             Df Sum Sq Mean Sq F value  Pr(>F)   
## type         2  70.43   35.21   6.903 0.00379 **
## Residuals   27 137.74    5.10                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1