流行病學資料的分析

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

成組病例對照研究配對病例對照研究

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

流行病學資料的分析第3頁成組病例對照研究分析

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

2檢驗等價流行病學資料的分析第4頁成組病例對照分析cc病例變量

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

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

2=28.94，P=0.0000

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

2檢驗

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

OR同質(zhì)性檢驗：

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

流行病學資料的分析第14頁分層病例對照研究目標：按混雜原因分層，控制混雜原因影響。方法：Mantel-Haenszel分層分析

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

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

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

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

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

流行病學資料的分析第24頁流行病學資料分析病例對照研究資料

成組病例對照研究配對病例對照研究

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

流行病學資料的分析第26頁隊列研究分析

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

暴露變量[fw=頻數(shù)變量]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抗體滴度高低與嬰兒呼吸道疾患發(fā)生率相關(guān)，滴度較低者其嬰兒較易患呼吸道疾病。流行病學資料的分析第30頁隊列研究資料---聯(lián)絡(luò)強度指標率差RD=P1-P0相對危險度RR=P1/P0歸因危險度

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