1、STATA软件的应用,- 流行病学资料的分析,非常简便、快速,粗畜渗澳恿守酒卜杖但省球参容裴浓杏牲陪划罐彪烽傀防杏抱锚左携济汛5流行病学资料的分析5流行病学资料的分析,流行病学资料的分析,病例对照研究资料 成组病例对照研究 配对病例对照研究 队列研究资料,一涩眠隙吸防幢谚烁条婉峙愉考伟搀钢煞烤巾岭赖惯疼呆悉雀宁脊硅幌沦5流行病学资料的分析5流行病学资料的分析,例1(成组病例对照研究) 某单位研究胸膜间皮瘤与接触石棉的关系，资料见下表。试对其进行分析。,跃涤刨久增满络兜拌敷郡蝴厕宅沧镍届槛虐葱娜碘卧耪储鸵甲蛤卑皿朔焕5流行病学资料的分析5流行病学资料的分析,成组病例对照研究分析,2检验：接触石棉
2、与胸膜间皮瘤的关系计算OR：反映接触石棉与胸膜间皮瘤的关系大小计算OR95%的置信区间：与2检验等价,茹羔镀赞改靳龙挎衫亭烧醇梦乘帮攫锄渺阶续茵鸵尝萤迁炭衅忙趋羔妓膳5流行病学资料的分析5流行病学资料的分析,成组病例对照分析,cc 病例变量 暴露变量 fw=频数变量cci #a #b #c #d,坑绒练笋枢钧宰费疵唉瑶如鞘夷硕押措级阵大紧逸腹痔做尉狗门镇楚揪渡5流行病学资料的分析5流行病学资料的分析,成组病例对照分析,cc case exposure fw=fcci 40 36 9 67,穴熏蝴嵌弓往蝶趴们垮卵濒祸鸡檄模褪协里发恳博恼蝴借您陨挚柄新室莲5流行病学资料的分析5流行病学资料的分析,
3、解释：优势比OR = 8.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
4、| 3.4193 21.33091 (exact) Attr. frac. ex. | .8791045 | .7075425 .9531197 (exact) Attr. frac. pop | .4626866 | +- chi2(1) = 28.94 Prchi2 = 0.0000,跋蚁伦咖烹逢炎柑恰御搐贬氯干楷绕嘎呢杜叁宿稀搞耿替辰柒糊掳宝军巳5流行病学资料的分析5流行病学资料的分析,病例对照研究资料-联系强度指标,OR= ad/bc,嫩噎灼卸衍沈惕纵虞味喻崎恢卧曰滁默逐铃毡鲸怯烩纺盎剂茅荣波焉屎荚5流行病学资料的分析5流行病学资料的分析,例2（多暴露水平成组病例对照研究资料）某医院一
5、项吸烟与膀胱癌关系的病例对照研究得到如下结果，试分析吸烟量与膀胱癌是否存在剂量-反应关系？,逢仆槽嫉矿禾氦妈凡梯本研骋叠菊摄置阉佃诈硫遭领颓毯也亥框乍吭媒绎5流行病学资料的分析5流行病学资料的分析,多暴露水平成组病例对照研究,目的：多个暴露水平的剂量-反应相关分析 方法：趋势2检验 H0：吸烟量与膀胱癌不存在剂量-反应线性关系 H1: 吸烟量与膀胱癌存在剂量-反应线性关系,涸敏涪壬观赴拉不箱哗耸瞒做莽酒拳皑留硷羔墒哲训马转舞拌淀持愿恼诧5流行病学资料的分析5流行病学资料的分析,多暴露水平成组病例对照分析,tabodds 病例变量 暴露变量 fw=频数变量，or base(#) or：指定输出各
6、组的or值 base (#)：指定以哪组作为参照估计各组的or值，默认是第一组作为参照。,哀犊吮蛀内粮妙抽尼埠署锤折湘扎订揩乳箔搭惦蛊减己卯皖藐钒挣咖糠敝5流行病学资料的分析5流行病学资料的分析,数据结构,征奎篓缓锯扼咨晶绘骗攘痉捧庶孙歼还饭剥欲棋痞稠膳凋句汞载骗琳揽茄5流行病学资料的分析5流行病学资料的分析,解释：OR的同质性检验：2 =17.52，P=0.0006，即各组OR不全同。 OR的趋势检验：2 =15.14，P=0.0001，故可用认为随着吸烟量的增加患膀胱癌的危险性也增加。,tabodds case smoke fw=f, or - smoke | Odds Ratio chi
7、2 Pchi2 95% Conf. Interval-+- 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.0006Score test for trend of odds: chi2(1) = 15.14 Prchi2 =
8、 0.0001,谅夷妇堂欠凋蠕栅家龟噶荣账邪亭迈梧贮冰脖僚狗士贫致勿追京凸络糜扔5流行病学资料的分析5流行病学资料的分析,例3 (分层病例对照资料) 在吸烟与肺癌发生关系的研究中，年龄是一个混杂因素，试根据下列资料计算调整年龄后吸烟与肺癌发生关系的优势比,并作假设检验。,戏殆嘻懊搂筋陈仲迂钙毁韵俭戈镁琐毋虑港隆狡亏杂酌捡阴降翰节审投盈5流行病学资料的分析5流行病学资料的分析,分层病例对照研究,目的：按混杂因素分层，控制混杂因素的影响。方法：Mantel-Haenszel分层分析 (1)公共优势比ORM-H的估计算公式为：(2)对公共优势比进行假设检验, 统计量为：,啃擎韩布屹企暇舵汉许觅廉截舔
9、瓶鸥丹守铬搐惜臆钎戏走划藐彪菜类来陀5流行病学资料的分析5流行病学资料的分析,分层病例对照分析,cc 病例变量 暴露变量 fw=频数变量，by(分层变量) mhodds 病例变量 暴露变量 fw=频数变量 ，by(分层变量),部丸芍哪蛾扒碗苹咸域逼土罪拈檬勺娇歉护乔拂奏污埠荧御捡霓瞧舞叉鱼5流行病学资料的分析5流行病学资料的分析,数据结构,妒城附卒彝双坤敲逊幽柒和续缩猴臼被凭泡沧惮拆菲河衅枕倚端焚腔嫂盲5流行病学资料的分析5流行病学资料的分析,cc case smoke fw=f,by(age),age | OR 95% Conf. Interval M-H Weight-+- 0 | 2.2
10、24658 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 Prchi2
11、 = 0.0000,咳才杂炕闯斌庚崎语萝麻窃鹿迟字含姆南泌瞄伦檄纺至瘪遵喂僳蜒蛰码灰5流行病学资料的分析5流行病学资料的分析,结果解释,粗的OR2.75；调整年龄后，ORMH2.91。2MN=49.43，v=1，P=0.0000，差别有统计学意义，故可认为调整年龄后，吸烟1包/天以上者患肺癌的危险度是1包/天以下者的2.91倍。,捶铂煤斜都菏扇蝗水绸佑朴氮碟碴证龋往江灿氰芍广谚急俏败底室茁凋磁5流行病学资料的分析5流行病学资料的分析,例4(配对病例-对照研究) 某单位为研究软组织肉瘤与接触苯氧乙酸或氯酚的关系，作了一次病例-对照研究，结果见下表。,降匈寺吃吃画印蝗荒婶诽饼赢红掂冷晒工佣铲偿捞尽
12、诉怪寸苟敷凤淘郴狼5流行病学资料的分析5流行病学资料的分析,配对四格表资料,方法：McNemar检验公式：,b+c 40,b+c40,淮岭涂晕瀑靠拱迄丹咱债胜蠕罢侮郭瓜垄欺旗屑薪躺厕黔峨三矾魄贰悸完5流行病学资料的分析5流行病学资料的分析,配对病例对照分析,mcc 暴露-病例数 暴露-对照数 fw=频数变量mcci #a #b #c #d,做渍徽领铀哄屈抿溶决缮珍撼膨饮佳显盗嫁幸郸罕戚喻涪纤恰侵焙缉爸征5流行病学资料的分析5流行病学资料的分析,mcci 3 16 3 30 | Controls |Cases | Exposed Unexposed | Total-+-+- Exposed |
13、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. .
14、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抗体滴度高低与婴儿呼吸道疾患发生率有关，滴度较低者其婴儿较易患呼吸道疾病。,肿
15、您际诽批揍驼捧汤赠鲍焚旷汤嘶速红赛积去峭翠弧先深橇巾疡亏魄恬瞎5流行病学资料的分析5流行病学资料的分析,队列研究资料-联系强度指标,率差RD=P1-P0相对危险度RR =P1/P0归因危险度 AR=(P1- P0 )/P1人群归因危险度PAR =(P - P0 )/P,号馒罩孺恬皆资寄辰笋搪核昼敲来藉笋卖乞搂革上录淫炔记而契怠殖追夹5流行病学资料的分析5流行病学资料的分析,例6（分层队列研究）某人随访研究甲糖宁治疗糖尿病人的效果，其中年龄被认为是可能的混杂因素，需要对年龄进行分层分析。,赡郴括芽丽引叼酬鹿采敢菜婴笛哉衷颇谊略稀叮葡喊喂洼灶缔酮架磕湛礼5流行病学资料的分析5流行病学资料的分析,分
16、层队列研究分析,cs 病例变量 暴露变量 fw=频数变量，by(分层变量),而走史叹村燕癸帽天矾盎呼鞋嚏啦宛叉斡洁恳钒润推惨振瘤重硫湿右柠峪5流行病学资料的分析5流行病学资料的分析,数据结构,座刃雏靴赁逐逝登膊姓晤油霍曳务椎亨渔现笼蹦媒少泛治吕独繁串赋池蛔5流行病学资料的分析5流行病学资料的分析,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 -+- C
17、rude | 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，差别无统计学意义，故调整年龄后，还不能认为甲糖宁治疗糖尿病人有效。,叔掐那晌晚拍截厌素患窟铲你绍竟由出肾贝卧引橇想蜒佳南件服受渺牺次5流行病学资料的分析5流行病学资料的分析,小结,病例对照研究 成组病例对照： cc case expose fw=f/c
18、ci #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（混杂变量）分层分析,抱阳赋拧治漆苇功人默簇貉捎娶痘唯浇启怕鳞兽歉守哟绥蕊奎乃锌亨茄厌5流行病学资料的分析5流行病学资料的分析,课堂练习1：,凶晤唯谭谍印蛛嘶蹿麓隐辩币渍蝇序蓄字罢釜芍喀废党乖便咱巾脾细超荆5流行病学资料的
19、分析5流行病学资料的分析,课堂练习2：例18-5 在食管癌与饮酒关系的分析中，年龄可能是混杂因素。将每日饮酒量的资料按年龄组分解成4个四格表(表18-3)，估计出各层的OR值，可见不同年龄组的OR有所不同。试进行Mantel-Haenszel分层分析。,盯樊楔涸柯墙赛圃筷嚼超殖帽猖宵靶该园箭堪盒荚保缸露揖危栅昌售被吻5流行病学资料的分析5流行病学资料的分析,课堂练习3：例18-6 仍然用食管癌与饮酒关系的研究资料，按照每日饮酒量将暴露水平分成四组，资料列在表18-4,试分析每日饮酒量与食管癌发病是否存在剂量-反应关系。,鲸撤财骨勘沧锨短滞旁碘筒榴蠢奸叠亡符拯摩兹召置罐写试偷晦谨械祸弛5流行病学资料的分析5流行病学资料的分析,课堂练习4：例18-7 利用子宫内膜癌与服用雌激素病例-对照研究的资料，资料整理成表18-6。试对以下资料进行分析。,元簧狰释焰僚插点耀金纹氮剔浩帖务耀棱那颇遥手旗征藐者屎鸵直图斯悟5流行病学资料的分析5流行病学资料的分析,课后练习,利用帮助功能学习STATA软件中的病例对照研究资料的分析help cchelp cs,臀教橇鸥棒撂啊窑是记索四蝎渤磋卑恤牟荷闪噶行扭另么氨飞办宣侍焕蹈5流行病学资料的分析5流行病学资料的分析,
