流行病学资料的分析

STATA软件应用---流行病学资料分析非常简便、快速流行病学资料的分析第1页流行病学资料分析病例对照研究资料

成组病例对照研究配对病例对照研究

队列研究资料流行病学资料的分析第2页例1(成组病例对照研究)某单位研究胸膜间皮瘤与接触石棉关系，资料见下表。试对其进行分析。

流行病学资料的分析第3页成组病例对照研究分析

2检验：接触石棉与胸膜间皮瘤关系计算OR：反应接触石棉与胸膜间皮瘤关系大小计算OR95%置信区间：与

2检验等价流行病学资料的分析第4页成组病例对照分析cc病例变量

暴露变量[fw=频数变量]cci#a#b#c#d流行病学资料的分析第5页成组病例对照分析cccaseexposure[fw=f]

cci4036967流行病学资料的分析第6页解释：优势比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.94Pr>chi2=0.0000流行病学资料的分析第7页病例对照研究资料---联络强度指标OR=ad/bc流行病学资料的分析第8页例2（多暴露珠平成组病例对照研究资料）某医院一项吸烟与膀胱癌关系病例对照研究得到以下结果，试分析吸烟量与膀胱癌是否存在剂量-反应关系？流行病学资料的分析第9页多暴露珠平成组病例对照研究目标：多个暴露珠平剂量-反应相关分析方法：趋势

2检验

H0：吸烟量与膀胱癌不存在剂量-反应线性关系H1:吸烟量与膀胱癌存在剂量-反应线性关系流行病学资料的分析第10页多暴露珠平成组病例对照分析tabodds病例变量暴露变量[fw=频数变量]，orbase(#)or：指定输出各组or值base(#)：指定以哪组作为参考预计各组or值，默认是第一组作为参考。流行病学资料的分析第11页数据结构流行病学资料的分析第12页解释：

OR同质性检验：

2=17.52，P=0.0006，即各组OR不全同。

OR趋势检验：

2=15.14，P=0.0001，故可用认为伴随吸烟量增加患膀胱癌危险性也增加。

taboddscasesmoke[fw=f],or

---------------------------------------------------------------------------smoke|OddsRatiochi2P>chi2[95%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.52Pr>chi2=0.0006Scoretestfortrendofodds:chi2(1)=15.14Pr>chi2=0.0001流行病学资料的分析第13页例3(分层病例对照资料)在吸烟与肺癌发生关系研究中，年纪是一个混杂原因，试依据以下资料计算调整年纪后吸烟与肺癌发生关系优势比,并作假设检验。

流行病学资料的分析第14页分层病例对照研究目标：按混杂原因分层，控制混杂原因影响。方法：Mantel-Haenszel分层分析

(1)公共优势比ORM-H预计算公式为：(2)对公共优势比进行假设检验,

统计量为：流行病学资料的分析第15页分层病例对照分析cc病例变量暴露变量[fw=频数变量]，by(分层变量)mhodds病例变量暴露变量[fw=频数变量]，by(分层变量)流行病学资料的分析第16页数据结构流行病学资料的分析第17页cccasesmoke[fw=f],by(age)

age|OR[95%Conf.Interval]M-HWeight-----------------+-------------------------------------------------0|2.2246581.4373823.43067214.28571(exact)1|4.1749072.5546596.8213947.819588(exact)-----------------+-------------------------------------------------Crude|2.7474562.0116883.743611(exact)M-Hcombined|2.914544

2.1433083.963297-------------------------------------------------------------------Testofhomogeneity(M-H)chi2(1)=3.94Pr>chi2=0.0470TestthatcombinedOR=1:

Mantel-Haenszelchi2(1)=49.43

Pr>chi2=0.0000流行病学资料的分析第18页结果解释粗OR＝2.75；调整年纪后，ORMH＝2.91。

2MN=49.43，v=1，P=0.0000，差异有统计学意义，故可认为调整年纪后，吸烟1包/天以上者患肺癌危险度是1包/天以下者2.91倍。流行病学资料的分析第19页例4(配对病例-对照研究)某单位为研究软组织肉瘤与接触苯氧乙酸或氯酚关系，作了一次病例-对照研究，结果见下表。

流行病学资料的分析第20页配对四格表资料方法：McNemar检验公式：b+c≤40b+c>40流行病学资料的分析第21页配对病例对照分析mcc暴露-病例数暴露-对照数[fw=频数变量]mcci#a#b#c#d流行病学资料的分析第22页mcci316330

|Controls|Cases|ExposedUnexposed|Total-----------------+------------------------+----------Exposed|316|19Unexposed|330|33-----------------+------------------------+----------Total|646|52McNemar'schi2(1)=8.89Prob>chi2=0.0029ExactMcNemarsignificanceprobability=0.0044ProportionwithfactorCases.3653846Controls.1153846[95%Conf.Interval]-----------------------------difference.25.0811853.4188147ratio3.1666671.4226597.048617rel.diff..2826087.1253028.4399146oddsratio5.3333331.52662728.56284(exact)b+c<40b+c≥40流行病学资料的分析第23页结果解释OR=5.33，其95%CI为(1.53,28.56)。

本例样本例数较少，应计算确切概率。P=0.0044，差异有统计学意义，故可认为接触苯氧乙酸或氯酚者发生软组织肉瘤危险性为不接触者5.33倍。

流行病学资料的分析第24页流行病学资料分析病例对照研究资料

成组病例对照研究配对病例对照研究

队列研究资料流行病学资料的分析第25页例5（队列研究资料）就下表资料分析母亲乳汁中IgG抗体滴度高低与6个月以上婴儿患呼吸道疾病关系。

流行病学资料的分析第26页队列研究分析

2检验：分析母亲乳汁中IgG抗体滴度与呼吸道疾病关系。计算RR：分析母亲乳汁中IgG抗体滴度与呼吸道疾病关系大小。计算RR95%置信区间流行病学资料的分析第27页队列研究（Cohortstudies）分析cs病例变量

暴露变量[fw=频数变量]csi#a#b#c#d流行病学资料的分析第28页队列研究分析cscaseexposure[fw=f]

csi165710流行病学资料的分析第29页

|ExposedUnexposed|Total

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

Cases|165|21

Noncases|710|17

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

Total|2315|38

||

Risk|.6956522.3333333|.5526316

|Pointestimate|[95%Conf.Interval]

|------------------------+----------------------

Riskdifference|.3623188|.0585558.6660819

Riskratio|2.086957|.97110584.484977

Attr.frac.ex.|.5208333|-.0297539.7770334

Attr.frac.pop|.3968254|

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

chi2(1)=4.82Pr>chi2=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.397

2=4.82，P=0.0281，可认为母亲乳汁中IgG抗体滴度高低与婴儿呼吸道疾患发生率相关，滴度较低者其婴儿较易患呼吸道疾病。流行病学资料的分析第30页队列研究资料---联络强度指标率差RD=P1-P0相对危险度RR=P1/P0归因危险度

AR=(P1-P0)/P1人群归因危险度PAR=(P-P0)/P流行病学资料的分析第31页例6（分层队列研究）某人随访研究甲糖宁治疗糖尿病人效果，其中年纪被认为是可能混杂原因，需要对年纪进行分层分析。流行病学资料的分析第32页分层队列研究分析cs病例变量暴露变量[fw=频数变量]，by(分层变量)流行病学资料的分析第33页数据结构流行病学资料的分析第34页csdexpose[fw=f],by(age)age|RR[95%Conf.Interval]M-HWeight-----------------+-------------------------------------------------0|1.811321.61120445.3678982.3451331|1.332908.74715142.377898.124352-----------------+-------------------------------------------------

Crude|1.505602.89179242.54189

M-Hcombined|1.440071.86379142.400817-------------------------------------------------------------------Testofhomoge