第十一章无序分类资料的统计分析的Stata实现

第十一章无序分类资料的统计分析的Stata实现本章使用的Stata的命令：2检验tab变量1变量2,allchi2exactexpectedrow2检验即时命令tabi#11#12[...]\#21#22[...][\...]2Exactcci#a#b#c#d（见Stata7附加程序）Yates’校正配对四格表mcci#a#b#c#d[,tblevel(#)]例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．建立检验假设，确定检验水准H0：12，即两种药物治疗慢性支气管炎的疗效相同H1：12，即两种药物治疗慢性支气管炎的疗效不同0.05Stata数据如下：gx11211222Stata命令：w11682435tabgx[weight=w],chirow结果：+----------------+|Key||----------------||frequency||rowpercentage|+----------------+|xg|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.00Pearsonchi2(1)=30.4463Pr=0.000Pearsonchi2(1)=30.4463Pr=0.000，P0.05，按=0.05水准拒绝H0，差别有统计学意义，可认为慢支口服液II号治疗慢性支气管炎有效率高于消咳喘。例11-3为评价中西结合治疗抑郁发作的疗效。将187例患者随机分为2组，两组患者均选用阿咪替林西医综合治疗，中西医结合组在上述治疗的同时，再配合中医辨证治疗，根据中医辨证分型采用不同的方剂，治疗结果见表11-5，问两种治疗 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 的疗效有无差别？表11-5试验组与对照组疗效组别有效无效合计有效率（%）中西医结合组92（88.973）2（5.027）9497.87西医组85（88.027）8（4.973）9391.40注；括号内为理论频数例11-3假设检验步骤1．建立检验假设，确定检验水准H0：12，即两种治疗方案疗效相同H1：12，即两种治疗方案疗效不同0.05Stata命令：exactcci922858结果：Proportion|ExposedUnexposed|TotalExposed-----------------+------------------------+------------------------Cases|922|940.9787Controls|858|930.9140-----------------+------------------------+------------------------Total|17710|1870.9465|||Pointestimate|[95%Conf.Interval]|------------------------+----------------------||Cornfield'slimitsOddsratio|4.329412|.8179455.Adjusted||1.004989.UnadjustedAttr.frac.ex.||.7690217|-.2225754|.0049646..AdjustedUnadjustedAttr.frac.pop|.7526596|+-----------------------------------------------chi2(1)=3.87Pr>chi2=0.0491Yates'adjustedchi2(1)=2.70Pr>chi2=0.1005本例需要用校正卡方，p=0.1005，两种治疗方案疗效的差异无统计学意义。例11-4为了解国产紫外线瞬间消毒器与进口高压蒸汽消毒机对牙科手机消毒灭菌的效果，将刚去腐揭卡过垢等待处理的牙科手机29个随机分为A、B两组，A组为紫外线消毒组，B组为高压蒸汽组。消毒前细菌培养均为阳性，消毒后细菌培养结果见下表。问两种消毒法消毒后细菌培养阳性率有无差别？表11-6两种 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 消毒后细菌培养结果组别阳性阴性合计A10515B11314合计111829由于总频数29小于40，对两组阳性率的比较宜采用Fisher精确概率检验，假设检验步骤如下：1．建立检验假设，确定检验水准H0：12，即两种方法消毒后细菌培养阳性率相同H1：12，即两种方法消毒后细菌培养阳性率不同0.05Stata数据：gx11211222w101513Stata命令：tabgx[weight=w],chirowexpectexact+--------------------+|Key||--------------------||frequency||expectedfrequency||rowpercentage|+--------------------+|xg|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.00Pearsonchi2(1)=10.8977Pr=0.001Fisher'sexact=0.0021-sidedFisher'sexact=0.001本例需要用Fisher's确切概率法，p=0.002，两种方法消毒后细菌培养阳性率不同。例11-5为探讨埃兹蛋白（Ezrin）在胃癌组织中的表达情况，采用免疫组化法检测50例胃癌组织、25例胃粘膜不典型增生和25例正常胃粘膜中Ezrin的表达，结果见表11-9。问不同胃组织Ezrin表达阳性率是否相同？表11-9Ezrin在不同胃组织中的表达组别观测例数阳性例数阴性例数阳性率（%）正常胃粘膜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．建立检验假设，确定检验水准H0：123，即3种不同胃组织Ezrin表达阳性率相等H1：1、2、3不全相等，即3种不同胃组织Ezrin表达阳性率不全相等0.05Stata数据如下：abx1172111314312182214327Stata命令为：tabab[weight=x],chi结果：|ba|12|Total-----------+----------------------+----------1|718|252|1114|253|437|50-----------+----------------------+----------Total|6139|100Pearsonchi2(2)=27.6166Pr=0.000P0.05，按=0.05水准拒绝H0，差别有统计学意义，可认为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．建立检验假设，确定检验水准H0：A、B受试者疾病类型总体构成相同H1：A、B受试者疾病类型总体构成不同0.05Stata数据为：abx115122113211420155217221823222421253Stata命令为：tabab[weight=x],chi结果为：|ba|12345|Total-----------+-------------------------------------------------------+----------1|52121205|722|71822213|71-----------+-------------------------------------------------------+----------Total|123943418|143Pearsonchi2(4)=1.1048Pr=0.894P=0.894，A、B受试者疾病类型总体构成相同。例11-7将100份样品一分为二，分别用含血培养基与无血培养基接种培养，观察弯曲菌检出情况，结果如表11-12所示。试问：两种培养基接种培养弯曲菌的阳性率是否相等？两种培养基培养结果间是否有关联性？将表11-12整理为表11-13形式表11-13两种培养基弯曲菌检出结果含血培养基无血培养基合计+-+521769-82331合计6040100H0：两种培养基接种培养弯曲菌的阳性率相同H1：两种培养基接种培养弯曲菌的阳性率不同0.05利用Stata的即时命令mcci5217823结果为：|Controls|Cases|ExposedUnexposed|Total-----------------+------------------------+----------Exposed|5217|69Unexposed|823|31-----------------+------------------------+----------Total|6040|100McNemar'schi2(1)=3.24Prob>chi2=0.0719ExactMcNemarsignificanceprobability=0.10780.05，按=0.05水准不拒绝H0，尚不能认为两种培养基接种培养弯曲菌的阳性率不相同。H0：两种培养基培养结果之间无关联性H1：两种培养基培养结果之间有关联性0.05Stata命令为：tabi5217\823,chi结果为：|colrow|12|Total-----------+----------------------+----------1|5217|692|823|31-----------+----------------------+----------Total|6040|100Pearsonchi2(1)=21.8872Pr=0.0000.05，按=0.05水准拒绝H0，可认为两种培养基接种培养弯曲菌结果之间存在关联性。