第10章单因素方差分析

第10章 单因素方差分析

单因素方差分析(0ne-Way ANOV A)，又称一维方差分析，它能够对单因素多个独立样本 的均数进行比较，可以用10种检验方法对变量间的均数进行两两比较(即多重比较检验)并给出方差分析表，还可以作出5种类型图形(Type of plots)和2种均数图形(Means plot options)

10.1 单因素方差分析的计量资料

[例10—1] 某社区随机抽取了30名糖尿病患者、IGT 异常人和正常人进行载脂蛋白 (mg ／dL)测定，结果示于表10—1。试问3组人群的载脂蛋白测定结果含量是否相同？(倪宗瓒．卫生统计学．第4版，北京：人民卫生出版社，2001.50) 组别（B ）

载脂蛋白测定 糖尿病(1)

85.7 105.2 109.5 96.0 115.2 95.3 110.0 100.0 125.6 111.0 106.5 96.0 124.5 105.1 76.4 95.3 110.0 95.2 99.0 120.0

144.0 117.0 110.0 109.0 103.0 123.0 127.0 121.0 159.0 115.0 IGT 异常(2) 正常人(3)

本例是一个完全随机设计的单因素方差分析。已建立SAS 数据集文件并保存Sasuser.onewav4。

(1)进入SAS ／Win(v8)系统，单击Solutions －Analysis －Analyst ，得到分析家窗口。

(2)单击File-open By SAS Name —Sasuser-0neway4—0K ，调入数据文件。

(3)在“分析家”窗口单击Statistics-ANOV A-One way ANOV A ，得到图10—1所示对话框。本例因变量(Dependent)为A(载脂蛋白)，单击A —Dependent 。自变量(1ndependent)： B(3种人的组别)，单击B —Independent 。

图10.1 0ne —way ANOV A ：0neway4(单因素方差分析)对话框

(4)单击Tests 按钮，得到图10—2所示对话框。在此对话框的ANOV A(F —检验)选项 中可进行如下设置。

Analysis of variance ，方差分析。

Welch ’s variance-weighted ANOV A ，威尔奇方差—权重方差分析。

Tests for equal variance ，相等方差检验，即方差齐性检验。

Barlett ’s test ，巴特尼特检验。

Brown-Forsythe test ，布朗—福塞斯检验。

Levene ’s test ，列文检验。本例以上都选。

图10—2 ANOV A(方差分析)选项卡

(5)单击图10。2所示对话框的Power Analysis标签，得到图10—3所示的Power Analysis选项卡。本例选择Perform Power analysis(执行功效分析)，。(A1phas)值取默认值0.05

图10—3 Power Analysis(功效分析)选项卡

(6)单击0K按钮，返回图10—1所示对话框。单击Means按钮，得到图10—4所示对话框。此对话框的Comparisons选项卡的Comparison method选项区域，列出了10种待选多重均数比较方法。

Bonferroni t-test，修正最小显著差异t检验法(本例选择此项)。

Tukey’s HSD，’图基可靠显著性差异检验法。

Duncan’s multiple-range test，邓肯多重范围检验法。

Dunnett’s t·test，邓尼特t检验法。

Fisher’s LSD，费雪尔最小显著差异(LSD，least·significant—difference)法。

Gabriel’s multiple—comparison procedure，嘉百列多重比较过程。

Student—Newman·Keels multiple—rangete8t，SNK多重范围检验法。

Waller—Duncan k—ratios。test，娃尔—邓肯A—比率t检验法。

Schaffer’s multiple comparison procedure，谢弗多重比较过程。

Ryan-Einot—Gabriel—Wel6ch multiple—rangete8t，R—E—G。W多重范围检验法。

图10-4 One-way ANOV A：Means(待选多重均数比较方法)对话框

(7)单击Bonferroni t-test，得到图10。5所示对话框。在comparisons选项卡中可进行相关设置。本例显著性水平(Signillcancelevel)为0．贴。也可以选择其他值。主效应(Main effects)为B，加选(Add)效应／方法(effect／method)为B／Bonferronit。test，单击B—Add，完成相应设置。

图10—5 One—way ANOV A：Means(已选均数比较)对话框

(8)单击Breakdown标签，得到图10—6所示对话框的Breakdown选项卡。数学变量统计量全选。

图10-6 Breakdown(均数细分)选项卡

(9)单击0K按钮，返回图10—1所示对话框。单击P10ts按钮，得到图10—7所示对话框。

图10-7 one-Way ANOV A：P1ots(作图)对话框

(10)单击0K按钮，返回图10—1所示对话框。单击0K按钮，得到如下数值结果。

图形结果如图10-8到图10-10所示。

图10—8 Box plot of A by B(箱形图)