流行病学的分析.ppt

STATA软件的应用,-流行病学资料的分析,非常简便、快速,流行病学资料的分析,病例对照研究资料成组病例对照研究配对病例对照研究队列研究资料,例1(成组病例对照研究)某单位研究胸膜间皮瘤与接触石棉的关系，资料见下表。试对其进行分析。,成组病例对照研究分析,2检验：接触石棉与胸膜间皮瘤的关系计算OR：反映接触石棉与胸膜间皮瘤的关系大小计算OR95%的置信区间：与2检验等价,成组病例对照分析,cc病例变量暴露变量fw=频数变量cci#a#b#c#d,成组病例对照分析,cccaseexposurefw=fcci4036967,解释：优势比OR=8.27,其95%CI（3.42，21.33）2=28.94，P=0.0000由此可见，接触石棉者发生间皮瘤的危险性为未接触者的8.27倍。,Proportion|ExposedUnexposed|TotalExposed-+-+-Cases|4036|760.5263Controls|967|760.1184-+-+-Total|49103|1520.3224|Pointestimate|95%Conf.Interval|-+-Oddsratio|8.271605|3.419321.33091(exact)Attr.frac.ex.|.8791045|.7075425.9531197(exact)Attr.frac.pop|.4626866|+-chi2(1)=28.94Prchi2=0.0000,病例对照研究资料-联系强度指标,OR=ad/bc,例2（多暴露水平成组病例对照研究资料）某医院一项吸烟与膀胱癌关系的病例对照研究得到如下结果，试分析吸烟量与膀胱癌是否存在剂量-反应关系？,多暴露水平成组病例对照研究,目的：多个暴露水平的剂量-反应相关分析方法：趋势2检验H0：吸烟量与膀胱癌不存在剂量-反应线性关系H1:吸烟量与膀胱癌存在剂量-反应线性关系,多暴露水平成组病例对照分析,tabodds病例变量暴露变量fw=频数变量，orbase(#)or：指定输出各组的or值base(#)：指定以哪组作为参照估计各组的or值，默认是第一组作为参照。,数据结构,解释：OR的同质性检验：2=17.52，P=0.0006，即各组OR不全同。OR的趋势检验：2=15.14，P=0.0001，故可用认为随着吸烟量的增加患膀胱癌的危险性也增加。,taboddscasesmokefw=f,or-smoke|OddsRatiochi2Pchi295%Conf.Interval-+-0|1.000000.1|0.8541670.110.73580.3416232.1356872|1.6340581.550.21300.7485253.5672103|3.54807713.400.00031.7194707.321353-Testofhomogeneity(equalodds):chi2(3)=17.52Prchi2=0.0006Scoretestfortrendofodds:chi2(1)=15.14Prchi2=0.0001,例3(分层病例对照资料)在吸烟与肺癌发生关系的研究中，年龄是一个混杂因素，试根据下列资料计算调整年龄后吸烟与肺癌发生关系的优势比,并作假设检验。,分层病例对照研究,目的：按混杂因素分层，控制混杂因素的影响。方法：Mantel-Haenszel分层分析(1)公共优势比ORM-H的估计算公式为：(2)对公共优势比进行假设检验,统计量为：,分层病例对照分析,cc病例变量暴露变量fw=频数变量，by(分层变量)mhodds病例变量暴露变量fw=频数变量，by(分层变量),数据结构,cccasesmokefw=f,by(age),age|OR95%Conf.IntervalM-HWeight-+-0|2.2246581.4373823.43067214.28571(exact)1|4.1749072.5546596.8213947.819588(exact)-+-Crude|2.7474562.0116883.743611(exact)M-Hcombined|2.9145442.1433083.963297-Testofhomogeneity(M-H)chi2(1)=3.94Prchi2=0.0470TestthatcombinedOR=1:Mantel-Haenszelchi2(1)=49.43Prchi2=0.0000,结果解释,粗的OR2.75；调整年龄后，ORMH2.91。2MN=49.43，v=1，P=0.0000，差别有统计学意义，故可认为调整年龄后，吸烟1包/天以上者患肺癌的危险度是1包/天以下者的2.91倍。,例4(配对病例-对照研究)某单位为研究软组织肉瘤与接触苯氧乙酸或氯酚的关系，作了一次病例-对照研究，结果见下表。,配对四格表资料,方法：McNemar检验公式：,b+c40,b+c40,配对病例对照分析,mcc暴露-病例数暴露-对照数fw=频数变量mcci#a#b#c#d,mcci316330|Controls|Cases|ExposedUnexposed|Total-+-+-Exposed|316|19Unexposed|330|33-+-+-Total|646|52McNemarschi2(1)=8.89Probchi2=0.0029ExactMcNemarsignificanceprobability=0.0044ProportionwithfactorCases.3653846Controls.115384695%Conf.Interval-difference.25.0811853.4188147ratio3.1666671.4226597.048617rel.diff.2826087.1253028.4399146oddsratio5.3333331.52662728.56284(exact),b+cchi2=0.0281,解释：率差RD=0.3623，95%CI:0.059,0.666)相对危险度RR=2.086，95%CI:(0.971,4.485)归因危险度AR=0.521，其95%CI（-0.030，0.777）人群归因危险度PAR=0.3972=4.82，P=0.0281，可认为母亲乳汁中IgG抗体滴度高低与婴儿呼吸道疾患发生率有关，滴度较低者其婴儿较易患呼吸道疾病。,队列研究资料-联系强度指标,率差RD=P1-P0相对危险度RR=P1/P0归因危险度AR=(P1-P0)/P1人群归因危险度PAR=(P-P0)/P,例6（分层队列研究）某人随访研究甲糖宁治疗糖尿病人的效果，其中年龄被认为是可能的混杂因素，需要对年龄进行分层分析。,分层队列研究分析,cs病例变量暴露变量fw=频数变量，by(分层变量),数据结构,csdexposefw=f,by(age),age|RR95%Conf.IntervalM-HWeight-+-0|1.811321.61120445.3678982.3451331|1.332908.74715142.377898.124352-+-Crude|1.505602.89179242.54189M-Hcombined|1.440071.86379142.400817-Testofhomogeneity(M-H)chi2(1)=0.240Prchi2=0.6244,粗的RR1.51；调整年龄后，RRMH1.44，其95%CI：0.86-2.40，差别无统计学意义，故调整年龄后，还不能认为甲糖宁治疗糖尿病人有效。,小结,病例对照研究成组病例对照：cccaseexposefw=f/cci#a#b#c#dcccaseexposefw=f，by（混杂变量）分层分析taboddscaseexposefw=f，orbase(#)多暴露水平配对病例对照：mccrcfw=f/mcci#a#b#c#d队列研究cscaseexposefw=f/csi#a#b#c#dcscaseexposefw=f，by（混杂变量）分层分析,课堂练习1：,课堂练习2：例18-5在食管癌与饮酒关系的分析中，年龄可能是混杂因素。将每日饮酒量的资料按年龄组分解成4个四格表(表18-3)，估计出各层的OR值，可见不同年龄组的OR有所不同。试进行Mantel-Haenszel分层分析。,课堂练