流行病學(xué)資料的分析

STATA軟件應(yīng)用---流行病學(xué)資料分析非常簡(jiǎn)便、快速流行病學(xué)資料的分析第1頁(yè)流行病學(xué)資料分析病例對(duì)照研究資料

成組病例對(duì)照研究配對(duì)病例對(duì)照研究

隊(duì)列研究資料流行病學(xué)資料的分析第2頁(yè)例1(成組病例對(duì)照研究)某單位研究胸膜間皮瘤與接觸石棉關(guān)系，資料見下表。試對(duì)其進(jìn)行分析。

流行病學(xué)資料的分析第3頁(yè)成組病例對(duì)照研究分析

2檢驗(yàn)：接觸石棉與胸膜間皮瘤關(guān)系計(jì)算OR：反應(yīng)接觸石棉與胸膜間皮瘤關(guān)系大小計(jì)算OR95%置信區(qū)間：與

2檢驗(yàn)等價(jià)流行病學(xué)資料的分析第4頁(yè)成組病例對(duì)照分析cc病例變量

暴露變量[fw=頻數(shù)變量]cci#a#b#c#d流行病學(xué)資料的分析第5頁(yè)成組病例對(duì)照分析cccaseexposure[fw=f]

cci4036967流行病學(xué)資料的分析第6頁(yè)解釋：優(yōu)勢(shì)比OR=8.27,其95%CI（3.42，21.33）

2=28.94，P=0.0000

由此可見，接觸石棉者發(fā)生間皮瘤危險(xiǎn)性為未接觸者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流行病學(xué)資料的分析第7頁(yè)病例對(duì)照研究資料---聯(lián)絡(luò)強(qiáng)度指標(biāo)OR=ad/bc流行病學(xué)資料的分析第8頁(yè)例2（多暴露珠平成組病例對(duì)照研究資料）某醫(yī)院一項(xiàng)吸煙與膀胱癌關(guān)系病例對(duì)照研究得到以下結(jié)果，試分析吸煙量與膀胱癌是否存在劑量-反應(yīng)關(guān)系？流行病學(xué)資料的分析第9頁(yè)多暴露珠平成組病例對(duì)照研究目標(biāo)：多個(gè)暴露珠平劑量-反應(yīng)相關(guān)分析方法：趨勢(shì)

2檢驗(yàn)

H0：吸煙量與膀胱癌不存在劑量-反應(yīng)線性關(guān)系H1:吸煙量與膀胱癌存在劑量-反應(yīng)線性關(guān)系流行病學(xué)資料的分析第10頁(yè)多暴露珠平成組病例對(duì)照分析tabodds病例變量暴露變量[fw=頻數(shù)變量]，orbase(#)or：指定輸出各組or值base(#)：指定以哪組作為參考預(yù)計(jì)各組or值，默認(rèn)是第一組作為參考。流行病學(xué)資料的分析第11頁(yè)數(shù)據(jù)結(jié)構(gòu)流行病學(xué)資料的分析第12頁(yè)解釋：

OR同質(zhì)性檢驗(yàn)：

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

OR趨勢(shì)檢驗(yàn)：

2=15.14，P=0.0001，故可用認(rèn)為伴隨吸煙量增加患膀胱癌危險(xiǎn)性也增加。

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流行病學(xué)資料的分析第13頁(yè)例3(分層病例對(duì)照資料)在吸煙與肺癌發(fā)生關(guān)系研究中，年紀(jì)是一個(gè)混雜原因，試依據(jù)以下資料計(jì)算調(diào)整年紀(jì)后吸煙與肺癌發(fā)生關(guān)系優(yōu)勢(shì)比,并作假設(shè)檢驗(yàn)。

流行病學(xué)資料的分析第14頁(yè)分層病例對(duì)照研究目標(biāo)：按混雜原因分層，控制混雜原因影響。方法：Mantel-Haenszel分層分析

(1)公共優(yōu)勢(shì)比ORM-H預(yù)計(jì)算公式為：(2)對(duì)公共優(yōu)勢(shì)比進(jìn)行假設(shè)檢驗(yàn),

統(tǒng)計(jì)量為：流行病學(xué)資料的分析第15頁(yè)分層病例對(duì)照分析cc病例變量暴露變量[fw=頻數(shù)變量]，by(分層變量)mhodds病例變量暴露變量[fw=頻數(shù)變量]，by(分層變量)流行病學(xué)資料的分析第16頁(yè)數(shù)據(jù)結(jié)構(gòu)流行病學(xué)資料的分析第17頁(yè)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流行病學(xué)資料的分析第18頁(yè)結(jié)果解釋粗OR＝2.75；調(diào)整年紀(jì)后，ORMH＝2.91。

2MN=49.43，v=1，P=0.0000，差異有統(tǒng)計(jì)學(xué)意義，故可認(rèn)為調(diào)整年紀(jì)后，吸煙1包/天以上者患肺癌危險(xiǎn)度是1包/天以下者2.91倍。流行病學(xué)資料的分析第19頁(yè)例4(配對(duì)病例-對(duì)照研究)某單位為研究軟組織肉瘤與接觸苯氧乙酸或氯酚關(guān)系，作了一次病例-對(duì)照研究，結(jié)果見下表。

流行病學(xué)資料的分析第20頁(yè)配對(duì)四格表資料方法：McNemar檢驗(yàn)公式：b+c≤40b+c>40流行病學(xué)資料的分析第21頁(yè)配對(duì)病例對(duì)照分析mcc暴露-病例數(shù)暴露-對(duì)照數(shù)[fw=頻數(shù)變量]mcci#a#b#c#d流行病學(xué)資料的分析第22頁(yè)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流行病學(xué)資料的分析第23頁(yè)結(jié)果解釋OR=5.33，其95%CI為(1.53,28.56)。

本例樣本例數(shù)較少，應(yīng)計(jì)算確切概率。P=0.0044，差異有統(tǒng)計(jì)學(xué)意義，故可認(rèn)為接觸苯氧乙酸或氯酚者發(fā)生軟組織肉瘤危險(xiǎn)性為不接觸者5.33倍。

流行病學(xué)資料的分析第24頁(yè)流行病學(xué)資料分析病例對(duì)照研究資料

成組病例對(duì)照研究配對(duì)病例對(duì)照研究

隊(duì)列研究資料流行病學(xué)資料的分析第25頁(yè)例5（隊(duì)列研究資料）就下表資料分析母親乳汁中IgG抗體滴度高低與6個(gè)月以上嬰兒患呼吸道疾病關(guān)系。

流行病學(xué)資料的分析第26頁(yè)隊(duì)列研究分析

2檢驗(yàn)：分析母親乳汁中IgG抗體滴度與呼吸道疾病關(guān)系。計(jì)算RR：分析母親乳汁中IgG抗體滴度與呼吸道疾病關(guān)系大小。計(jì)算RR95%置信區(qū)間流行病學(xué)資料的分析第27頁(yè)隊(duì)列研究（Cohortstudies）分析cs病例變量

暴露變量[fw=頻數(shù)變量]csi#a#b#c#d流行病學(xué)資料的分析第28頁(yè)隊(duì)列研究分析cscaseexposure[fw=f]

csi165710流行病學(xué)資料的分析第29頁(yè)

|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)

相對(duì)危險(xiǎn)度RR=2.086，95%CI:(0.971,4.485)

歸因危險(xiǎn)度

AR=0.521，其95%CI（-0.030，0.777）

人群歸因危險(xiǎn)度PAR=0.397

2=4.82，P=0.0281，可認(rèn)為母親乳汁中IgG抗體滴度高低與嬰兒呼吸道疾患發(fā)生率相關(guān)，滴度較低者其嬰兒較易患呼吸道疾病。流行病學(xué)資料的分析第30頁(yè)隊(duì)列研究資料---聯(lián)絡(luò)強(qiáng)度指標(biāo)率差RD=P1-P0相對(duì)危險(xiǎn)度RR=P1/P0歸因危險(xiǎn)度

AR=(P1-P0)/P1人群歸因危險(xiǎn)度PAR=(P-P0)/P流行病學(xué)資料的分析第31頁(yè)例6（分層隊(duì)列研究）某人隨訪研究甲糖寧治療糖尿病人效果，其中年紀(jì)被認(rèn)為是可能混雜原因，需要對(duì)年紀(jì)進(jìn)行分層分析。流行病學(xué)資料的分析第32頁(yè)分層隊(duì)列研究分析cs病例變量暴露變量[fw=頻數(shù)變量]，by(分層變量)流行病學(xué)資料的分析第33頁(yè)數(shù)據(jù)結(jié)構(gòu)流行病學(xué)資料的分析第34頁(yè)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