第十一章无序分类资料的统计分析的Stata实现.docx
《第十一章无序分类资料的统计分析的Stata实现.docx》由会员分享，可在线阅读，更多相关《第十一章无序分类资料的统计分析的Stata实现.docx（15页珍藏版）》请在冰点文库上搜索。
第十一章无序分类资料的统计分析的Stata实现
第十一章无序分类资料的统计分析的Stata实现
本章使用的Stata的命令：
检验
tab变量1变量2,allchi2exactexpectedrow
检验即时命令
tabi#11#12[...]\#21#22[...][\...]
Yates’校正
Exactcci#a#b#c#d（见Stata7附加程序）
配对四格表
mcci#a#b#c#d[,tblevel（#）]
例11-1根据某地区的血型普查结果可知，该地区人群中血型为O的占30%，血型为A的占25%，血型为B的占35％，血型为AB的占10%。
研究者在邻近该地区的一个山区人群中进行一个血型的流行病调查，在该山区人群中随机抽样调查了200人，检测这些对象的血型，表11-1给出了血型检测的结果。
问该山区人群与这个地区人群的血型分布是否一致？
表11-1山区人群血型抽样调查结果
血型
O
A
B
AB
合计
人数
50
70
50
30
200
例11-2某研究小组为研究慢支口服液II号对慢性支气管炎治疗效果，以口服消咳喘为对照进行了临床试验，试验组120人、对照组117人（两组受试者病程、病情等均衡），疗程2周，疗效见表11-3。
问慢支口服液II号与消咳喘治疗慢性支气管炎的疗效是否相同？
表11-3试验组与对照组疗效
组别
有效
无效
合计
有效率（%）
试验组
116
4
120
96.67
对照组
82
35
117
70.09
合计
198
39
237
83.54
1．建立检验假设，确定检验水准
：
，即两种药物治疗慢性支气管炎的疗效相同
：
，即两种药物治疗慢性支气管炎的疗效不同
Stata数据如下：
g
x
w
1
1
116
2
1
82
1
2
4
2
2
35
Stata命令：
tabgx[weight=w],chirow
结果：
+----------------+
|Key|
|----------------|
|frequency|
|rowpercentage|
+----------------+
|x
g|12|Total
-----------+----------------------+----------
1|1164|120
|96.673.33|100.00
-----------+----------------------+----------
2|8235|117
|70.0929.91|100.00
-----------+----------------------+----------
Total|19839|237
|83.5416.46|100.00
Pearsonchi2
（1）=30.4463Pr=0.000
Pearsonchi2
（1）=30.4463Pr=0.000，
，按
=0.05水准拒绝
，差别有统计学意义，可认为慢支口服液II号治疗慢性支气管炎有效率高于消咳喘。
例11-3为评价中西结合治疗抑郁发作的疗效。
将187例患者随机分为2组，两组患者均选用阿咪替林西医综合治疗，中西医结合组在上述治疗的同时，再配合中医辨证治疗，根据中医辨证分型采用不同的方剂，治疗结果见表11-5，问两种治疗方案的疗效有无差别？
表11-5试验组与对照组疗效
组别
有效
无效
合计
有效率（%）
中西医结合组
92（88.973）
2（5.027）
94
97.87
西医组
85（88.027）
8（4.973）
93
91.40
注；括号内为理论频数
例11-3假设检验步骤
1．建立检验假设，确定检验水准
：
，即两种治疗方案疗效相同
：
，即两种治疗方案疗效不同
Stata命令：
exactcci922858
结果：
Proportion
|ExposedUnexposed|TotalExposed
-----------------+------------------------+------------------------
Cases|922|940.9787
Controls|858|930.9140
-----------------+------------------------+------------------------
Total|17710|1870.9465
||
|Pointestimate|[95%Conf.Interval]
|------------------------+----------------------
||Cornfield'slimits
Oddsratio|4.329412|.8179455.Adjusted
||1.004989.Unadjusted
Attr.frac.ex.|.7690217|-.2225754.Adjusted
||.0049646.Unadjusted
Attr.frac.pop|.7526596|
+-----------------------------------------------
chi2
（1）=3.87Pr>chi2=0.0491
Yates'adjustedchi2
（1）=2.70Pr>chi2=0.1005
本例需要用校正卡方，p=0.1005，两种治疗方案疗效的差异无统计学意义。
例11-4为了解国产紫外线瞬间消毒器与进口高压蒸汽消毒机对牙科手机消毒灭菌的效果，将刚去腐揭卡过垢等待处理的牙科手机29个随机分为A、B两组，A组为紫外线消毒组，B组为高压蒸汽组。
消毒前细菌培养均为阳性，消毒后细菌培养结果见下表。
问两种消毒法消毒后细菌培养阳性率有无差别？
表11-6两种方法消毒后细菌培养结果
组别
阳性
阴性
合计
A
10
5
15
B
1
13
14
合计
11
18
29
由于总频数29小于40，对两组阳性率的比较宜采用Fisher精确概率检验，假设检验步骤如下：
1．建立检验假设，确定检验水准
：
，即两种方法消毒后细菌培养阳性率相同
：
，即两种方法消毒后细菌培养阳性率不同
Stata数据：
g
x
w
1
1
10
2
1
1
1
2
5
2
2
13
Stata命令：
tabgx[weight=w],chirowexpectexact
+--------------------+
|Key|
|--------------------|
|frequency|
|expectedfrequency|
|rowpercentage|
+--------------------+
|x
g|12|Total
-----------+----------------------+----------
1|105|15
|5.79.3|15.0
|66.6733.33|100.00
-----------+----------------------+----------
2|113|14
|5.38.7|14.0
|7.1492.86|100.00
-----------+----------------------+----------
Total|1118|29
|11.018.0|29.0
|37.9362.07|100.00
Pearsonchi2
（1）=10.8977Pr=0.001
Fisher'sexact=0.002
1-sidedFisher'sexact=0.001
本例需要用Fisher's确切概率法，p=0.002，两种方法消毒后细菌培养阳性率不同。
例11-5为探讨埃兹蛋白（Ezrin）在胃癌组织中的表达情况，采用免疫组化法检测50例胃癌组织、25例胃粘膜不典型增生和25例正常胃粘膜中Ezrin的表达，结果见表11-9。
问不同胃组织Ezrin表达阳性率是否相同？
表11-9Ezrin在不同胃组织中的表达
组别
观测例数
阳性例数
阴性例数
阳性率（%）
正常胃粘膜
25
7（15.250）
18（9.750）
28.0
不典型增生
25
11（15.250）
14（9.750）
44.0
胃癌组织
50
43（30.500）
7（19.500）
86.0
合计
100
61
39
61.0
注：
括号内为理论频数
这是一个3个样本率的比较问题，假设检验步骤为：
1．建立检验假设，确定检验水准
：
，即3种不同胃组织Ezrin表达阳性率相等
：
、
、
不全相等，即3种不同胃组织Ezrin表达阳性率不全相等
Stata数据如下：
a
b
x
1
1
7
2
1
11
3
1
43
1
2
18
2
2
14
3
2
7
Stata命令为：
tabab[weight=x],chi
结果：
|b
a|12|Total
-----------+----------------------+----------
1|718|25
2|1114|25
3|437|50
-----------+----------------------+----------
Total|6139|100
Pearsonchi2
（2）=27.6166Pr=0.000
，按
=0.05水准拒绝
，差别有统计学意义，可认为3种不同胃组织Ezrin表达阳性率不全相等。
例11-6为评价国产注射用头孢美唑钠（A）治疗中、重度呼吸系统细菌性感染性疾病的临床有效性及安全性，以先锋美他醇（B）为对照进行临床试验，入组受试者疾病类型构成情况见表11-10。
问A、B两组受试者疾病类型总体构成有无差别？
表11-10两组受试者疾病类型
组别
急性扁桃体炎
肺炎
急支炎
慢支炎急发
支扩伴感染
A
5（6.042）
21（19.636）
21（21.650）
20（20.643）
5（4.028）
B
7（5.958）
18（19.364）
22（21.350）
21（20.357）
3（3.972）
注：
括号内为理论频数
这是一个2组构成比比较的问题，其假设检验步骤为：
1．建立检验假设，确定检验水准
：
A、B受试者疾病类型总体构成相同
：
A、B受试者疾病类型总体构成不同
Stata数据为：
a
b
x
1
1
5
1
2
21
1
3
21
1
4
20
1
5
5
2
1
7
2
2
18
2
3
22
2
4
21
2
5
3
Stata命令为：
tabab[weight=x],chi
结果为：
|b
a|12345|Total
-----------+-------------------------------------------------------+----------
1|52121205|72
2|71822213|71
-----------+-------------------------------------------------------+----------
Total|123943418|143
Pearsonchi2（4）=1.1048Pr=0.894
P=0.894，A、B受试者疾病类型总体构成相同。
例11-7将100份样品一分为二，分别用含血培养基与无血培养基接种培养，观察弯曲菌检出情况，结果如表11-12所示。
试问：
两种培养基接种培养弯曲菌的阳性率是否相等？
两种培养基培养结果间是否有关联性？
将表11-12整理为表11-13形式
表11-13两种培养基弯曲菌检出结果
无血培养基
含血培养基
合计
+
-
+
52
17
69
-
8
23
31
合计
60
40
100
：
两种培养基接种培养弯曲菌的阳性率相同
：
两种培养基接种培养弯曲菌的阳性率不同
利用Stata的即时命令
mcci5217823
结果为：
|Controls|
Cases|ExposedUnexposed|Total
-----------------+------------------------+----------
Exposed|5217|69
Unexposed|823|31
-----------------+------------------------+----------
Total|6040|100
McNemar'schi2
（1）=3.24Prob>chi2=0.0719
ExactMcNemarsignificanceprobability=0.1078
，按
=0.05水准不拒绝
，尚不能认为两种培养基接种培养弯曲菌的阳性率不相同。
：
两种培养基培养结果之间无关联性
：
两种培养基培养结果之间有关联性
Stata命令为：
tabi5217\823,chi
结果为：
|col
row|12|Total
-----------+----------------------+----------
1|5217|69
2|823|31
-----------+----------------------+----------
Total|6040|100
Pearsonchi2
（1）=21.8872Pr=0.000
，按
=0.05水准拒绝
，可认为两种培养基接种培养弯曲菌结果之间存在关联性。
升级会员 