用r语言如何做线性回归分析（用R语言做数据分析）

例子：利用四种不同配方的材料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值还是从箱线图都可以说明：四种材料生产出的元件的平均寿命无显著差异。