例子:利用四种不同配方的材料A1、A2、A3、A4生产出来的元件,测得其使用寿命如下图所示。问:四种不同配方下元件的使用寿命有无显著差异?
R语言中的aov()函数提供了方差分析的计算,aov()函数的使用方式如下:
aov(formula,data=NULL)
其中formula是方差分析的公式,data是数据框。另外,可使用summary()列出方差分析表的详细信息。
接下来用数据框的格式输入数据,调用aov()函数计算方差分析,用summary()提取方差分析的信息。
> lamp<-data.frame(
X=c(1600, 1610, 1650, 1680, 1700, 1700, 1780, 1500, 1640, 1400, 1700, 1750, 1640, 1550, 1600, 1620, 1640, 1600, 1740, 1800, 1510, 1520, 1530, 1570, 1640, 1600),
A=factor(c(rep(1,7),rep(2,5),rep(3,8),rep(4,6)))
)
> lamp.aov<-aov(X~A,data = lamp)
> summary(lamp.aov)
Df Sum Sq Mean Sq F value Pr(>F)
A 3 49212 16404 2.166 0.121
绘制箱线图
>plot(lamp$X~lamp$A)
无论是从p值还是从箱线图都可以说明:四种材料生产出的元件的平均寿命无显著差异。
,