方差分析公式.docx

方差分析公式.docx

《方差分析公式.docx》由会员分享，可在线阅读，更多相关《方差分析公式.docx（11页珍藏版）》请在冰点文库上搜索。

方差分析公式

方差分析公式

（2012-06-2611:

03:

09）

标签：

分类：

统计方法

杂谈

方差分析

方差分析（analysisofvariance,简写为ANOV或ANOVA）可用于两个或两个以上样本均数的比较。

应用时要求各样本是相互独立的随机样本；各样本来自正态分布总体且各总体方差相等。

方差分析的基本思想是按实验设计和分析目的把全部观察值之间的总变异分为两部分或更多部分，然后再作分析。

常用的设计有完全随机设计和随机区组设计的多个样本均数的1：

匕较。

—、完全随机设计的多个样本均数的1：

匕较

又称单因素方差分析。

把总变异分解为组间（处理间）变异和组变异（误差）两部分。

目的是推断k个样本所分别代表的pl,m2,……pk是否相等，以便比较多个处理的差别有无统计学意义。

其计算公式见表19-6.

表19-6完全随机设计的多个样本均数比较的方差分析公式

变异来源

藹均差平方和SS

目由度V

均方MS

F

总

ZX2-C*

N-1

组间（处理组间）

k-1

SStgiDi/v型间

MSt§iPi/MS型间

」—C

b

组（误差）

SS总-SS朋

N-k

SSj§/v组

*C=（XX）2/N=Xni;k为处理组数

表19-7F值.P值与统计结论

a

F值

P值

统计结论

0.05

>0.05

不拒绝Hor差别无统计学意义

0.05

2Fo.05（vl.V2）

<0.05

拒绝Ho，接受Hi,差别有统计学意义

0.01

»Fo.01（vl.V2）

<0.01

拒绝Ho，接受H-差别有肓度统计学意义

方差分析计算的统计量为Fz按表19-7所示关系作判断。

例19.9某湖水不同季节氯化物含量测量值如表19-8,问不同季节氯化物含量

有无差别？

表19-8某湖水不同季节氯化物含量（mg/L）

冷

春

夏

秋

冬

22.6

19.1

18.9

19.0

22.8

22.8

13.6

16.9

21.0

24.5

17.2

17.6

16.9

18.0

15.1

14.8

20.0

15.2

16.6

13.1

21.9

18.4

14.2

16.9

21.5

20.1

16.7

16.2

21.2

21.2

19.6

14.8

ZXij

■J

167.9

159.3

131.9

129.3

588.4（IX）

ns

8

8

8

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

8

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

F

P

141.170

3

47.057

9.38

<0.01

组

140.465

28

5.017

总

281.635

31

二、随机区组（配伍组）设计的多个样本均数1：

匕较

又称两因素方差分析。

把总变异分解为处理间变异、区组间变异及误差三部分。

除推断k个样本所代表的总体均数,pl,p2,......pk是否相等外,还要推断b个区组所代表的总体均数是否相等。

也就是说,除比较多个处理的差别有无统计学意义外，还要比较区组间的差别有无统计学意义。

该设计考虑了个体变异对处理的影响,故可提高检验效率。

表19-10随机区组设计的多个样本均数比较的方差分析公式

变异来源

离均差平方和SS

目由

度V

均

方

MS

F

总

ZX2-C

N-1

处

理

间

j）_1C

k-1

SS

fcts

/v处

理

MS

/M

Sig

差

区组间

$£弓（弓丿c（44.84『十一=1.84）+（M.€

b3

-K9X.456CS=2.90418

b-1

ss

区组

/v区

组

MS

区组

MS

误

差

SS总・SSya-SS区型

V总

s

g-v

ss逞差

/v«

s

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

0

45

90

135

1

5.27

5.27

4.94

4.61

20.09

2

5.27

5.22

4.88

4.66

20.03

3

5.88

5.83

5.38

5.00

22.09

4

5.44

5.38

5.27

5.00

21.09

5

5.66

5.44

5.38

4.88

21.36

6

6.22

6.22

5.61

5.22

23.27

7

5.83

5.72

5.38

4.88

21.81

8

5.27

5.11

5.00

4.44

19.82

工Xij

•J

44.84

44.19

41.84

38.69

169.56（ZX）

Nj

8

8

8

8

32（N）

Xi

5.6050

5.5238

5.2300

4.8363

•

J

252.1996

245.0671

219.2962

187.5585

904.1214（IX2）

Hl:

四个总体均数不等或不全相等

c（二0.05

区组间：

H0:

八个区组的总体均数相等，即口二血二……p8

Hl:

八个区组的总体均数不等或不全相等

c（二0.05

先作表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=

b

（44•84严+（41•19尸一"】・84）・+（38,69?

8

—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资料的方差分析表

变异来源

SS

V

MS

F

P

处理间

2.90438

3

0.9681

77.448

<0.01

区组间

2.49800

7

03569

28.552

<0.01

0.26297

21

0.0125

总

5.66538

31

三、多个样本均数间的两两比较的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资料作两两比较

H0:

任两个季节的湖水氯化物含量的总体均数相等，即pA二pB

Hl:

任两总体均数不等，即pA/pB

表19-13|q|值.P值与统计结论

a

|q|

P值

统计结论

0.05

>0.05

不拒绝Ho■差^椒计学息义

0.05

>qo.O5（v.a）

<0.05

拒绝Ho.接受Hi,差别有统计学意义

0.01

nqo.Ol（v.a）

<0.01

拒绝Ho,接受H-差别有髙度统计学意义

a=0.051•将四个样本的均数由大到小排列编秩，注明处理组。

Xi

167.9

159.3

131.9

129.3

处理组

春

夏

秋

餐

秩次

1

2

3

4

2.计算sxA-xB本例各处理组的样本含量nl相等，按式（19,14）计算两均数

差的标准误。

已知MS组二5.017,n二8

5.017

=0.7f）!

0

3.列两两比较的q检验计

算表（表19-14）

表19-14两两比较的q检验计算表

A与B

Xa-Xb

（2）

M.a

q值

qo.os（v^）

QO.01（v.a）

P值

（1）

（3）

（4）=

（2）/0.7919

（5）

（6）

（7）

⑴与（4）

38.6

4

4&744

3.85

4.80

<0.01

⑴与（3）

36.0

3

45.460

3.49

4.45

<0.01

⑴与

（2）

8.6

2

10.860

2.89

3.89

<0.01

（2）与（4）

30.0

3

37.884

3.49

4.45

<0.01

（2）与（3）

27.4

2

34.600

2.89

3.89

<0.01

（2）与（4）

2.6

2

3.283

2.89

3.89

<0.05

表中第

（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值应该是一样的。