Stata基本操作和数据分析入门第11章无序分类资料的统计分析Stata实现.doc

第十一章 无序分类资料的统计分析本章使用的STATA的命令：卡方检验tabulate varname1 varname2 weight, all chi2 exact expected卡方检验即时命令tabi #11 #12 . #21 #22 . .Yates校正卡方Exactcci #a #b #c #d配对四格表mcci #a #b #c #d , tb level(#) 例11-1 根据某地区的血型普查结果可知，该地区人群中血型为O的占30%，血型为A的占25%，血型为B的占35，血型为AB的占10%。研究者在邻近该地区的一个山区人群中进行一个血型的流行病调查，在该山区人群中随机抽样调查了200人，检测这些对象的血型，表11-1给出了血型检测的结果。问该山区人群与这个地区人群的血型分布是否一致？表11-1 山区人群血型抽样调查结果血型OABAB合计人数50705030200例11-2 某研究小组为研究慢支口服液II号对慢性支气管炎治疗效果，以口服消咳喘为对照进行了临床试验，试验组120人、对照组117人（两组受试者病程、病情等均衡），疗程2周，疗效见表11-3。问慢支口服液II号与消咳喘治疗慢性支气管炎的疗效是否相同？表11-3 试验组与对照组疗效组别有效无效合计有效率（%）试验组116412096.67对照组823511770.09合计1983923783.541建立检验假设，确定检验水准：，即两种药物治疗慢性支气管炎的疗效相同 ：，即两种药物治疗慢性支气管炎的疗效不同解：STATA数据如下：gxw1111621821242235STATA命令：tab g x weight=w,chi row结果：+-+| Key |-| frequency | row percentage |+-+ | x g | 1 2 | Total-+-+- 1 | 116 4 | 120 | 96.67 3.33 | 100.00 -+-+- 2 | 82 35 | 117 | 70.09 29.91 | 100.00 -+-+- Total | 198 39 | 237 | 83.54 16.46 | 100.00 Pearson chi2(1) = 30.4463 Pr = 0.000Pearson chi2(1) = 30.4463 Pr = 0.000，按=0.05水准拒绝，差别有统计学意义，可认为慢支口服液II号治疗慢性支气管炎有效率高于消咳喘。例11-3 为评价中西结合治疗抑郁发作的疗效。将187例患者随机分为2组，两组患者均选用阿咪替林西医综合治疗，中西医结合组在上述治疗的同时，再配合中医辨证治疗，根据中医辨证分型采用不同的方剂，治疗结果见表11-5，问两种治疗方案的疗效有无差别？ 表11-5 试验组与对照组疗效组别有效无效合计有效率（%）中西医结合组92（88.973）2（5.027）9497.87西医组85（88.027）8（4.973）9391.40注 ；括号内为理论频数 假设检验步骤：1建立检验假设，确定检验水准：，即两种治疗方案疗效相同 ：，即两种治疗方案疗效不同STATA命令：exactcci 92 2 85 8结果： Proportion | Exposed Unexposed | Total Exposed-+-+- Cases | 92 2 | 94 0.9787 Controls | 85 8 | 93 0.9140-+-+- Total | 177 10 | 187 0.9465 | | | Point estimate | 95% Conf. Interval |-+- | | Cornfields limits Odds ratio | 4.329412 | .8179455 . Adjusted | | 1.004989 . Unadjusted Attr. frac. ex. | .7690217 | -.2225754 . Adjusted | | .0049646 . Unadjusted Attr. frac. pop | .7526596 | +- chi2(1) = 3.87 Prchi2 = 0.0491 Yates adjusted chi2(1) = 2.70 Prchi2 = 0.1005本例需要用校正卡方，p=0.1005，两种治疗方案疗效相同。例11-4 为了解国产紫外线瞬间消毒器与进口高压蒸汽消毒机对牙科手机消毒灭菌的效果，将刚去腐揭卡过垢等待处理的牙科手机29个随机分为A、B两组，A组为紫外线消毒组，B组为高压蒸汽组。消毒前细菌培养均为阳性，消毒后细菌培养结果见下表。问两种消毒法消毒后细菌培养阳性率有无差别？表11-6 两种方法消毒后细菌培养结果组别阳性阴性合计A10515B11314合计111829由于总频数29小于40，对两组阳性率的比较宜采用Fisher精确概率检验，假设检验步骤如下：1建立检验假设，确定检验水准：，即两种方法消毒后细菌培养阳性率相同 ：，即两种方法消毒后细菌培养阳性率不同解：STATA数据：gxw11102111252213STATA命令：tab g x weight=w,chi row expect exact+-+| Key |-| frequency | expected frequency | row percentage |+-+ | x g | 1 2 | Total-+-+- 1 | 10 5 | 15 | 5.7 9.3 | 15.0 | 66.67 33.33 | 100.00 -+-+- 2 | 1 13 | 14 | 5.3 8.7 | 14.0 | 7.14 92.86 | 100.00 -+-+- Total | 11 18 | 29 | 11.0 18.0 | 29.0 | 37.93 62.07 | 100.00 Pearson chi2(1) = 10.8977 Pr = 0.001 Fishers exact = 0.002 1-sided Fishers exact = 0.001本例需要用Fishers确切概率法，p=0.002，两种方法消毒后细菌培养阳性率不同。例11-5 为探讨埃兹蛋白（Ezrin）在胃癌组织中的表达情况，采用免疫组化法检测50例胃癌组织、25例胃粘膜不典型增生和25例正常胃粘膜中Ezrin的表达，结果见表11-9。问不同胃组织Ezrin表达阳性率是否相同？表11-9 Ezrin在不同胃组织中的表达组别观测例数阳性例数阴性例数阳性率（%）正常胃粘膜257（15.250）18（9.750）28.0不典型增生2511（15.250）14（9.750）44.0胃癌组织5043（30.500）7（19.500）86.0合计100613961.0注 ：括号内为理论频数这是一个3个样本率的比较问题，假设检验步骤为：1建立检验假设，确定检验水准：，即3种不同胃组织Ezrin表达阳性率相等 ：、不全相等，即3种不同胃组织Ezrin表达阳性率不全相等解：STATA数据如下：ab x1172111314312182214327STATA命令为：tab a b weight=x,chi结果： | b a | 1 2 | Total-+-+- 1 | 7 18 | 25 2 | 11 14 | 25 3 | 43 7 | 50 -+-+- Total | 61 39 | 100 Pearson chi2(2) = 27.6166 Pr = 0.000，按=0.05水准拒绝，差别有统计学意义，可认为3种不同胃组织Ezrin表达阳性率不全相等。例11-6 为评价国产注射用头孢美唑钠（A）治疗中、重度呼吸系统细菌性感染性疾病的临床有效性及安全性，以先锋美他醇（B）为对照进行临床试验，入组受试者疾病类型构成情况见表11-10。问A、B两组受试者疾病类型总体构成有无差别？表11-10 两组受试者疾病类型组别急性扁桃体炎肺炎急支炎慢支炎急发支扩伴感染A5（6.042）21（19.636）21（21.650）20（20.643）5（4.028）B7（5.958）18（19.364）22（21.350）21（20.357）3（3.972）注 ：括号内为理论频数这是一个2组构成比比较的问题，其假设检验步骤为：1建立检验假设，确定检验水准：A、B受试者疾病类型总体构成相同：A、B受试者疾病类型总体构成不同解：STATA数据为：ab x115122113211420155217221823222421253STATA命令为：tab a b weight=x,chi结果为： | b a | 1 2 3 4 5 | Total-+-+- 1 | 5 21 21 20 5 | 72 2 | 7 18 22 21 3 | 71 -+-+- Total | 12 39 43 41 8 | 143 Pearson chi2(4) = 1.1048 Pr = 0.894P=0.894，A、B受试者疾病类型总体构成相同。例11-7 将100份样品一分为二，分别用含血培养基与无血培养基接种培养，观察弯曲菌检出情况，结果如表11-12所示。试问：两种培养基接种培养弯曲菌的阳性率是否相等？两种培养基培养结果间是否有关联性？将表11-12整理为表11-13形式表11-13 两种培养基弯曲菌检出结果无血培养基含血培养基合计+-[phone]合计6040100解：两种培养基接种培养弯曲菌的阳性率相同：两种培养基接种培养弯曲菌的阳性率不同利用STATA的即时命令mcci 52 17 8 23结果为： | Controls |Cases | Exposed Unexposed | Total-+-+- Exposed | 52 17 | 69 Unexposed | 8 23 | 31-+-+- Total | 60 40 | 100McNemars chi2(1) = 3.24 Prob chi2 = 0.0719Exact McNemar significance probability = 0.1078，按=0.05水准不拒绝，尚不能认为两种培养基接种培养弯曲菌的阳性率不相同。：两种培养基培养结果之间无关联性：两种培养基培养结果之间有关联性STATA命令为：tabi 52 178 23,chi结果为： | col row | 1 2 | Total-+-+- 1 | 52 17 | 69 2 | 8 23 | 31 -+-+- Total | 60 40 | 100 Pearson chi2(1) = 21.8872 Pr = 0.000，按=0.05水准拒绝，可认为两种培养基接种培养弯曲菌结果之间存在关联性。