流行病学资料的分析课件

1、流行病学资料的分析1stata软件的应用- - 流行病学资料的分析流行病学资料的分析非常简便、非常简便、快速快速流行病学资料的分析2流行病学资料的分析流行病学资料的分析病例对照研究资料 成组病例对照研究成组病例对照研究 配对病例对照研究配对病例对照研究 队列研究资料流行病学资料的分析3 例例1(成组病例对照研究成组病例对照研究) 某单位研究胸膜间皮瘤与接触某单位研究胸膜间皮瘤与接触石棉的关系，资料见下表。试对其进行分析。石棉的关系，资料见下表。试对其进行分析。 流行病学资料的分析4成组病例对照研究分析检验：接触石棉与胸膜间皮瘤的关系计算or：反映接触石棉与胸膜间皮瘤的关系大小计算or95%的置

2、信区间：与检验等价流行病学资料的分析5成组病例对照分析cc 病例变量病例变量 暴露变量暴露变量 fw=fw=频数变量频数变量 cci #a #b #c #d流行病学资料的分析6成组病例对照分析cc case exposure fw=fcci 40 36 9 67流行病学资料的分析7解释：优势比优势比or = 8.27, 其其95% ci（3.42，21.33） 2 = 28.94 ， p = 0.0000 由此可见，接触石棉者发生间皮瘤的危险性为未接触者的由此可见，接触石棉者发生间皮瘤的危险性为未接触者的8.278.27倍。倍。 proportion proportion | exposed

3、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 | 49 103 | 152 0.3224 total | 49 103 | 152 0.3224 | | | | | point estimate | 95% conf. interval | point es

4、timate | 95% conf. interval |-+- |-+- odds ratio | 8.271605 | 3.4193 21.33091odds ratio | 8.271605 | 3.4193 21.33091 (exact) (exact) attr. frac. ex. | .8791045 | .7075425 .9531197 (exact) attr. frac. ex. | .8791045 | .7075425 .9531197 (exact) attr. frac. pop | .4626866 | attr. frac. pop | .4626866 |

5、 +- +- chi2(1) = 28.94 prchi2 = 0.0000chi2(1) = 28.94 prchi2 = 0.0000流行病学资料的分析8病例对照研究资料-联系强度指标or= ad/bc流行病学资料的分析9例2（多暴露水平成组病例对照研究资料）某医院一项吸烟与膀胱癌关系的病例对照研究得到如下结果，试分析吸烟量与膀胱癌是否存在剂量-反应关系？流行病学资料的分析10多暴露水平成组病例对照研究目的：多个暴露水平的剂量目的：多个暴露水平的剂量-反应相关分析反应相关分析 方法：方法：趋势趋势检验检验 h h0 0：吸烟量与膀胱癌不存在剂量：吸烟量与膀胱癌不存在剂量- -反应线性关系反

6、应线性关系 h h1 1: : 吸烟量与膀胱癌存在剂量吸烟量与膀胱癌存在剂量- -反应线性关系反应线性关系1)()()(0202013022vxmxmnnnnntaxkikiiiiikiiii流行病学资料的分析11多暴露水平成组病例对照分析otabodds 病例变量 暴露变量 fw=频数变量，or base(#) or：指定输出各组的or值 base (#)：指定以哪组作为参照估计各组的or值，默认是第一组作为参照。流行病学资料的分析12数据结构流行病学资料的分析13解释：or的同质性检验：的同质性检验： =17.52，p=0.0006，即各组，即各组or不全同。不全同。 or的趋势检验：的趋

7、势检验： =15.14，p=0.0001，故可用认为随着吸烟量的，故可用认为随着吸烟量的增加患膀胱癌的危险性也增加。增加患膀胱癌的危险性也增加。 tabodds case smoke fw=f, or tabodds case smoke fw=f, or - smoke | odds ratio chi2 pchi2 95% conf. interval smoke | odds ratio chi2 pchi2 95% conf. interval-+-+- 0 | 1.000000 . . . . 0 | 1.000000 . . . . 1 | 0.854167 0.11 0.7358

8、 0.341623 2.135687 1 | 0.854167 0.11 0.7358 0.341623 2.135687 2 | 1.634058 1.55 0.2130 0.748525 3.567210 2 | 1.634058 1.55 0.2130 0.748525 3.567210 3 | 3.548077 13.40 0.0003 1.719470 7.321353 3 | 3.548077 13.40 0.0003 1.719470 7.321353-test of homogeneity (equal odds):test of homogeneity (equal odds

9、): chi2(3) = 17.52 chi2(3) = 17.52 prchi2 = 0.0006 prchi2 = 0.0006score test for trend of odds:score test for trend of odds: chi2(1) = 15.14 chi2(1) = 15.14 prchi2 = 0.0001 prchi2 = 0.0001流行病学资料的分析14例3 (分层病例对照资料) 在吸烟与肺癌发生关系的研究中，年龄是一个混杂因素，试根据下列资料计算调整年龄后吸烟与肺癌发生关系的优势比,并作假设检验。 流行病学资料的分析15分层病例对照研究目的：按混杂因

10、素分层，控制混杂因素的影响。目的：按混杂因素分层，控制混杂因素的影响。方法方法：mantel-haenszel分层分析分层分析 (1)(1)公共优势比公共优势比ororm-hm-h的估计算公式为：的估计算公式为：(2)(2)对公共优势比进行假设检验对公共优势比进行假设检验, , 统计量为：统计量为：iiiiiihmncbndaor1)(22vvtaiiihmiiiinbat iiiiiiinnmmnnv30101流行病学资料的分析16分层病例对照分析occ 病例变量 暴露变量 fw=频数变量，by(分层变量) omhodds 病例变量 暴露变量 fw=频数变量 ，by(分层变量)流行病学资料的

11、分析17数据结构流行病学资料的分析18cc case smoke fw=f,by(age) age | or 95% conf. interval m-h weightage | or 95% conf. interval m-h weight-+-+- 0 | 2.224658 1.437382 3.430672 14.28571 (exact) 0 | 2.224658 1.437382 3.430672 14.28571 (exact) 1 | 4.174907 2.554659 6.821394 7.819588 (exact) 1 | 4.174907 2.554659 6.8213

12、94 7.819588 (exact)-+-+- crude | crude | 2.7474562.747456 2.011688 3.743611 (exact) 2.011688 3.743611 (exact) m-h combined | m-h combined | 2.9145442.914544 2.143308 3.963297 2.143308 3.963297 -test of homogeneity (m-h) chi2(1) = 3.94 prchi2 = 0.0470test of homogeneity (m-h) chi2(1) = 3.94 prchi2 =

13、0.0470 test that combined or = 1: test that combined or = 1: mantel-haenszel chi2(1)mantel-haenszel chi2(1) = 49.43 = 49.43 prchi2 prchi2 = 0.0000 = 0.0000流行病学资料的分析19结果解释粗的or2.75；调整年龄后，ormh2.91。2mn=49.43，v=1，p=0.0000，差别有统计学意义，故可认为调整年龄后，吸烟1包/天以上者患肺癌的危险度是1包/天以下者的2.91倍。流行病学资料的分析20例4(配对病例-对照研究) 某单位为研究软组

14、织肉瘤与接触苯氧乙酸或氯酚的关系，作了一次病例-对照研究，结果见下表。 流行病学资料的分析21配对四格表资料方法：方法：mcnemar检验检验公式：公式：)()(22cbcb22(1)()cbcbcb+c 40b+c40流行病学资料的分析22配对病例对照分析mcc 暴露暴露- -病例数病例数 暴露暴露- -对照数对照数 fw=fw=频数变量频数变量 mcci #a #b #c #d流行病学资料的分析23mcci 3 16 3 30mcci 3 16 3 30 | controls | | controls |cases | exposed unexposed | totalcases | ex

15、posed unexposed | total-+-+-+-+- exposed | 3 16 | 19 exposed | 3 16 | 19 unexposed | 3 30 | 33 unexposed | 3 30 | 33-+-+-+-+- total | 6 46 | 52 total | 6 46 | 52mcnemars chi2(1)mcnemars chi2(1) = 8.89 prob chi2 = 0.0029 = 8.89 prob chi2 = 0.0029exact mcnemarexact mcnemar significance probability = 0

16、.0044 significance probability = 0.0044proportion with factorproportion with factor cases .3653846 cases .3653846 controls .1153846 95% conf. interval controls .1153846 95% conf. interval - - - - difference .25 .0811853 .4188147 difference .25 .0811853 .4188147 ratio 3.166667 1.422659 7.048617 ratio

17、 3.166667 1.422659 7.048617 rel. diff. .2826087 .1253028 .4399146 rel. diff. .2826087 .1253028 .4399146 odds ratio 5.333333 1.526627 28.56284 (exact) odds ratio 5.333333 1.526627 28.56284 (exact)b+c chi2 = 0.0281 chi2(1) = 4.82 prchi2 = 0.0281解释：率差率差rd= 0.3623， 95% ci: 0.059, 0.666) 相对危险度相对危险度rr = 2

18、.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 ，可认为母亲乳汁中可认为母亲乳汁中iggigg抗体滴度高低与婴儿呼抗体滴度高低与婴儿呼吸道疾患发生率有关，吸道疾患发生率有关，滴度较低滴度较低者其婴儿较易患呼吸道疾病者其婴儿较易患呼吸道疾病。流行病学资料的分析31队列研究资料-联系强度指标率差率差rd=p1-p0相对危险度相对危险度rr =p1/p0归因危险度归因危险度 ar=(p1- p0 )/p1人

19、群归因危险度人群归因危险度par =(p - p0 )/p流行病学资料的分析32例6（分层队列研究）某人随访研究甲糖宁治疗糖尿病人的效果，其中年龄被认为是可能的混杂因素，需要对年龄进行分层分析。流行病学资料的分析33分层队列研究分析ocs 病例变量 暴露变量 fw=频数变量，by(分层变量) 流行病学资料的分析34数据结构流行病学资料的分析35cs 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 .747151

20、4 2.37789 8.124352 -+- crude | 1.505602 .8917924 2.54189 m-h combined | 1.440071 .8637914 2.400817-test of homogeneity (m-h) chi2(1) = 0.240 prchi2 = 0.6244o粗的rr1.51；调整年龄后，rrmh1.44，其95%ci：0.86-2.40，差别无统计学意义，故调整年龄后，还不能认为甲糖宁治疗糖尿病人有效。流行病学资料的分析36小结病例对照研究 成组病例对照： cc case expose fw=f/cci #a #b #c #dcase expose fw=f/cci #a #b #c #d cc case expose fw=fcase expose fw=f，byby（混杂变量）（混杂变量）分层分析分层分析 tabodds case expose case expose fw=f，or base(#)多暴露水平 配对病例对照：mcc r r c fw=f/mc fw=