方差分析公式文档格式.docx
《方差分析公式文档格式.docx》由会员分享,可在线阅读,更多相关《方差分析公式文档格式.docx(11页珍藏版)》请在冰点文库上搜索。
N-k
SSj§
/v组
*C=(XX)2/N=Xni;
k为处理组数
表19-7F值.P值与统计结论
a
F值
P值
统计结论
0.05
<
Fo.O5(v1.V2)
>
不拒绝Hor差别无统计学意义
2Fo.05(vl.V2)
拒绝Ho,接受Hi,差别有统计学意义
0.01
»
Fo.01(vl.V2)
拒绝Ho,接受H-差别有肓度统计学意义
方差分析计算的统计量为Fz按表19-7所示关系作判断。
例19.9某湖水不同季节氯化物含量测量值如表19-8,问不同季节氯化物含量
有无差别?
表19-8某湖水不同季节氯化物含量(mg/L)
冷
春
夏
秋
冬
22.6
19.1
18.9
19.0
22.8
13.6
16.9
21.0
24.5
17.2
17.6
18.0
15.1
14.8
20.0
15.2
16.6
13.1
21.9
18.4
14.2
21.5
20.1
16.7
16.2
21.2
19.6
ZXij
■J
167.9
159.3
131.9
129.3
588.4(IX)
ns
8
32(N)
Xi
20.99
19.91
16.49
16.16
3548.51
3231.95
2206.27
2114.11
11100.84(IX2)
HO:
湖水四个季节氯化物含量的总体均数相等,即==
Hl:
四个总体均数不等或不全相等
c(二0.05
先作表19-8下半部分的基础计算。
C二(8)2/N二(588.4)2/32二10819.205
SS总二》X2-C二11100.84-10819.205=281.635
S5加f«
]
-:
b?
iy.2C5=141.170
V总二N-1二31
(167.9)24(I関.3)B+(13].9j8+(l25-3)a
V组间二k-1二4-1二3
SS组二SS总-SS组间二281.635-141.107=140.465
V组二N-k二32-4二28
MS组间二SS组间/v组间二141.107/3二47.057
MS组二SS组/v组二140.465/28二5.017
F二MS组间/MS组二47.057/5.017二9.380
以vl(即组间自由度)二3,v2(即组自由度)二28查附表19-2,F界值表,
得F0.05(3,28)=2.95,F0.01(3,28)=4.57.$例算得的F二9.380>
F0.01
(3,28),P<
0.01,按g0.05检验水准拒绝H0,接受H1,可认为湖水不
同季节的氯化物含量不等或不全相等。
必要时可进一步和两两比较的q检验,
以确定是否任两总体均数间不等。
资料分析时,常把上述计算结果列入方差分析表,如表19-9.
表19-9例19.9资料的方差分析表
SS
V
MS
P
141.170
3
47.057
9.38
组
140.465
28
5.017
281.635
31
二、随机区组(配伍组)设计的多个样本均数1:
又称两因素方差分析。
把总变异分解为处理间变异、区组间变异及误差三部分。
除推断k个样本所代表的总体均数,pl,p2,......pk是否相等外,还要推断b个区组所代表的总体均数是否相等。
也就是说,除比较多个处理的差别有无统计学意义外,还要比较区组间的差别有无统计学意义。
该设计考虑了个体变异对处理的影响,故可提高检验效率。
表19-10随机区组设计的多个样本均数比较的方差分析公式
变异来源
离均差平方和SS
目由
度V
均
方
ZX2-C
处
理
间
j)_1C
fcts
/v处
/M
Sig
差
区组间
$£
弓(弓丿c(44.84『十一=1.84)+(M.€
b3
-K9X.456CS=2.90418
b-1
ss
区组
/v区
误
SS总・SSya-SS区型
V总
s
g-v
ss逞差
/v«
C、k、N的意义同表19-6zb为区组数
例19.10为研究酵解作用对血糖浓度的影响,从8名健康人中抽血并制成血滤液。
每个受试者的血滤液被分成4份,再随机地把4份血滤液分别放置0,45,
90,135分钟,测定其血:
唐浓度(表19-11),试问放置不同时间的血糖浓度有无差别?
处理间:
H0:
四个不同时间血糖浓度的总体均数相等,即口二卩2二二p4
表19-11血滤放置不同时间的血糖浓度(mmol/L)
区组号
放置时间(分)
受试者小计
IXij
•J
45
90
135
1
5.27
4.94
4.61
20.09
2
5.22
4.88
4.66
20.03
5.88
5.83
5.38
5.00
22.09
4
5.44
21.09
5
5.66
21.36
6
6.22
5.61
23.27
7
5.72
21.81
5.11
4.44
19.82
工Xij
44.84
44.19
41.84
38.69
169.56(ZX)
Nj
5.6050
5.5238
5.2300
4.8363
•
J
252.1996
245.0671
219.2962
187.5585
904.1214(IX2)
区组间:
八个区组的总体均数相等,即口二血二……p8
八个区组的总体均数不等或不全相等
先作表19-11下半部分和右侧一栏的基本计算。
C=(XX)2/N二(169.56)2/32二898.45605
SS总二》X2-C二904.1214-898.45605二5.66535
V总二N-1二32-1二31
SS处理==:
>
_(7=
(44•84严+(41•19尸一"
】・84)・+(38,69?
—898・456砧=2-00438
V处理=k-l=4-l=3
+(19,82户
SS嗣.(20.09)^(20,03)-^
k4
45635=2.49EC.»
V区组二b-1二8-1二7
SS误差二SS总-SS处理-SS区组二5.66535-2.90438-2.49800二0.26297
V误差二(k-1)(b-1)=3x7=21
MS处理二SS处理/v处理二2.90438/3二0.9681
MS区组二SS区组/v区组二2.49800/7二0.3569
MS误差二SS误差/v误差二0.26297/21二0.0125
F处理二MS处理/MS误差二0.9681/0.0125=77.448
F区组二MS区组/MS误差二0.3569/0.0125=28.552
推断处理间的差别,按vl二3,v2=21查F界值表,得F0.005(3,21)=3.07,F0.01(3,21)=4.87,P<
0.01;
推断区组间的差别,按vl=7fv2=21查F界值表,得F0.05(7,21)=2.49,F0.01(7,21)=3.64,P<
0.01按c(二0.05检验水准皆拒绝H0,接受H1,可认为放置时间长短会影响血糖浓度且不同受试者的血糖浓度亦有差别。
但尚不能认为任两个不同放置时间的血糖浓度总体均数皆有差别,必要时可进一步作两两t出交的q检验。
表19-12例19.10资料的方差分析表
处理间
2.90438
0.9681
77.448
区组间
2.49800
03569
28.552
0.26297
21
0.0125
5.66538
三、多个样本均数间的两两比较的q检验
经方差分析后,若按c(二0.05检验水准不拒绝H0,通常就不再作进一步分析;
若按c(二0.05甚至c(二0.01检验水准拒绝H0,且需了解任两个总体均数间是否都存在差别,可进一步作多个样本均数间的两两比较。
两两比较的方法较多,在此仅介绍较常用的q检验(Newman-Keuls法)
(各组ni相等)公式(19.14)
“V2\rzn3)
(各组ni不等)公式(19.15)
式中,xA-xB为两两对比中,任两个对比组A、B的样本均数之差;
sxA-xB为两样本均数差的标准误;
ni为各处理组的样本含量;
nA,nB分别为A、B两对比组的样本含量;
MS误差为单因素方差分析中的组均方(MS组)或两因素方差分析中的误差均方(MS误差)。
计算的统计量为q,按表19-13所示关系作判断。
例19.11对例19.9资料作两两比较
任两个季节的湖水氯化物含量的总体均数相等,即pA二pB
任两总体均数不等,即pA/pB
表19-13|q|值.P值与统计结论
|q|
QO.05(v.a)
不拒绝Ho■差^椒计学息义
qo.O5(v.a)
拒绝Ho.接受Hi,差别有统计学意义
nqo.Ol(v.a)
拒绝Ho,接受H-差别有髙度统计学意义
a=0.051•将四个样本的均数由大到小排列编秩,注明处理组。
处理组
餐
秩次
2.计算sxA-xB本例各处理组的样本含量nl相等,按式(19,14)计算两均数
差的标准误。
已知MS组二5.017,n二8
=0.7f)!
3.列两两比较的q检验计
算表(表19-14)
表19-14两两比较的q检验计算表
A与B
Xa-Xb
(2)
M.a
q值
qo.os(v^)
QO.01(v.a)
(1)
(3)
(4)=
(2)/0.7919
(5)
(6)
(7)
⑴与(4)
38.6
4&
744
3.85
4.80
⑴与(3)
36.0
45.460
3.49
4.45
⑴与
(2)
8.6
10.860
2.89
3.89
(2)与(4)
30.0
37.884
(2)与(3)
27.4
34.600
2.6
3.283
表中第
(1)栏为各对比组,如第一行1与4,指A为第1组,B为第4组。
第
(2)栏为两对比组均数之差,如第一行为XI与X4之差,余类推。
第(3)栏为四个样本均数按大小排列时,A、B两对比组围所包含的组数a,如第一"
1与4"
围包含4个组,故a二4.第(4)栏是按式(19.13)计算的统计量q值,式中的分母0.7919是按式(19.14)计算出来的SXA-XBM(5)、(6)栏是根据误差自由度v与组数a查附表19-3q界值表所得的q界值本例v误差二28,因q界值表中自由度一栏无28,可用近似值30或用插法得出q界值,本例用近似值30查表f当a二4时,q0.05(30,4)=3.85,qO.Ol(30f4)=4.80f余类推。
第(7)栏是按表19-13判定的。
4.结论由表19-14可见,除秋季与冬季为P<
0.05外,其它任两对比组皆为P
0.01,按g0.05检验水准均拒绝H0,接受H1,可认为不同季节的湖水氯化物含量皆不同,春季氯化物含量最高,冬季含量最低。
PS:
逬行方差分析前必须要做方差齐性检验和正态分布检验,至于如何做,方法很多了/常见的正态性检验有Kolmogorov-Smirnov检验和Shapiro-Wilk检验。
方差齐性常采用Bartlett检验。
相同的数据,不同的软件,采用相同的方法给岀的p值应该是一样的。