SPSS多组率间两两比较

:多个样本率的卡方检验及两两仳

伴们都掌握了吗如果不止两个分组,又该如何进行卡方检验以及之后的两两比较呢来

某医生拟探讨药物以外的其他方法是否可降低患者的胆固醇浓度,如增强体育锻炼、

减少体重及改善饮食习惯等

位高胆固醇、生活习惯差的受试者,

组给予降胆固醇药物一组给予飲食干预,另一组给予运动干预经过

该医生重新测量受试者的胆固醇浓度,分为高和正常两类

该医生收集了受试者接受的干预方法(

)和试验结束时胆固醇的风险程度

)等变量信息,并按照分类汇总整理部分数据如下:

我们使用非参数检验之多个独立樣本的秩和检验的例子:为调查A、B、C三种治疗措施对患者谷丙转氨酶(ALT)的影响某科室将45名患者随机等分为三组,每组15人分别采取A、B、C三种治疗措施。血常规测得治疗后ALT水平(U/L)如下试问应用三种治疗措施后,患者的ALT水平是否有差异

表1. 三组患者治疗后的ALT水平(U/L)

在仩一篇文章()中,Kruskal-Wallis检验结果为:H=34.155P<0.001,按照α=0.05的检验标准拒绝H0,即三组的ALT水平有差别但是Kruskal-Wallis不能检验出三组之间的两两比较情况,在SPSS軟件的22.0版本中已经能够通过软件直接得出两两比较的结果了。以下内容演示如何在SPSS软件22.0版本中进行多组之间的两两比较

1. 将数据录入SPSS,將分组情况填写在第一列(1=A组2=B组,3=C组)ALT水平填写在第二列。

首先结果中报告了Kruskal-Wallis检验结果,P<0.001拒绝H0假设。双击此结果表格出现详细嘚两两比较检验结果如下图。

左边为Kruskal-Wallis检验结果右边为各分组之间的对比情况,在View(视图)选项中选择Pairwise Comparisons(两两比较)即可出现两两比较嘚结果。

我们把两两比较的结果调出来仔细分析报告的统计量有Test Statistic(检验统计量)、Std.Error(标准误)、Std.Text Statistic(标准检验统计量)、Sig.(P值)、Adj.Sig.(调整後P值),结果中我们看到报告了两个P值①直接得出的P值是未经过矫正的,本次计算的两两比较次数为3次(即1组和3组、1组和2组、2组和3组比較)所以使用原始的P值与显著水平0.05比较时假阳性率会升高。因此软件自动对原始P值进行矫正将P值乘以比较次数得到②调整后P值,调整後P值即可与显著水平0.05直接进行比较得出结论

Kruskal-Wallis检验结果为:H=34.155,P<0.001按照α=0.05的检验标准,拒绝H0可以认为患者经过A、B、C三种治疗措施后血常規中ALT水平有差别。经过两两比较A组与B组之间患者的ALT水平有差别(调整后P=0.002)。A组与C组之间患者的ALT水平有差别(调整后P<0.001)B组与C组之间患鍺的ALT水平无差别(调整后P=0.054)。

(更多内容可关注“医咖会”微信公众号:传播医学知识和研究进展探讨临床研究方法学。)

在日常数据处理中我们经常会遇到多个样本率比较的情况,即R*C表格样式的数据在我们进行Χ2检验后,推断的结论如果是拒绝了H0,接受了H1那么我们只能认为各总体率之間是具有差异性的,但是对于任两个总体率是否具有差异我们不得而知。此时如果运用四格表资料的Χ2检验将会增大犯Ⅰ类错误的概率因此,多个样本率的多重比较不能简单地这样分析我们已知,比较的方法有Χ2分割法(partitions ofΧ2 method)、Scheffe’可信区间法和SNK法今天我们仅介绍Χ2汾割法。

在我们分割完成后我们必须重新规定检验水准,其目的是保证检验假设中Ⅰ型错误α的概率不变,通常有两种情况:

多个实验組间的两两比较

实验组与同一个对照组的比较

接下来我们以一个例子来实际操作一下:

1、多个实验组间的两两比较制作值标签、加权与㈣格表的spssΧ2检验的过程一致,这里不再赘述我们以A组和B组的比较为例,由上述公式计算可得α’=0.0125

表可得ν=1P=0.0125时,Χ2=6.24<6.76按α’=0.0125水准,A组与B組的比较拒绝H0接受H1,即可以认为A法与B法的有效率有差别

2、各实验组与同一个对照组的比较

我们以B组为对照组,其他两组为试验组spss操莋过程与上面一致,但检验水准α’计算公式不同,由上公式可得α’=0.0125使用spss计算Χ2值,并查表可得:

按α’=0.0125的检验水准A组与B组的比绝拒絕H0,接受H1即可以认为A法与B法的有效率有差别,C组与B组不拒绝H0即还不能认为两组的有效率具有差异。

我们再来以文字的方式总结一下步驟:

1. 根据实验内容计算出相应的α’;

2. 分别将分割好的实验组与对照组的数据录入注意频数列要进行加权,并且要在变量视图中设置好對应的值标签;

3. 分析→描述统计→交叉表注意行与列的选择,见下图;

4. 在统计量里选择卡方(见下图)继续后确定,即可得出结果

我要回帖

 

随机推荐