- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
平有显著性差异.
33
二、有交互作用的双 因素试验方差分析
34
在无交互作用时,对因素 A , B 各水
平的每种组合只进行一次试验,即
t 1 . 当 要 考虑 因 素 间的 交 互作 用 A B 时,在各水平组合下需要做重复试 验.设每种水平组合下试验次数均为 t
t 1 .此时相对应的数学模型就是前述
j 1
0,
i 1
ij
j 1
0.
下面我们分两种情况来讨论双因素试验方差分析.
13
一、无交互作用的双 因素试验方差分析
14
如果因素 A 与因素 B 之间不存在交互作用,则
ij 0 , i 1, 2, , r ; j 1, 2, , s ,
于是
ij i j
即每种水平组合 Ai , Bj 下的总体平均值 ij 可以看成是总
25
为清楚起见,将上述分析结果汇总成下表,称下表为无交互作用的双因素试验方差分析表.
表 无交互作用的双因素试验方差分析表
方差来源 平方和
自由度
均方
F值
临界值
显著性
因素 A
SA
r 1
SA
SA r 1
F r 1, r 1s 1
FA
因素 B
SB
s 1
SB
SB s 1
F s 1, r 1s 1
FB
误差
SE
Xijk ~ N ij , 2 i 1, 2, , r ; j 1, 2, , s ,
可将上式改写成如下形式:
Xij ij ijk i 1, 2, , r ; j 1, 2, , s ; k 1, 2, , t
ijk ~ N 0, 2 各ijk相互独立
7
为进行统计分析,需将均值 ij 作适当分解,为此,令
2 ij
T2 rs
1135.42 109.82 43
130.75 ,
30
SA
1 s
r
Ti
2
i 1
T2 rs
1 3207 .74 3
1 12
109 .82
64.5767
,
SB
1 r
s
Tj2
j 1
T2 rs
1 4261.64 1 109.82
4
12
60.74 ,
SE ST SA SB 130.75 64.5767 60.74 5.4333 ,
2,
L,
s
r
s
i
1
i
0,
j
j 1
0.
上式就是无交互作用的双因素试验方差分析的数学模型16.
由上式可知,为判断 A 对试验指标的影响是否显
著,即等价于检验假设
H0 A : 1 2 r 0 . 类似地,判断因素 B 对试验指标的影响是否显著,
即等价于检验假设
H0B : 1 2 s 0 .
1 rs
r i 1
s
ij ,
j 1
i
1 s
s
ij
j 1
,
j
1 r
r i 1
ij
8
i i , i 1, 2, , r j j , j 1, 2, , s
ij ij i j ,
其中 称为总平均, i 称为因素 A 的第 i 个水平 Ai
的效应, j 称为因素 B 的第 j 个水平 B j 的效应.
r 1s 1;
24
为此,选取 FA , FB 分别作为检验假设 H 0 A , H 0B 的统计 量.按照假设检验的程序,对显著性水平 ,确定临界值
F r 1, r 1s 1 ,F s 1, r 1s 1 .当 FA F r 1, r 1s 1 时,拒绝 H0A , FB F s 1, r 1s 1时,拒绝 H0B ,
① SE ~ 2r 1s 1;
2
②
当
H
0
A
为真时,
S
A 2
~ 2
r 1
,而且 S A 与 SE 相互独立,从而
SA
FA
r 1 SE
~ Fr 1,
r 1s 1
r 1s 1 ;
③
当
H
0
B
为真时,
SB
2
~ 2
s 1
,而且 SB 与 SE 相互独立,从而
SB
FB
r 1 SE
~ Fs 1,
r 1s 1
31
S A 自由度: r 1 3 ,
SB 自由度: s 1 2 ,
SE 自由度: r 1s 1 6,
于是,
SA 64.5767
FA SE 3
3 5.4333
23.7707 ,
6
6
SB 64.5767
FB SE 2
2 5.4333
33.5376 .
6
6
32
表 试验数据方差分析表
方差来源
i 1 j 1
r
SA s
Xi X 2
i 1
s
SB r
X j X 2
j 1
21
r
由于 Xi 是水平 Ai 下的所有观察值的平均,所以 X i X 2
i 1
反映了
X1, X2, , Xr 之间的差异程度.这种差异是由于因素 A 的不同水平所引起的, 因此 SA 称为因素 A 的效应平方和,简称为因素 A 的平方和.同样
将 ST 分解为
r s
ST
X ij X 2
i1 j 1
r s
Xij Xi X j X Xi X X j X 2
i1 j 1
r
s
X ij X i X j X 2 s r
Xi X 2 r s
Xj X 2
i1 j 1
i 1
j 1
r s
17
为构造检验统计量,我们仿造单因素试验方差分析 的做法,记
X
1 rs
r i 1
s
X ij ,
j 1
X i
1 s
s
X ij
j 1
,
X
j
1 r
r i 1
X ij
,
18
r s
ST
X ij X 2
i 1 j 1
其中 ST 称为总离差平方和,简称为总平方和,
也称为总变差平方和.
19
T
2 j
676
1281.64
2304
4261.64
29
解:
H0 A :试验温度的各个水平无显著性差异;
H0B :铜含量的各个水平无显著性差异.
rs
rs
X
2 ij
1135
.42
,T
Xij 109 .8 , r 4 , s 3 .
i1 j 1
i1 j 1
所以, ST
r i 1
s j 1
X
2
设在某项试验中有两个因素 A , B 在变化.因素 A 有 r
个不同的水平
A1, A2, , Ar , 因素 B 有 s 个不同水平
B1, B2, , Bs .
在水平组合 Ai , Bj 下的试验结果用 X ij 表示.
3
我们假定
X ij i 1, 2, , r ; j 1, 2, , s
r 1s 1
SE
r
SE
1s
1
总和
ST
rs 1
26
SA
SB
其中 FA
r 1 SE
SA SE
, FB
s 1 SE
SB . SE
r 1s 1
r 1s 1
为计算方便,常采用下列公式计算各偏差平方和.
ST
r i1
s j 1
X
2 ij
T2 rs
S A
1 s
r
Ti
2
i 1
T2 rs
SB
1 r
s
Tj2
X 221, X 222, , X 22t
X 2s1, X 2s2 , , X 2st
Ar
X r11, X r12, , X r1t
X r21, X r22, , X r2t
X rs1, X rs2 , , X rst
6
由于 X ijk k 1, 2, L , t 是取自总体 X ij 中的样本,则有
看作是取自正态总体 Xij ~ N ij , 2 中的容
量为 t 的样本.将这些数据列成下表
5
B 因素 各水平 B1
A 因素 各水平
A1
X111, X112, , X11t
A2
X 211, X 212, , X 21t
B2
Bs
X121, X122, , X12t
X1s1, X1s2 , , X1st
rs
2
Xij Xi X j X Xi X 2
X ij X i X j X X j X
i1 j 1
i1 j 1
r s
2
Xi X X j X .
i1 j 1
20
可以证明,上述平方和分解中交叉项均为 0.所以
其中
ST SE SA SB ,
r s
SE
X ij X i X j X 2
的道理, SB 称为因素 B 的效应平方和,简称为因素 B 的平方和.
22
又由于
SE ST S A SB 这表明 S E 是从总离差平方和 ST 中扣 除因素 A ,B 的效应平方和 S A 和 S B
之后的残量,这一残量反映了随机误
差因素的影响,因此 SE 称为误差平
方和.
23
与单因素试验的方差分析的讨论相类似,可以证明以下结论: