数据读取

bac<-read.table("D:\\Rdownload\\lianxi\\first\\test1.txt", header = T)

将ba数据框中的type转换为因子（factor），进行方差分析

bac$type<-as.factor(bac$type)
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

结论：

###“type”这一项的自由度为2，平方和为70.43，均方为35.21，F值为6.903，对应的P值为0.00379，标记了两个星号（**），表明该结果是高度显著的。

###残差（Residuals）的自由度为27，平方和为137.74，均方为5.10。

###表格底部还给出了显著性代码的含义：0代表p < 0.001，* * 代表 p < 0.01，*代表 p < 0.05，·代表 p < 0.1，空白则代表 p > 0.1。

###总结来说，这个结果表明“type”变量对模型中的响应变量具有显著影响。

绘制箱线图

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