流行病學(xué)資料的分析(ppt 42頁).ppt

1、STATA軟件的應(yīng)用,流行病學(xué)資料的分析,非常簡便、 快速,流行病學(xué)資料的分析,病例對(duì)照研究資料 成組病例對(duì)照研究 配對(duì)病例對(duì)照研究 隊(duì)列研究資料,例1(成組病例對(duì)照研究) 某單位研究胸膜間皮瘤與接觸石棉的關(guān)系，資料見下表。試對(duì)其進(jìn)行分析,成組病例對(duì)照研究分析,2檢驗(yàn)：接觸石棉與胸膜間皮瘤的關(guān)系 計(jì)算OR：反映接觸石棉與胸膜間皮瘤的關(guān)系大小 計(jì)算OR95%的置信區(qū)間：與2檢驗(yàn)等價(jià),成組病例對(duì)照分析,cc 病例變量 暴露變量 fw=頻數(shù)變量 cci #a #b #c #d,成組病例對(duì)照分析,cc case exposure fw=f cci 40 36 9 67,解釋：優(yōu)勢(shì)比OR = 8.27,

2、 其95% CI（3.42，21.33） 2 = 28.94 ， P = 0.0000 由此可見，接觸石棉者發(fā)生間皮瘤的危險(xiǎn)性為未接觸者的8.27倍,Proportion | Exposed Unexposed | Total Exposed -+-+- Cases | 40 36 | 76 0.5263 Controls | 9 67 | 76 0.1184 -+-+- Total | 49 103 | 152 0.3224 | | | Point estimate | 95% Conf. Interval |-+- Odds ratio | 8.271605 | 3.4193 21.330

3、91 (exact) Attr. frac. ex. | .8791045 | .7075425 .9531197 (exact) Attr. frac. pop | .4626866 | +- chi2(1) = 28.94 Prchi2 = 0.0000,病例對(duì)照研究資料-聯(lián)系強(qiáng)度指標(biāo),OR= ad/bc,例2（多暴露水平成組病例對(duì)照研究資料）某醫(yī)院一項(xiàng)吸煙與膀胱癌關(guān)系的病例對(duì)照研究得到如下結(jié)果，試分析吸煙量與膀胱癌是否存在劑量-反應(yīng)關(guān)系,多暴露水平成組病例對(duì)照研究,目的：多個(gè)暴露水平的劑量-反應(yīng)相關(guān)分析 方法：趨勢(shì)2檢驗(yàn) H0：吸煙量與膀胱癌不存在劑量-反應(yīng)線性關(guān)系 H1: 吸煙量與膀胱

4、癌存在劑量-反應(yīng)線性關(guān)系,多暴露水平成組病例對(duì)照分析,tabodds 病例變量 暴露變量 fw=頻數(shù)變量，or base(#) or：指定輸出各組的or值 base (#)：指定以哪組作為參照估計(jì)各組的or值，默認(rèn)是第一組作為參照,數(shù)據(jù)結(jié)構(gòu),解釋：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)性也增加,tabodds case smoke fw=f, or - smoke | Odds Ratio chi2 Pchi2 95% Conf. Interval -+- 0 | 1

5、.000000 . . . . 1 | 0.854167 0.11 0.7358 0.341623 2.135687 2 | 1.634058 1.55 0.2130 0.748525 3.567210 3 | 3.548077 13.40 0.0003 1.719470 7.321353 - Test of homogeneity (equal odds): chi2(3) = 17.52 Prchi2 = 0.0006 Score test for trend of odds: chi2(1) = 15.14 Prchi2 = 0.0001,例3 (分層病例對(duì)照資料) 在吸煙與肺癌發(fā)生關(guān)系

6、的研究中，年齡是一個(gè)混雜因素，試根據(jù)下列資料計(jì)算調(diào)整年齡后吸煙與肺癌發(fā)生關(guān)系的優(yōu)勢(shì)比,并作假設(shè)檢驗(yàn),分層病例對(duì)照研究,目的：按混雜因素分層，控制混雜因素的影響。 方法：Mantel-Haenszel分層分析 (1)公共優(yōu)勢(shì)比ORM-H的估計(jì)算公式為： (2)對(duì)公共優(yōu)勢(shì)比進(jìn)行假設(shè)檢驗(yàn), 統(tǒng)計(jì)量為,分層病例對(duì)照分析,cc 病例變量 暴露變量 fw=頻數(shù)變量，by(分層變量) mhodds 病例變量 暴露變量 fw=頻數(shù)變量 ，by(分層變量,數(shù)據(jù)結(jié)構(gòu),cc case smoke fw=f,by(age,age | OR 95% Conf. Interval M-H Weight -+- 0 | 2

7、.224658 1.437382 3.430672 14.28571 (exact) 1 | 4.174907 2.554659 6.821394 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 Pr

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

9、posed Unexposed | Total -+-+- Exposed | 3 16 | 19 Unexposed | 3 30 | 33 -+-+- Total | 6 46 | 52 McNemars chi2(1) = 8.89 Prob chi2 = 0.0029 Exact McNemar significance probability = 0.0044 Proportion with factor Cases .3653846 Controls .1153846 95% Conf. Interval - - difference .25 .0811853 .4188147 r

10、atio 3.166667 1.422659 7.048617 rel. diff. .2826087 .1253028 .4399146 odds ratio 5.333333 1.526627 28.56284 (exact,b+c 40,b+c40,結(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é)資料的分析,病例對(duì)照研究資料 成組病例對(duì)照研究 配對(duì)病例對(duì)照研究 隊(duì)列研究資料,例5（隊(duì)列研究資料） 就下表資料分析

11、母親乳汁中IgG抗體滴度高低與6個(gè)月以上嬰兒患呼吸道疾病的關(guān)系,隊(duì)列研究分析,2檢驗(yàn)：分析母親乳汁中IgG抗體滴度與呼吸道疾病的關(guān)系。 計(jì)算RR：分析母親乳汁中IgG抗體滴度與呼吸道疾病的關(guān)系大小。 計(jì)算RR95%的置信區(qū)間,隊(duì)列研究（Cohort studies）分析,cs 病例變量 暴露變量 fw=頻數(shù)變量 csi #a #b #c #d,隊(duì)列研究分析,cs case exposure fw=f csi 16 5 7 10, Exposed Unexposed | Total-+-+- Cases | 16 5 | 21 Noncases | 7 10 | 17-+-+- Total |

12、23 15 | 38 | | Risk | .6956522 .3333333 | .5526316 | Point estimate | 95% Conf. Interval |-+- Risk difference | .3623188 | .0585558 .6660819 Risk ratio | 2.086957 | .9711058 4.484977 Attr. frac. ex. | .5208333 | -.0297539 .7770334 Attr. frac. pop | .3968254 | +- chi2(1) = 4.82 Prchi2 = 0.0281,解釋：率差R

13、D= 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)，滴度較低者其嬰兒較易患呼吸道疾病,隊(duì)列研究資料-聯(lián)系強(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,例6（分層隊(duì)列研究）某人

14、隨訪研究甲糖寧治療糖尿病人的效果，其中年齡被認(rèn)為是可能的混雜因素，需要對(duì)年齡進(jìn)行分層分析,分層隊(duì)列研究分析,cs 病例變量 暴露變量 fw=頻數(shù)變量，by(分層變量,數(shù)據(jù)結(jié)構(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 .8917924 2.54189 M-H combined | 1.440071 .863

15、7914 2.400817 - Test of homogeneity (M-H) chi2(1) = 0.240 Prchi2 = 0.6244,粗的RR1.51；調(diào)整年齡后，RRMH1.44，其95%CI：0.86-2.40，差別無統(tǒng)計(jì)學(xué)意義，故調(diào)整年齡后，還不能認(rèn)為甲糖寧治療糖尿病人有效,小結(jié),病例對(duì)照研究 成組病例對(duì)照： cc case expose fw=f/cci #a #b #c #d cc case expose fw=f，by（混雜變量）分層分析 tabodds case expose fw=f，or base(#)多暴露水平 配對(duì)病例對(duì)照：mcc r c fw=f/mcci #a #b #c #d 隊(duì)列研究 cs case expose fw=f /csi #a #b #c #d cs case expose fw=f，by（混雜變量）分層分析,課堂練習(xí)1,課堂練習(xí)2：例18-5 在食管癌與飲酒關(guān)系的分析中，年齡可能是混雜因素。將每日飲酒量的資料按年齡組分解成4個(gè)四格表(