流行病學(xué)的分析

1、STATA軟件的應(yīng)用,--- 流行病學(xué)資料的分析,非常簡便、快速,流行病學(xué)資料的分析,病例對照研究資料 成組病例對照研究 配對病例對照研究 隊列研究資料,例1(成組病例對照研究) 某單位研究胸膜間皮瘤與接觸石棉的關(guān)系，資料見下表。試對其進行分析。,成組病例對照研究分析,?2檢驗：接觸石棉與胸膜間皮瘤的關(guān)系計算OR：反映接觸石棉與胸膜間皮瘤的關(guān)系大小計算OR95%的置信區(qū)間：與?2檢驗等價,成組病例對照分析

2、,cc 病例變量 暴露變量 [fw=頻數(shù)變量]cci #a #b #c #d,成組病例對照分析,cc case exposure [fw=f]cci 40 36 9 67,解釋：優(yōu)勢比OR = 8.27, 其95% CI（3.42，21.33） ?2 = 28.94 ， P = 0.0000 由此可見，接觸石棉者發(fā)生間皮瘤的危險性為未接觸者的8.27倍。

3、,Proportion | Exposed Unexposed | Total Exposed-----------------+------------------------+---------------------- Cases | 40 36 | 76 0.5263 Co

4、ntrols | 9 67 | 76 0.1184-----------------+------------------------+---------------------- Total | 49 103 | 152 0.3224 |

5、 | | Point estimate | [95% Conf. Interval] |------------------------+---------------------- Odds ratio | 8.271605 | 3.4193 21.33091 (exact)

6、Attr. frac. ex. | .8791045 | .7075425 .9531197 (exact) Attr. frac. pop | .4626866 | +----------------------------------------------- chi2(1)

7、 = 28.94 Pr>chi2 = 0.0000,病例對照研究資料---聯(lián)系強度指標(biāo),OR= ad/bc,例2（多暴露水平成組病例對照研究資料）某醫(yī)院一項吸煙與膀胱癌關(guān)系的病例對照研究得到如下結(jié)果，試分析吸煙量與膀胱癌是否存在劑量-反應(yīng)關(guān)系？,多暴露水平成組病例對照研究,目的：多個暴露水平的劑量-反應(yīng)相關(guān)分析 方法：趨勢?2檢驗 H0：吸煙量與膀胱癌不存在劑量-反應(yīng)線性關(guān)系 H1: 吸煙量

8、與膀胱癌存在劑量-反應(yīng)線性關(guān)系,,多暴露水平成組病例對照分析,tabodds 病例變量 暴露變量 [fw=頻數(shù)變量]，or base(#) or：指定輸出各組的or值 base (#)：指定以哪組作為參照估計各組的or值，默認(rèn)是第一組作為參照。,數(shù)據(jù)結(jié)構(gòu),解釋：OR的同質(zhì)性檢驗：?2 =17.52，P=0.0006，即各組OR不全同。 OR的趨勢檢驗：?2 =15.14，P=0.0001，故可用認(rèn)為隨

9、著吸煙量的增加患膀胱癌的危險性也增加。,tabodds case smoke [fw=f], or --------------------------------------------------------------------------- smoke | Odds Ratio chi2 P>chi2 [95% Conf. Interval]-------------

10、+------------------------------------------------------------- 0 | 1.000000 . . . . 1 | 0.854167 0.11 0.7358 0.341623 2.135687

11、 2 | 1.634058 1.55 0.2130 0.748525 3.567210 3 | 3.548077 13.40 0.0003 1.719470 7.321353-------------------------------------------------------------------------

12、--Test of homogeneity (equal odds): chi2(3) = 17.52 Pr>chi2 = 0.0006Score test for trend of odds: chi2(1) = 15.14 Pr>chi2 = 0

13、.0001,例3 (分層病例對照資料) 在吸煙與肺癌發(fā)生關(guān)系的研究中，年齡是一個混雜因素，試根據(jù)下列資料計算調(diào)整年齡后吸煙與肺癌發(fā)生關(guān)系的優(yōu)勢比,并作假設(shè)檢驗。,分層病例對照研究,目的：按混雜因素分層，控制混雜因素的影響。方法：Mantel-Haenszel分層分析 (1)公共優(yōu)勢比ORM-H的估計算公式為：(2)對公共優(yōu)勢比進行假設(shè)檢驗, 統(tǒng)計量為：,,,,,分層病例對照分析,cc 病例變量 暴露變量 [fw=頻數(shù)變

14、量]，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-----------------+-------------------------------------------------

15、 0 | 2.224658 1.437382 3.430672 14.28571 (exact) 1 | 4.174907 2.554659 6.821394 7.819588 (exact)-----------------+-------------------------------------------------

16、 Crude | 2.747456 2.011688 3.743611 (exact) M-H combined | 2.914544 2.143308 3.963297 -------------------------------------------------------------------Test of homogenei

17、ty (M-H) chi2(1) = 3.94 Pr>chi2 = 0.0470 Test that combined OR = 1: Mantel-Haenszel chi2(1) = 49.43 Pr&

18、gt;chi2 = 0.0000,結(jié)果解釋,粗的OR＝2.75；調(diào)整年齡后，ORMH＝2.91。?2MN=49.43，v=1，P=0.0000，差別有統(tǒng)計學(xué)意義，故可認(rèn)為調(diào)整年齡后，吸煙1包/天以上者患肺癌的危險度是1包/天以下者的2.91倍。,例4(配對病例-對照研究) 某單位為研究軟組織肉瘤與接觸苯氧乙酸或氯酚的關(guān)系，作了一次病例-對照研究，結(jié)果見下表。,配對四格表資料,方法：McNemar檢驗公式：,b+c ≤ 40,b

19、+c>40,配對病例對照分析,mcc 暴露-病例數(shù) 暴露-對照數(shù) [fw=頻數(shù)變量]mcci #a #b #c #d,mcci 3 16 3 30 | Controls |Cases | Exposed Unexposed | Total-----------------+---------------

20、---------+---------- Exposed | 3 16 | 19 Unexposed | 3 30 | 33-----------------+------------------------+---------- Total | 6

21、 46 | 52McNemar's chi2(1) = 8.89 Prob > chi2 = 0.0029Exact McNemar significance probability = 0.0044Proportion with factor Cases .3653846 Controls .1153846 [

22、95% Conf. Interval] --------- -------------------- difference .25 .0811853 .4188147 ratio 3.166667 1.422659 7.048617 rel. diff. .2826087 .125

23、3028 .4399146 odds ratio 5.333333 1.526627 28.56284 (exact),b+c < 40,b+c≥40,結(jié)果解釋,OR=5.33，其95% CI為(1.53, 28.56)。 本例樣本例數(shù)較少，應(yīng)計算確切概率。P = 0.0044，差別有統(tǒng)計學(xué)意義，故可認(rèn)為接觸苯氧乙酸或氯酚者發(fā)生軟組織肉瘤的危險性為不接觸者的5.33倍。,流行病

24、學(xué)資料的分析,病例對照研究資料 成組病例對照研究 配對病例對照研究 隊列研究資料,例5（隊列研究資料） 就下表資料分析母親乳汁中IgG抗體滴度高低與6個月以上嬰兒患呼吸道疾病的關(guān)系。,隊列研究分析,?2檢驗：分析母親乳汁中IgG抗體滴度與呼吸道疾病的關(guān)系。計算RR：分析母親乳汁中IgG抗體滴度與呼吸道疾病的關(guān)系大小。計算RR95%的置信區(qū)間,隊列研究（Cohort studies）分析,cs 病例變量

25、暴露變量 [fw=頻數(shù)變量]csi #a #b #c #d,隊列研究分析,cs case exposure [fw=f]csi 16 5 7 10,| Exposed Unexposed | Total-----------------+------------------------+---------- Cases | 16

26、5 | 21 Noncases | 7 10 | 17-----------------+------------------------+---------- Total | 23 15 | 38 | |

27、 Risk | .6956522 .3333333 | .5526316 | Point estimate | [95% Conf. Interval] |------------------------+---------------------- Risk difference | .36231

28、88 | .0585558 .6660819 Risk ratio | 2.086957 | .9711058 4.484977 Attr. frac. ex. | .5208333 | -.0297539 .7770334 Attr. frac. pop | .3968254 |

29、 +----------------------------------------------- chi2(1) = 4.82 Pr>chi2 = 0.0281,解釋：率差RD= 0.3623， 95% CI: 0.059, 0.666) 相對危險度RR = 2.086， 95% CI: (0.971,

30、4.485) 歸因危險度 AR = 0.521， 其95% CI（-0.030 ，0.777） 人群歸因危險度PAR = 0.397 ?2 = 4.82，P = 0.0281 ，可認(rèn)為母親乳汁中IgG抗體滴度高低與嬰兒呼吸道疾患發(fā)生率有關(guān)，滴度較低者其嬰兒較易患呼吸道疾病。,隊列研究資料---聯(lián)系強度指標(biāo),率差RD=P1-P0相對危險度RR =P1/P0

31、歸因危險度 AR=(P1- P0 )/P1人群歸因危險度PAR =(P - P0 )/P,例6（分層隊列研究）某人隨訪研究甲糖寧治療糖尿病人的效果，其中年齡被認(rèn)為是可能的混雜因素，需要對年齡進行分層分析。,分層隊列研究分析,cs 病例變量 暴露變量 [fw=頻數(shù)變量]，by(分層變量),數(shù)據(jù)結(jié)構(gòu),cs d expose [fw=f],by(age),age | RR [95% Conf. Inter

32、val] M-H Weight-----------------+------------------------------------------------- 0 | 1.811321 .6112044 5.367898 2.345133 1 | 1.332908 .7471514 2.37789 8.124

33、352 -----------------+------------------------------------------------- Crude | 1.505602 .8917924 2.54189 M-H combined | 1.440071 .8637914 2.400817------------------------

34、-------------------------------------------Test of homogeneity (M-H) chi2(1) = 0.240 Pr>chi2 = 0.6244,粗的RR＝1.51；調(diào)整年齡后，RRMH＝1.44，其95%CI：0.86-2.40，差別無統(tǒng)計學(xué)意義，故調(diào)整年齡后，還不能認(rèn)為甲糖寧治療糖尿病人有效。,小結(jié),病例對照研究 成組病例對照： cc

35、 case expose [fw=f]/cci #a #b #c #d cc case expose [fw=f]，by（混雜變量）分層分析 tabodds case expose [fw=f]，or base(#)多暴露水平 配對病例對照：mcc r c [fw=f]/mcci #a #b #c #d隊列研究 cs case expose [fw=f] /csi #a #b

36、 #c #d cs case expose [fw=f]，by（混雜變量）分層分析,課堂練習(xí)1：,課堂練習(xí)2：例18-5 在食管癌與飲酒關(guān)系的分析中，年齡可能是混雜因素。將每日飲酒量的資料按年齡組分解成4個四格表(表18-3)，估計出各層的OR值，可見不同年齡組的OR有所不同。試進行Mantel-Haenszel分層分析。,課堂練習(xí)3：例18-6 仍然用食管癌與飲酒關(guān)系的研究資料，按照每日飲酒量將暴露水平分成四組，資料