现代预防医学 2009 年第 36 卷第 2 期 Modern Preventive Medicine, 2009, Vol.36, NO.2现代预防医学 2009 年第 36 卷第 2 期 Modern Preventive Medicine, 2009, Vol.36, NO.2
� 1 ����������
��
��n�×10���
����� ��� � ��
�
����
� ����������� ������������
���
� ����������� ������������
���� � �������� �!������χ2
文章编号: 1003-8507(2009)02-0217-03 中图分类号: R195.1 文献标识码: A 【流行病与统计方法】
SPSS软件进行随机区组设计非参数检验的多重比较
米术斌, 张雷, 段一娜, 王英杰
摘要: [ 目的 ] 随机区组设计多组样本资料非参数秩和检验后, 在 SPSS 软件下实现各组样本的多重比较。 [ 方
法 ] 采用秩转换技术结合随机区组设计的方差分析进行多重比较。 [ 结果 ] 对同一数据资料进行处理, 所得结论相
同。 [ 结论 ] SPSS 软件可实现非参数检验的多重比较。
关键词: 随机区组设计; 多重比较; 秩和检验; SPSS
MULTIPLE COMPARISON OF RANDOMIZED BLOCK DESIGN OF NON-PARAMETRIC TEST
BY SPSS MI Shu -bin, ZHANG Lei, DUAN Yi -na, et al. ( Department of Preventive Medicine,
Chengde Medical Collage, Chengde 067000, China)
Abstract: [Objective] To perform multiple comparisons for multiple groups with randomized block design after nonpara-
metric test by SPSS software. [Methods] Multiple comparisons were used by changing order combination with the variance
analysis of randomized block design. [Results] We got the same conclusion by performing analysis for the same data. [Con-
clusion] Multiple comparisons can be realized with SPSS software.
Key words: Randomized block Design; Multiple comparison; Nonparametric test; SPSS
秩和检验是医学科研工作中较为常用的一种非参数统计方
法, 该方法对资料分布没有特殊要求, 可用于计量资料、 开口
资料和等级资料等的统计分析。 目前国内外权威的统计软件
如: SPSS、 SAS 等统计软件都介绍了秩和检验的方法, 比如推
断随机区组设计多个样本所来自的总体分布位置是否相同, 用
Friedman M 检验, 当 P﹤α (α 检验水准, 常为 0.05) 即 H0 假
设不成立, 接受 H1, 认为多个总体分布位置不同或不全相同。
若要进一步推断是哪两两总体分布位置不同, 可用 q 检验 [1],
但是常用的 SAS、 SPSS 等权威软件均未提供该功能 [2]。 因此,
本文以计量资料为例, 推荐随机区组设计多样本多重比较的方法。
1 资料与方法
为研究不同频率声音刺激的反应率, 随机抽取 8 名受试对
象在相同实验条件下分别接受 4 种不同频率声音的刺激, 反应
率 (%) 资料见表 1。 问 4 种频率声音刺激的反应率是否有差
别[1]?
� 1 ����� 4���
��
����
�×10���n = 8�
�� ��
�
�� ��
�
�� ��
�� ��
�
�
����
� ����
� ����
� ����
�
�� � ���� � �� � ���� � �� � ���� � ���� ����� � ����
�� ����� � �� ����� � �� ����� � ���� ����� � ����
�� � ���� � �� � ���� � �� ����� � ���� � ���� � ����
�� � ���� � �� � ���� � �� � ���� � ���� ����� � ����
�� � ���� � �� � ���� � �� � ���� � ���� � ���� � ����
�� � ���� � �� � ���� � �� � ���� � ���� ����� � ����
�� � ���� � �� � ���� � �� ����� � ���� ����� � ����
�� � ���� � �� � ���� � �� � ���� � ���� ����� � ����
Ri �� ��� �� ��� �� ����� �� �����
�
1.2 SPSS 软件方法
1.2.1 建立数据库 打开 SPSS Data Edition 窗口, 点击 Vari-
able View 标签, 输入变量名及定义变量的类型、 宽度、 小数
位数等。 其格式见表 2。
作者简介: 米术斌 (1976-), 男, 研究生学历, 讲师
作者单位: 承德医学院预防医学教研室, 承德, 067000
输入完毕后再点击 Data View 标签, 录入数据。 见表 3。
1.2.2 数据分析 操作步骤 Analyze→Nonparametric Tests→K
Related Samples…, Test Variables: 框选入: 频率 a、 频率 b、
频率 c、 频率 d, Test Type 中默认选择 Friedman 检验, 然后点
击 OK 即可。 SPSS 分析结果: 由 Friedman Test 检验知, Chi-
square 卡方统计 χ2 =15.152, P = 0.002, 表明 4 种频率声音刺
激的反应率总体分布位置不全相同。 此时, 研究者更想知道具
217· ·
现代预防医学 2009 年第 36 卷第 2 期 Modern Preventive Medicine, 2009, Vol.36, NO.2现代预防医学 2009 年第 36 卷第 2 期 Modern Preventive Medicine, 2009, Vol.36, NO.2
体是哪两种频率声音刺激的反应率不相同, 需进行多重比较,
可用多个相关样本两两比较的 q 检验[1], 但多重比较还有一定
争议, 包括 SAS、 SPSS在内的权威统计软件均未提供该功能[2]。
本文推荐一种经秩转换后用 SPSS 软件实现多重比较的新方法,
如下。
1.2.3 两两比较 采用秩转换技术结合随机区组设计的方差分
析, 该方法利用样本求出的秩次代替原变量进行参数分析, 并
进行多组样本的两两比较方法, 完成分析。 步骤 l, 秩转换:
对表 1 的资料按每一区组数据由小到大分别编秩作为新变量
秩变量, 建立数据库。 Variable View 窗口如表 4。
步骤 2, 采用方差分析: Analyze→ General Linear Model→
univariate…, 在 Univariate 对话框中, 将左侧待选变量 “秩变
� 5 ������
��� ����� ����
�� ��� ����� ����
�� ��� ����� ����
��
�� �� �
���� ���� �� �� � ��
�� ���� �� �� � ����� ����
�� �� � ����� ���� �� �� � ��
�� ���� �� �� ������ ����
�� �� � ��
�� ���� �� �� �
���� ���� �� �� �
���� ����
�� �� ������ ���� �� �� � ����� ���� �� �� � ����� ����
�� �� ������ ���� �� �� ������ ���� �� �� ������ ����
�� �� ������ ���� �� �� �
���� ���� �� �� ������ ����
�� �� ���
�� ���� �� �� �
���� ����
� �� � ��
�� ����
�� �� ������ ���� �� �� �
���� ����
� �� �
���� ����
�� �� � ����� ���� �� �� �
���� ����
� �� �
���� ����
�� �� � ����� ���� �� �� �
���� ����
� �� ���
�� ����
�� �� ������ ���� �� �� � ��
�� ���� � � � �
�
� 8 SNK�������
�������
�� � n
� ��
��
�
������ �
��
� ������� �
��
� � �������
��
� � ���
���
�� �� � ��
�
�� ���
���
�
� 3 ������
�� �� �� �� �� �� �� ��
�
�� � �
�� � � ��� �� ���
��
�� �� ��� �� ��� �� ���
�
�� � � ��� ��
�� � � ���
� ��� � � ��� � � ��� �� ���
� ��� � � ��� � �
�� � �
��
�
�� � � ��� � �
�� ��
��
� ��� � � ��� ��
�� ��
��
� ��� � � ��� � � ��� �� ���
�
� 2 �������
���� ����� �����
�� ���
� ����
� ���� ������� ������� �����
�� �� ������ � !� "� � �#$�� �#$�� !� %�&’(� ) �*��
�� +� ������ � !� "� � �#$�� �#$�� !� %�&’(� ) �*��
�� � ������ � !� "� � �#$�� �#$�� !� %�&’(� ) �*��
�� ,� ������ � !� "� � �#$�� �#$�� !� %�&’(� ) �*��
�
� �� �������� �
���� ����� �����
�� ���
� ����
� ���� ������� ������� �����
��� �������� � !� � �"#�� �"#�� � $�%&’� (�)*��
+,-�� �������� � !� � ./0+,1� �"#�� � $�%&’� (�)*��
23,� �������� � /� � �"#�� �"#�� � $�%&’� (�)*��
4��� �������� � /� � �"#�� �"#�� � $�%&’� (�)*��
� 6 �������
��
����� �����
� ��
� ���
� �� �� �� ��
������ � ���������� ��� � � ������ � � �� �!� ����!�
"#� �������� � �� �������� �������� ������
$%&������ � � ������ � ’� � � ������ � � ������ ������
$%&������ � ������� � �� � � ������ � �����’� ������
( � � ��� !�� ��� � � ��!��� � �
)*� ���� ��� ��� � � �
���)*� � ��� ��� ��� � � �
�
� 7 LSD�������
���� ���
�� � ��
� ������
�� � ��
� ������
�� � �� �� ������
�� � ��
� ������
�� � ��
� ������
��
��
� ������
�
218· ·
现代预防医学 2009 年第 36 卷第 2 期 Modern Preventive Medicine, 2009, Vol.36, NO.2现代预防医学 2009 年第 36 卷第 2 期 Modern Preventive Medicine, 2009, Vol.36, NO.2
(上接第216 页)
cm 为食管起始处, 进镜 15~25 cm 为上段, 26~34 cm 为中段,
35~40 cm 处为下段。 本组病例所见, 中段食管癌所占的比例
高, 4 762 例, 占总数的 54.78%, 其次为下段, 共 3 239 例,
占 37.26%, 上段最少, 692 例, 占 7.96%, 食管癌发病部位以
中段居多, 下段次之, 上段最少, 与文献报道相 1 年来达 50%
左右[5], 本组资料 (食管鳞癌 93.47%) 及国内相关文献显示 [6],
我省食管癌仍以鳞癌为主。 食管腺癌与食管鳞癌的好发部位差
异有统计学意义 (χ2 = 295.08, P﹤0. 005), 食管腺癌好发于
食管下段 (71.02%), 而食管鳞癌好发于食管中上段 (65.08
%)。 食管腺癌好发于食管下段的原因, 与其危险因素 Barrett
食管好发于食管下段有关[7,8]。 另有文献报道 [9], 胃食管反流与
食管腺癌有显著的因果关系, 而胃食管反流病变主要发生在食
管下段。
食管癌的发生有明显的地域性分布特点, 不同地区和人群
的食管癌发病水平差别较大。 甘肃省地处西北, 东西狭长, 南
北跨度较大, 共有 12 个地级市、 2 个自治州, 人口众多, 各
地自然生态环境不同, 人口分布不均, 社会经济发展、 人民生
活习惯存在较大差别。 统计结果显示食管癌检出率河西地区明
显高于陇南 、 陇东地区 , 其中以武威地区最高 , 检出率
5.03%。 据文献报道, 武威属我国食管癌中发区, 年平均检出
率为 5.87%, 最高达 8.13%[10]。 当地经济尚落后, 人群喜食腌
菜、 嗜烈性酒, 长期饮用硝酸盐氮及亚硝酸盐氮显著高于地下
水的地表水, 可能是食管癌高发的主要外因, 而当地居民的遗
传易感性则是食管癌发生的内在条件[11]。
本课题所调查的 163 家医院为全省截至 2006 年 12 月所有
具有胃镜检查诊断能力的医院, 覆盖了全省常住人口, 故可真
实地反映甘肃省近 30 年来食管癌的检出情况和食管癌的基本
发病趋势。 结果显示, 30 年来食管癌的检出率基本呈缓慢的
上升趋势, 食管癌仍是威胁我省群众生命的恶性肿瘤, 食管癌
的防治仍应列入我省肿瘤重点防治目标中。 在高发区如武威市
积极开展控制食管癌危险因素的工作及大力推进食管癌一、 二
级预防意义重大。
参考文献:
姚汉青, 王贵齐. Barrett 食管和食管腺癌: 东西方的差异[J]. 胃
肠病学, 2006, 11 (9): 513-515.
李建成, 杨莉芬. 甘肃省恶性肿瘤死亡率变化趋势 [J]. 中国肿
瘤, 2000, 9 (4): 170.
崔国亭, 陈占勋. 消化性溃疡与上消化道肿瘤临床流行病学分析
[J]. 中国消化内镜杂志, 1996, 13 (2): 300.
陈灏珠. 实用内科学[M]. 北京: 人民卫生出版社, 2004. 1721.
贺菊香, 耿排力, 赵久达. 733 例食管癌流行病学分析[J]. 中国
公共卫生, 2005, 21 (1): 12.
李晋革, 于飞. 胃镜诊断食管癌贲门癌 676 例分析[J]. 中华消化
志, 1999, 16 (4): 236-237.
Charles J, Lightdale MD. Esophageal Cancer [J]. The American
Journal of Gastroenterology, 1999, 94 (1): 20-29.
Van Eyken P. Definition of Barrett’ s oesophagus [J]. Acta Gas-
troenterol Belg, 2000, 63 (1): 10-12.
Lagergren J, Bergstrom R, Lindgren A, et al. Symptomatic Gas-
troesophageal Reflux as a Pisk Factor for Esophageal Adenocarcino-
ma [J]. The New England Journal of Medicine, 1999, 340 (11):
825-831.
吴正奇, 张志镒, 米登海, 等. 武威 881 例食管癌分析[J]. 卫生
职业教育, 2004, 24 (2): 128-129.
李玉民, 石斌 . 甘肃省武威市胃癌危险因素病例对照研究[J].
兰州大学学报 (自然科学版), 2003, 5: 386.
(收稿日期: 2008-02-25)
[ 1 ]
[ 2 ]
[ 3 ]
[ 4 ]
[ 5 ]
[ 6 ]
[ 7 ]
[ 8 ]
[ 9 ]
[10 ]
[11 ]
量” 选入右侧 Dependent Variable: 框中, 将 “区组” 和 “频
率分组 ” 选入 Fixed Factor [S] 框中 , 在右侧标签栏中选择
“Model…”, 在出现的对话框中选择 “Custom” , 在 “Build
TermS (S)” 下拉列表框中选择 “Main effects”, 然后将左侧
Factor & Covariates: 待选变量框中的 “区组 (F)” 和 “频率分
组 (F)” 选入右侧 Model 模型框中, 点击 Continue。 然后选择
Post Hoc… 标签, 在出现的对话框中将 Factor (s) 框中的 “频
率分组” 选入 Post Hoc Tests for: 框中, 然后根据比较的需要
选择两两比较的方法, 如本例选择 LSD 和 SNK 法, 点击 Con-
tinue, 再点击 OK。
2 结果
不同声音频率的 F 统计量为 F=11.987, P = 0.000 (见表
6), 即拒绝 H0, 接受 H1, 认为多个总体分布位置不同或不全
相同 (与 K Related Samples…检验结果一致), 经 LSD、 SNK
法分析 (见表 7、 8), 可认为频率 A、 B 和频率 C、 D 声音刺
激的反应频率有差别; 尚不能认为频率 A 和频率 B、 频率 C
和频率 D 声音刺激反应率有差别。
以上结果与用多个相关样本两两比较的 q 检验所得结果一
致[1]。
3 讨论
参数检验后进行多重比较的方法一般统计软件中均有提
供, 而秩和检验的多重比较目前常用统计软件 SPSS、 SAS 中
均未提供 [2], 只能通过公式用计算器结合查表的方法进行检
验, 但此方法计算较为复杂, 还必需有 q 界值表, 针对此种情
况, 利用非参数检验特点、 秩转换技术和随机区组设计资料方
差分析多重比较方法完全可通过 SPSS 软件实现两两比较。 该
方法利用电脑和统计软件进行统计分析, 结论更准确和可靠,
方法更简便和快捷, 值得向广大科学研究工作者推广和使用。
参考文献:
[ 1 ] 孙振球, 徐勇勇 . 医学统计学 [M]. 北京 : 人民卫生出版社 ,
2004. 133-136.
[ 2 ] 张文彤. SPSS11 统计分析教程 [M]. 北京: 北京希望电子出版
社, 2002. 256-257.
(收稿日期: 2008-02-28)
� 8 SNK�������
�������
�� � n
� ��
��
�
������ �
��
� ������� �
��
� � �������
��
� � ���
���
�� �� � ��
�
�� ���
���
�
� 7 LSD�������
���� ���
�� � ��
� ������
�� � ��
� ������
�� � �� �� ������
�� � ��
� ������
�� � ��
� ������
��
��
� ������
�
219· ·
本文档为【SPSS软件进行随机区组设计非参数检验的多重比较】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。