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

1、STATA軟件的運(yùn)用- 流行病學(xué)資料的分析非常簡便、非常簡便、快速快速流行病學(xué)資料的分析流行病學(xué)資料的分析o病例對照研討資料o 成組病例對照研討o 配對病例對照研討 o隊(duì)列研討資料 例例1(成組病例對照研討成組病例對照研討) 某單位研討胸膜間皮瘤與接觸某單位研討胸膜間皮瘤與接觸石棉的關(guān)系，資料見下表。試對其進(jìn)展分析。石棉的關(guān)系，資料見下表。試對其進(jìn)展分析。 成組病例對照研討分析成組病例對照分析occ 病例變量 暴露變量 fw=頻數(shù)變量occi #a #b #c #d成組病例對照分析occ case exposure fw=foocci 40 36 9 67解釋：優(yōu)勢比OR = 8.27, 其9

2、5% CI3.42，21.33 2 = 28.94 ， P = 0.0000 由此可見，接觸石棉者發(fā)生間皮瘤的危險(xiǎn)性為未接觸者的8.27倍。 Proportion Proportion | Exposed Unexposed | Total Exposed | Exposed Unexposed | Total Exposed-+-+-+-+- Cases | 40 36 | 76 0.5263 Cases | 40 36 | 76 0.5263 Controls | 9 67 | 76 0.1184 Controls | 9 67 | 76 0.1184-+-+-+-+- Total | 4

3、9 103 | 152 0.3224 Total | 49 103 | 152 0.3224 | | | | | Point estimate | 95% Conf. Interval | Point estimate | 95% Conf. Interval |-+- |-+- Odds ratio | 8.271605 | 3.4193 21.33091 (exact) Odds ratio | 8.271605 | 3.4193 21.33091 (exact) Attr. frac. ex. | .8791045 | .7075425 .9531197 (exact) Attr. fr

4、ac. ex. | .8791045 | .7075425 .9531197 (exact) Attr. frac. pop | .4626866 | Attr. frac. pop | .4626866 | +- +- chi2(1) = 28.94 Prchi2 = 0.0000 chi2(1) = 28.94 Prchi2 = 0.0000病例對照研討資料-聯(lián)絡(luò)強(qiáng)度目的oOR= ad/bc例2多暴露程度成組病例對照研討資料某醫(yī)院一項(xiàng)吸煙與膀胱癌關(guān)系的病例對照研討得到如下結(jié)果，試分析吸煙量與膀胱癌能否存在劑量-反響關(guān)系？多暴露程度成組病例對照研討o目的：多個(gè)暴露程度的劑量目的：多個(gè)暴露程度

5、的劑量-反響相關(guān)分析反響相關(guān)分析 o方法：趨勢方法：趨勢2檢驗(yàn)檢驗(yàn)o H0：吸煙量與膀胱癌不存在劑量：吸煙量與膀胱癌不存在劑量-反響反響線性關(guān)系線性關(guān)系 o H1: 吸煙量與膀胱癌存在劑量吸煙量與膀胱癌存在劑量-反響線性反響線性關(guān)系關(guān)系1)()()(0202013022vxmxmnnnnnTaxkikiiiiikiiii多暴露程度成組病例對照分析otabodds 病例變量 暴露變量 fw=頻數(shù)變量，or base(#) o or：指定輸出各組的or值o base (#)：指定以哪組作為參照估計(jì)各組的or值，默許是第一組作為參照。數(shù)據(jù)構(gòu)造解釋：OR的同質(zhì)性檢驗(yàn)：2 =17.52，P=0.0006

6、，即各組OR不全同。 OR的趨勢檢驗(yàn)：2 =15.14，P=0.0001，故可用以為隨著吸煙量的添加患膀胱癌的危險(xiǎn)性也添加。 tabodds case smoke fw=f, or tabodds case smoke fw=f, or - smoke | Odds Ratio chi2 smoke | Odds Ratio chi2 Pchi2 95% Conf. IntervalPchi2 95% Conf. Interval-+-+- 0 | 0 | 1.000001.000000 . . . 0 . . . . . 1 | 0.854167 0.11 1 | 0.854167 0.11

7、 0.7358 0.341623 2.6870.7358 0.341623 2.687 2 | 1.634058 1.55 2 | 1.634058 1.55 0.2130 0.748525 3.5672100.2130 0.748525 3.567210 3 | 3.548077 13.40 3 | 3.548077 13.40 0.0003 1.719470 7.3230.0003 1.719470 7.323-Test of homogeneity (equal odds): chi2(3) Test of homogeneity (equal odds): chi2(3) = 17.5

8、2= 17.52 Prchi2 Prchi2 = 0.0006= 0.0006Score test for trend of odds: chi2(1) Score test for trend of odds: chi2(1) = 15.14= 15.14 Prchi2 Prchi2 = 0.0001= 0.0001例3 (分層病例對照資料) 在吸煙與肺癌發(fā)生關(guān)系的研討中，年齡是一個(gè)混雜要素，試根據(jù)以下資料計(jì)算調(diào)整年齡后吸煙與肺癌發(fā)生關(guān)系的優(yōu)勢比,并作假設(shè)檢驗(yàn)。 分層病例對照研討o目的：按混雜要素分層，控制混雜要素的影響。目的：按混雜要素分層，控制混雜要素的影響。o方法：方法：Mantel-

9、Haenszel分層分析分層分析 o(1)公共優(yōu)勢比公共優(yōu)勢比ORM-H的估計(jì)算公式為：的估計(jì)算公式為：o(2)對公共優(yōu)勢比進(jìn)展假設(shè)檢驗(yàn)對公共優(yōu)勢比進(jìn)展假設(shè)檢驗(yàn), 統(tǒng)計(jì)量為：統(tǒng)計(jì)量為：iiiiiiHMncbndaOR1)(22vVTaiiiHMiiiinbaT iiiiiiinnmmnnV30101分層病例對照分析occ 病例變量 暴露變量 fw=頻數(shù)變量，by(分層變量) omhodds 病例變量 暴露變量 fw=頻數(shù)變量 ，by(分層變量)數(shù)據(jù)構(gòu)造cc case smoke fw=f,by(age) age | OR 95% Conf. Interval M-H Weight-+- 0 |

10、 2.224658 1.437382 3.430672 14.28571 (exact) 1 | 4.174907 2.554659 6.824 7.819588 (exact)-+- Crude | 2.747456 2.011688 3.743611 (exact) M-H combined | 2.914544 2.143308 3.963297 -Test of homogeneity (M-H) chi2(1) = 3.94 Prchi2 = 0.0470 Test that combined OR = 1: Mantel-Haenszel chi2(1) = 49.43 Prchi

11、2 = 0.0000結(jié)果解釋o粗的OR2.75；調(diào)整年齡后，ORMH2.91。2MN=49.43，v=1，P=0.0000，差別有統(tǒng)計(jì)學(xué)意義，故可以為調(diào)整年齡后，吸煙1包/天以上者患肺癌的危險(xiǎn)度是1包/天以下者的2.91倍。例4(配對病例-對照研討) 某單位為研討軟組織肉瘤與接觸苯氧乙酸或氯酚的關(guān)系，作了一次病例-對照研討，結(jié)果見下表。 配對四格表資料o方法：方法：McNemar檢驗(yàn)檢驗(yàn)o公式：公式：)()(22cbcb22(1)()Cbcbcb+c 40b+c40配對病例對照分析omcc 暴露-病例數(shù) 暴露-對照數(shù) fw=頻數(shù)變量omcci #a #b #c #dmcci 3 16 3 30

12、mcci 3 16 3 30 | Controls | Controls | |Cases | Exposed Cases | Exposed Unexposed | TotalUnexposed | Total-+-+-+-+- Exposed | 3 Exposed | 3 16 | 1916 | 19 Unexposed | 3 Unexposed | 3 30 | 3330 | 33-+-+-+-+- Total | 6 Total | 6 46 | 5246 | 52McNemars chi2(1) = 8.89 Prob McNemars chi2(1) = 8.89 Prob c

13、hi2 = 0.0029chi2 = 0.0029Exact McNemar significance probability Exact McNemar significance probability = 0.0044= 0.0044Proportion with factorProportion with factor Cases .3653846 Cases .3653846 Controls .1153846 95% Controls .1153846 95% Conf. IntervalConf. Interval - - - - difference .25 .0811853 .

14、difference .25 .0811853 .41881474188147 ratio 3.166667 ratio 3.166667 1.422659 7.0486171.422659 7.048617 rel. rel. diff. .2826087 .1253028 .4399diff. .2826087 .1253028 .4399146146 odds ratio 5.333333 odds ratio 5.333333 1.526627 28.56284 (exact)1.526627 28.56284 (exact)b+c chi2 = 0.0281 chi2(1) = 4.

15、82 Prchi2 = 0.0281解釋：率差RD= 0.3623， 95% CI: 0.059, 0.666) 相對危險(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 ，可以為母親乳汁中IgG抗體滴度高低與嬰兒呼吸道疾患發(fā)生率有關(guān)，滴度較低者其嬰兒較易患呼吸道疾病。隊(duì)列研討資料-聯(lián)絡(luò)強(qiáng)度目的o率差率差RD=P1-P0o相對危險(xiǎn)度相對危險(xiǎn)度RR =P1/P0o歸因危險(xiǎn)度歸因危險(xiǎn)度 AR=(P1- P0 )/P1o人群

16、歸因危險(xiǎn)度人群歸因危險(xiǎn)度PAR =(P - P0 )/P例6分層隊(duì)列研討某人隨訪研討甲糖寧治療糖尿病人的效果，其中年齡被以為是能夠的混雜要素，需求對年齡進(jìn)展分層分析。分層隊(duì)列研討分析ocs 病例變量 暴露變量 fw=頻數(shù)變量，by(分層變量) 數(shù)據(jù)構(gòu)造cs d expose fw=f,by(age) age | RR 95% Conf. Interval M-H Weight-+- 0 | 1.811321 .6112044 5.367898 2.345133 1 | 1.332908 .7471514 2.37789 8.124352 -+- Crude | 1.505602 .891792

17、4 2.54189 M-H combined | 1.440071 .8637914 2.400817-Test of homogeneity (M-H) chi2(1) = 0.240 Prchi2 = 0.6244o粗的RR1.51；調(diào)整年齡后，RRMH1.44，其95%CI：0.86-2.40，差別無統(tǒng)計(jì)學(xué)意義，故調(diào)整年齡后，還不能以為甲糖寧治療糖尿病人有效。小結(jié)o病例對照研討o 成組病例對照：o cc case expose fw=f/cci #a #b #c #do cc case expose fw=f，by混雜變量分層分析o tabodds case expose fw=f，or

18、 base(#)多暴露程度 o 配對病例對照：mcc r c fw=f/mcci #a #b #c #do隊(duì)列研討o cs case expose fw=f /csi #a #b #c #do cs case expose fw=f，by混雜變量分層分析課堂練習(xí)課堂練習(xí)1：課堂練習(xí)課堂練習(xí)2：例：例18-5 在食管癌與飲酒關(guān)系的分析在食管癌與飲酒關(guān)系的分析中，年齡能夠是混雜要素。將每日飲酒量的資料按中，年齡能夠是混雜要素。將每日飲酒量的資料按年齡組分解成年齡組分解成4個(gè)四格表個(gè)四格表(表表18-3)，估計(jì)出各層的，估計(jì)出各層的OR值，可見不同年齡組的值，可見不同年齡組的OR有所不同。試進(jìn)展有所不同。試進(jìn)展Mantel-Haenszel分層分析。分層分析。 課堂練習(xí)課堂練習(xí)3：例：例18-6