流行病学的分析课件

1、STATA软件的应用- 流行病学资料的分析非常简便、快速流行病学资料的分析病例对照研究资料 成组病例对照研究 配对病例对照研究 队列研究资料 例1(成组病例对照研究) 某单位研究胸膜间皮瘤与接触石棉的关系，资料见下表。试对其进行分析。 成组病例对照研究分析2检验：接触石棉与胸膜间皮瘤的关系计算OR：反映接触石棉与胸膜间皮瘤的关系大小计算OR95%的置信区间：与2检验等价成组病例对照分析cc 病例变量 暴露变量 fw=频数变量cci #a #b #c #d成组病例对照分析cc case exposure fw=fcci 40 36 9 67解释：优势比OR = 8.27, 其95% CI（3.4

2、2，21.33） 2 = 28.94 ， P = 0.0000 由此可见，接触石棉者发生间皮瘤的危险性为未接触者的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.33091 (exact) At

3、tr. frac. ex. | .8791045 | .7075425 .9531197 (exact) Attr. frac. pop | .4626866 | +- chi2(1) = 28.94 Prchi2 = 0.0000病例对照研究资料-联系强度指标OR= ad/bc例2（多暴露水平成组病例对照研究资料）某医院一项吸烟与膀胱癌关系的病例对照研究得到如下结果，试分析吸烟量与膀胱癌是否存在剂量-反应关系？多暴露水平成组病例对照研究目的：多个暴露水平的剂量-反应相关分析 方法：趋势2检验 H0：吸烟量与膀胱癌不存在剂量-反应线性关系 H1: 吸烟量与膀胱癌存在剂量-反应线性关系多暴露水平

4、成组病例对照分析tabodds 病例变量 暴露变量 fw=频数变量，or base(#) or：指定输出各组的or值 base (#)：指定以哪组作为参照估计各组的or值，默认是第一组作为参照。数据结构解释：OR的同质性检验：2 =17.52，P=0.0006，即各组OR不全同。 OR的趋势检验：2 =15.14，P=0.0001，故可用认为随着吸烟量的增加患膀胱癌的危险性也增加。 tabodds case smoke fw=f, or - smoke | Odds Ratio chi2 Pchi2 95% Conf. Interval-+- 0 | 1.000000 . . . . 1 |

5、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.0006Score test for trend of odds: chi2(1) = 15.14 Prchi2 = 0.0001例3 (分层病例对照资料) 在吸烟与肺癌发生关系的研究中，年龄是一个混杂因素，试根据下列资料计算

6、调整年龄后吸烟与肺癌发生关系的优势比,并作假设检验。 分层病例对照研究目的：按混杂因素分层，控制混杂因素的影响。方法：Mantel-Haenszel分层分析 (1)公共优势比ORM-H的估计算公式为：(2)对公共优势比进行假设检验, 统计量为：分层病例对照分析cc 病例变量 暴露变量 fw=频数变量，by(分层变量) mhodds 病例变量 暴露变量 fw=频数变量 ，by(分层变量)数据结构cc case smoke fw=f,by(age) age | OR 95% Conf. Interval M-H Weight-+- 0 | 2.224658 1.437382 3.430672 14

7、.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 Prchi2 = 0.0000结果解释粗的OR2.75；调整年龄

8、后，ORMH2.91。2MN=49.43，v=1，P=0.0000，差别有统计学意义，故可认为调整年龄后，吸烟1包/天以上者患肺癌的危险度是1包/天以下者的2.91倍。例4(配对病例-对照研究) 某单位为研究软组织肉瘤与接触苯氧乙酸或氯酚的关系，作了一次病例-对照研究，结果见下表。 配对四格表资料方法：McNemar检验公式：b+c 40b+c40配对病例对照分析mcc 暴露-病例数 暴露-对照数 fw=频数变量mcci #a #b #c #dmcci 3 16 3 30 | Controls |Cases | Exposed Unexposed | Total-+-+- Exposed |

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

10、2826087 .1253028 .4399146 odds ratio 5.333333 1.526627 28.56284 (exact)b+c 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抗体滴度高低与婴儿呼吸道疾患发生率有关，滴度较低者其婴儿较易患呼吸道疾病。队列研究

11、资料-联系强度指标率差RD=P1-P0相对危险度RR =P1/P0归因危险度 AR=(P1- P0 )/P1人群归因危险度PAR =(P - P0 )/P例6（分层队列研究）某人随访研究甲糖宁治疗糖尿病人的效果，其中年龄被认为是可能的混杂因素，需要对年龄进行分层分析。分层队列研究分析cs 病例变量 暴露变量 fw=频数变量，by(分层变量) 数据结构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 .747

12、1514 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.6244粗的RR1.51；调整年龄后，RRMH1.44，其95%CI：0.86-2.40，差别无统计学意义，故调整年龄后，还不能认为甲糖宁治疗糖尿病人有效。小结病例对照研究 成组病例对照： cc 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 #c #d cs case expose fw=f，by（混杂变量）分层分析课堂练习1：课堂练习2：例18-5 在食管癌与饮酒关系的分析中，年龄可能是混杂因素。将每日饮酒量的资料按年龄组分解成4个四格表(表18-3)，估计出各层的OR值，可见不同年龄组的OR有所不同。试进行Mante