配对卡方检验(Paired χ² Test)——Stata软件实现

发布于 2022年1月9日 星期日 16:47:29 浏览:3887
原创不易,转载请注明来源,感谢!
附件下载:
配对卡方检验.zip 请勿重复点击,如无响应请耐心等待或稍后再试。

在前面文章中介绍了配对χ2检验(Pairedχ2 Test)理论,本篇文章将使用实例演示在Stata软件中实现配对χ2检验(Pairedχ2 Test)的操作步骤。

关键词:Stata; 配对四格表; 配对卡方检验; McNemar检验

一、案例介绍

某实验室分别用A方法和B方法对600名患者血清中某抗体进行测定,数据见图1,其中1为阳性,2为阴性。本文案例可从“附件下载”处下载。问两种方法的检测结果有无差别?

图1

二、问题分析

本案例的分析目的是比较两种方法的检测结果有无差别,属于配对设计,可以认为是在比较两种方法检测结果的一致性,或者比较两者阳性检测率的差异。每个数字代表的是一个对子,而不是单独的一个频数。对于本案例可采用2×2配对χ2检验(McNemar test)进行分析。但需要满足3个条件:

条件1:观察变量为二分类变量,本案例数据为检测方法的阴性和阳性,满足该条件。

条件2:试验方法或干预措施为两个,本案例为A和B两种检测方法,满足该条件。

条件3:研究设计为配对设计,即对同一批样本或研究对象进行两种方法的检测或干预。本案例对同一组患者使用两种抗体检测方法,为配对设计,满足该条件。

三、软件操作及结果解读

(一) 适用条件判断

根据研究设计,本案例数据满足条件1—条件3。

(二) 统计描述及推断

1. 软件操作

①描述性分析

tabi 110 130\30 330, row column

图2

②配对χ2检验

 mcci 110 130 30 330

图3

2. 结果解读

(1) 统计描述

图2中给出了A方法和B方法检测一致性的对子数和不一致的对子数,各自检测的阳/阴性频数和百分比。可知两种检验方法一致的有440对(同为阳性110对,同为阴性330对),不一致的有160对(A方法为阳性,B方法为阴性130对;A方法阴性,B方法阳性30对)。

(2) 统计学推断

对于2×2配对四格表,定义四个格子分别为a、b、c和d,a、d为结果一致的频数,b和c为结果不一致的频数。需要根据b+c的大小,选择分析方法:当b+c≥40时,选择McNemar检验(2×2配对χ2检验);当20≤b+c<40时,选择校正McNemar检验(Stata无法实现连续校正);当b+c<20时,选择确切概率法。由于本案例b+c=160≥40,所以使用配对χ2检验的结果。

图3给出了配对χ2检验及确切概率法的统计量、自由度和P值。由于本案例b+c=160≥40,所以使用配对χ2检验的结果。可知,χ2=62.50,P<0.001,提示两种方法检出的阳性率差异有统计学意义,即两种检测方法结果不同。

四、结论

本研究采用2×2 配对χ2检验(Paired Samples McNemar Tests)评价A、B两种方法检测患者血清中某抗体阳性率的一致性。结果显示检测不一致的有160对(A方法为阳性,B方法为阴性130对;A方法阴性,B方法阳性30对);配对χ2检验分析发现,两种检测方法结果不同(χ2=62.50,P<0.001)。

End
文章目录 沉浸式阅读