《[精选]流行病学资料的分析》由会员分享,可在线阅读,更多相关《[精选]流行病学资料的分析(42页珍藏版)》请在金锄头文库上搜索。
1、STATA软件的应用,- 流行病学资料的分析,非常简便、 快速,流行病学资料的分析,病例对照研究资料 成组病例对照研究 配对病例对照研究 队列研究资料,例1(成组病例对照研究) 某单位研究胸膜间皮瘤与接触石棉的关系,资料见下表。试对其进行分析。,成组病例对照研究分析,2检验:接触石棉与胸膜间皮瘤的关系 计算OR:反映接触石棉与胸膜间皮瘤的关系大小 计算OR95%的置信区间:与2检验等价,成组病例对照分析,cc 病例变量 暴露变量 fw=频数变量 cci #a #b #c #d,成组病例对照分析,cc case exposure fw=f cci 40 36 9 67,解释:优势比OR = 8.
2、27, 其95% CI(3.42,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
3、.33091 (exact) Attr. 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 -+
5、- 0 | 1.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 (分层病例对照资料) 在吸烟
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 Weigh
7、t -+- 0 | 2.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)
8、 = 49.43 Prchi2 = 0.0000,结果解释,粗的OR2.75;调整年龄后,ORMH2.91。2MN=49.43,v=1,P=0.0000,差别有统计学意义,故可认为调整年龄后,吸烟1包/天以上者患肺癌的危险度是1包/天以下者的2.91倍。,例4(配对病例-对照研究) 某单位为研究软组织肉瘤与接触苯氧乙酸或氯酚的关系,作了一次病例-对照研究,结果见下表。,配对四格表资料,方法:McNemar检验 公式:,b+c 40,b+c40,配对病例对照分析,mcc 暴露-病例数 暴露-对照数 fw=频数变量 mcci #a #b #c #d,mcci 3 16 3 30 | Control
9、s | Cases | Exposed 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 .0811
10、853 .4188147 ratio 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,结果解释,OR=5.33,其95% CI为(1.53, 28.56)。 本例样本例数较少,应计算确切概率。P = 0.0044,差别有统计学意义,故可认为接触苯氧乙酸或氯酚者发生软组织肉瘤的危险性为不接触者的5.33倍。,流行病学资料的分析,病例对照研究资料 成组病例对照研究 配对病例对照研究 队列研究资料,例5
11、(队列研究资料) 就下表资料分析母亲乳汁中IgG抗体滴度高低与6个月以上婴儿患呼吸道疾病的关系。,队列研究分析,2检验:分析母亲乳汁中IgG抗体滴度与呼吸道疾病的关系。 计算RR:分析母亲乳汁中IgG抗体滴度与呼吸道疾病的关系大小。 计算RR95%的置信区间,队列研究(Cohort studies)分析,cs 病例变量 暴露变量 fw=频数变量 csi #a #b #c #d,队列研究分析,cs case exposure fw=f csi 16 5 7 10,| Exposed Unexposed | Total-+-+- Cases | 16 5 | 21 Noncases | 7 10
12、| 17-+-+- Total | 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 Prch
13、i2 = 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抗体滴度高低与婴儿呼吸道疾患发生率有关,滴度较低者其婴儿较易患呼吸道疾病。,队列研究资料-联系强度指标,率差RD=P1-P0 相对危险度RR =P1/P0 归因危险度 AR=(P1- P0 )/P1 人群归因危险度PAR =(P -
14、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 .7471514 2.37789 8.124352 -+- Crude | 1.505602 .8917924 2.54189 M-H co
15、mbined | 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(#)多暴露水平 配
16、对病例对照: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有所不同。试进行Mantel-Haenszel分层分析。,课堂练习3:例18-6 仍然用食管癌与饮酒关系的研究资料,按照每日饮酒量将暴露水平分成四组,资料列在表18-4,试分析每日饮酒量与食管癌发病是否存在剂量-反应关系。,课堂练习4:例18-7 利用子宫内膜癌与服用雌激素病例-对照研究的资料,资料整理成表18-6。试对以下资料进行分析。,课后练习,利用帮助功能学习STATA软件中的病例对照研究资料的分析 help cc help cs,