卡方检验

Stata第六章卡方检验

本节STATA命令摘要

[by分层变量名:]tab2变量1变量2[,allchi2exactcellcolumnrow]

tabi#11#12[...]\[#21#22[...][\...][,allchi2exactcellcolumnrow]

•列联表分析

STATA命令：

[by分层变量:]tab2变量1变量2[，

allchi2lichi2exactcellcolumnrow]

上述命令中，变量1为行计数变量；变量2为列计数变量；all表示卡方(c2)检验，似然比(likelihoodratio)检验以及一些统计描述指标和检验，但不包括Fisher精确检验；exact表示Fisher精确检验；chi2表示c2检验；lichi2表示likelihoodratio检验；cell表示输出的列联表中显示每个观察计数值占该列联表总观察计数值的比例；row表示输出的列联表中显示每个观察计数值占该观察计数值所在行的各观察计数值总数的比例；coloumn表示输出的列联表中显示每个观察计数值占该观察计数值所在的列各观察计数值总数的比例。例：某地调查肝癌病人与健康人饮用“醋冷水”(一种以冷水和醋为主要成分的饮料)的习惯。用group=1表示肝癌组患者和group=2表示健康人；用custom=1表示经常饮用醋冷水；custom=2表示偶尔饮用醋冷水和custom=3表示从不饮用醋冷水。具体资料为：(摘自医学统计方法，金丕焕主编，p163)。

组别经常偶尔从不饮用合计肝癌组26442898

健康组28491794

合计549345192

显然这是一个病例对照研究，所以每组人数是人为确定的，因此只需计算各组"经常"，"偶而"和"从不饮用"占本组的频数以及检验患肝癌是否与饮水习惯有关。

tab2groupcustom,rowchi2

->tabulationofgroupbycustom

|custom

group|123|Total

-----------+--------------------------------------------+----------

1|①264428|98

|②26.5344.9028.57|100.00

-----------+--------------------------------------------+----------

2|③284917|94

|④29.7952.1318.09|100.00

-----------+--------------------------------------------+----------

Total|⑤549345|192

|⑥28.1248.4423.44|100.00

Pearsonchi2(2)=2.9497Pr=0.229

①该行表示第一组(肝癌组)的3个观察数；②该行表示第一组的各个观察数的

占第一组观察总数的百分比；③该行表示第二组(健康组)的3个观察数；④该行表示第二组的各个观察数的占第二组观察总数的百分比；⑤该行表示关于饮用醋冷水习惯的三个分类：“经常”，“偶尔”和“从不”的合计数；⑥该行表示上述三个合计数分别占总样本数的百分比。

从上述结果可知：卡方值为2.9497以及自由度为2，p为0.229，所以根据当前资料和结果显示：患肝癌病与是否饮用醋冷水无关。

•列联表分析也可以用立即命令执行：

tabi#11#12[...]\[#21#22[...][\...][,allchi2exactcellcolumnrow]

#11表示列联表中第一行第一列的观察数；#12表示列联表中第一行第二列的观察数；\表示换一行；#21表示列联表中第二行第一列的观察数，其它以此类推。子命令：allchi2exactcellcolumnrow与上述命令完全相同。以上例的资料为

例：相应的STATA命令为：

tabi264428\284917,chi2lichi2

|col

row|123|Total

-----------+----------------------------------+----------

1|264428|98

|26.5344.9028.57|100.00

-----------+----------------------------------+----------

2|284917|94

|29.7952.1318.09|100.00

-----------+----------------------------------+----------

Total|549345|192

|28.1248.4423.44|100.00

Pearsonchi2(2)=2.9497Pr=0.229

likelihood-ratiochi2(2)=2.9760Pr=0.226

两种方法结果对应相同。另外数学上可以证明：当大样本是，卡方检验(Pearsonchi2test)与似然比检验(likelihoodratiochi2test)趋向一致。•Fisher精确检验

例：用新旧两种药治疗某种疾病(资料如下表所示)，试问：两药对该病的治愈率是否不同？(资料摘自医学卫生统计，金丕焕主编，p165)组别未治愈治愈合计

旧药4(2.2)2(3.8)6

新药0(1.8)5(3.2)5

合计4711

表中括号中的数为理论数(即：两种药的疗效无差异的无效假设Ho 成立时的期望频数)。

由于理论数均小于5，故必须用Fisher精确检验法。相应的STATA命令为：

tabi42\05,exactcolrow