查看原文
其他

SPSS详细操作:多个独立样本列联表的卡方检验

2017-01-21 李延龙 医咖会

上一期我们聊到了独立样本四格表的卡方检验,实际上临床也经常会遇到分组超过2个的计数资料,这一期我们来一起搞定多个独立样本列联表的χ2检验。

一、问题与数据

以下是胃癌真菌病因研究中3种食物样品的真菌检出率,比较3个检出率有无差异。


表1  物种食物样品的真菌检出率


二、对数据结构的分析

本例是独立四格表χ2检验的拓展,由两个分组增加到多个分组,分析思路与四格表χ2检验类同,需要注意的是,这里我们不光想知道多个分组间有无差异,如果差异存在统计学意义,那么具体到组间两两比较是否均存在差异。

三、SPSS分析方法

多个独立样本列联表χ2检验的SPSS操作与四格表一样,这里不再赘述(忘记的小伙伴赶快戳SPSS详细操作:独立样本四格表的χ2检验,复习一下)。需要注意的是,不同于四格表χ2检验,SPSS对于R*C列联表χ2检验不会自动输出Fisher确切概率检验结果,如果样本例数较少,建议在Exact设置中勾选Exact(如下图)。

 


四、结果解读

表2  统计汇总


 

表3  卡方检验结果


 

多个独立样本列联表χ2检验的结果选择:


1、所有理论频数≥5,看Pearson Chi-Square的结果;


2、超过20%的理论频数<5或至少1个理论频数<1,看Fisher’s Exact Test结果(也可以考虑增加样本量或者依据专业判断适当合并行或列,再进行χ2检验)。


本例中SPSS提示没有理论频数小于5,且最小的理论频数为8.00,故直接选择Pearson Chi-Square结果,即χ2=22.841,P<0.001,提示三种食物中真菌检出率不同。

五、组间的两两比较

通过上述χ2检验,明确了三种食物的真菌检出率并不相同,此时我们还需要进一步考虑三种食物真菌检出率到底谁与谁之间的差异存在统计学意义,这里就需要用到“卡方分割”,通俗讲就是把R*C列联表拆分成若干个四个表分别进行χ2检验,进而判断不同组两两比较差异是否用统计学意义,但是,因为多组比较可能会增加犯I类错误概率,所以还需要对χ2检验的P值进行校正,这里主要介绍 Bonferroni校正。


本例中需要进行3次两两比较,校正的检验水准α=0.05/比较次数=0.05/3=0.0167。


到这里,有的小伙伴要问了,SPSS数据库中原来有3组,怎么才能方便地构造任意两组的“四格表”,进行χ2检验呢?这里教大家一个SPSS中比较实用的小技巧——选择特定对象进行统计分析。


A、菜单的Data中找到Select Cases



 

B、Select Cases中提供了多种用于选择研究对象的方式 ,这里我们将用到条件筛选(如下图)




C、条件筛选中提供了丰富的筛选公式,假如想选择1-大米和2-地瓜粉,可以做如下图设置,“食物=1|食物=2”,这里“|”代表“或者”,即数据库只要有1或者2都会被选中进行统计分析→Continue。




 

按照上面介绍的小技巧,我们就可以进行任意两组的四格表χ2检验(表4)


表4. 不同食物真菌检出率比较 



如上表,按照校正的检验水准α=0.0167,大米和地瓜粉,大米和豆酱之间的真菌检出率差异具有统计学意义,而地瓜粉和豆酱之间差异无统计学意义。

六、撰写结论

大米、地瓜粉和豆酱的真菌检出率并不相同(χ2=22.841,P<0.001),其中地瓜粉最高为96.7%,其次为豆酱为80.0%,大米最低为43.3%。大米的真菌检出率分别与地瓜粉和豆酱相比差异均有统计学意义(Bonferroni校正,P<0.0167),而地瓜粉和豆酱之间真菌检查率差异无统计学意义(Bonferroni校正,P>0.0167)。


PS: 多个独立样本的χ2检验除了包含上述R*2列联表卡方检验外,还包含R*C卡方检验,即我们考虑的指标变量为多分类(例如血型),其统计分析思路和SPSS操作分析与R*2列联表卡方检验一致,这里不再赘述。


统计学精彩回顾

1. SPSS详细操作:独立样本四格表的卡方检验

2. SPSS实例教程:有序多分类Logistic回归


研究进展精彩回顾

1. 一文读懂:1973-2016肺癌治疗进展

2. 柳叶刀: 2016年医疗大事件盘点

医咖会微信:medieco-ykh

关注医咖会,轻松学习统计

我们建了一个微信群,有临床研究设计或统计学方面的难题?快加小咖个人微信(xys2016ykf),拉你进群,和其他小伙伴们一起交流学习吧。


点击左下角“阅读原文”,看看医咖会既往推送了哪些研究设计或统计学文章。还可以到医咖会公众号下方的自定义菜单,点击“直接搜索”,查找你想了解的内容。

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存