现代统计图形 - 科学网—博客

5.19 因素效应图 109
1 > plot.design(warpbreaks, col = "blue")
mean of breaks
25 30 35
A
B
Factors
L
M
H
wool tension
图 5.26: 经纱断裂数据的因素效应图:每种羊毛(A、B)和拉力强度(L、
M、H)下断裂数目的均值。
x为包含自变量(分类变量)的数据框,它也可以包含因变量,这种情
况下第二个参数就不必提供了;y为因变量;fun为计算因变量水平的函数;
data为所有数据,包含自变量和因变量,当参数x为一个公式时,则会用到
本参数提取数值,否则不必提供本参数;其它参数用于调整图形外观或设
置图形标题、标记等,此处略去。
我们以经纱断裂数据warpbreaks为例,这个数据包含三个变量:经
纱断裂数目breaks、羊毛种类wool(A、B两种)和拉力强度tension(L、
M、H三种强度)。我们可以通过因素效应图看羊毛种类和拉力强度两个变
量分别对经纱断裂根数的影响,如图5.26。图中的刻度线展示出了不同羊毛
种类和拉力强度下断裂根数的均值。我们也可以通过方差分析模型进一步
看这两种因素的影响:
1 > summary(aov(breaks ~ wool + tension, data = warpbreaks))
Df Sum Sq Mean Sq F value Pr(>F)
wool 1 451 451 3.34 0.0736 .
tension 2 2034 1017 7.54 0.0014 **
Residuals 50 6748 135
---
Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1