1、流行病学的新方法在药物流行流行病学的新方法在药物流行病学的应用病学的应用施侣元施侣元 吕美霞吕美霞华中科技大学同济医学院华中科技大学同济医学院流行病与卫生统计学系流行病与卫生统计学系 自自20世纪世纪60年代以来,被誉为年代以来,被誉为“公共卫生公共卫生之母之母”的流行病学,无论在内容或研究的流行病学,无论在内容或研究方方法法上均取得令人瞩目的进展上均取得令人瞩目的进展 随着医学模式转变,由传统的单纯生物型向着生物随着医学模式转变,由传统的单纯生物型向着生物-心心理理-社会的新型模式转变。社会的新型模式转变。流行病学研究内容从所有与人类健康相关的卫生事件流行病学研究内容从所有与人类健康相关的卫
2、生事件(health events)的各个方面,甚至超出卫生事件的的各个方面,甚至超出卫生事件的自然自然 和社会问题的挑战,无一不是新世纪面对的新课和社会问题的挑战,无一不是新世纪面对的新课题,如,全球气候变暖、厄尔尼诺现象,超大型城市题,如,全球气候变暖、厄尔尼诺现象,超大型城市不断增加,社会动乱及犯罪等等,无一不能用流行病不断增加,社会动乱及犯罪等等,无一不能用流行病学方法去研究。学方法去研究。随着学科彼此间的交融,渗透使流行病学在研究方法随着学科彼此间的交融,渗透使流行病学在研究方法上出现了突破性的进展。上出现了突破性的进展。古语说:古语说:“工欲善其事,必先利其器工欲善其事,必先利其器
3、”一、常规流行病学方法一、常规流行病学方法 研究方法研究方法观观察察法法实实验验法法数数理理法法描述流行病学描述流行病学 分析流行病学分析流行病学 横断面调查横断面调查监测监测生态学研究生态学研究病例对照研究病例对照研究队列研究队列研究实验流行病学实验流行病学 临床试验临床试验现场试验现场试验社区干预项目社区干预项目理论流行病学理论流行病学 产生假设产生假设 检验假设检验假设 验证假设验证假设 二、多阶段复合设计研究二、多阶段复合设计研究 20世纪世纪70年代开始,提出混合型研究设计。年代开始,提出混合型研究设计。巢式病例对照研究(巢式病例对照研究(nested case-control st
4、udy)病例病例-队列研究(队列研究(case-cohort study)病例病例-病例研究(病例研究(case-only study)随访患病率研究(随访患病率研究(follow-up prevalence study)两阶段病例对照研究(两阶段病例对照研究(two-phase case-control study)病例交叉设计(病例交叉设计(case-crossover study)等等。等等。设计效率高、花费少、适用范围广,便于计算机处理和分析数据。设计效率高、花费少、适用范围广,便于计算机处理和分析数据。特点:特点:(一)巢式病例对照研究(一)巢式病例对照研究 nested case-
5、control study 1973年美国流行病学家年美国流行病学家Mantel提出的提出的综合式综合式病病例对照研究设计,例对照研究设计,1982年正式提出年正式提出巢式病例对巢式病例对照研究照研究。又称套叠式或嵌入式病例对照研究,也称队列又称套叠式或嵌入式病例对照研究,也称队列内病例对照研究(内病例对照研究(case-control study nested in a cohort),),巢式病例对照研究巢式病例对照研究概念:概念:是以队列中所有的病例作为病例组,再是以队列中所有的病例作为病例组,再根据病例发病时间,在研究队列的非病例根据病例发病时间,在研究队列的非病例中随机中随机匹配匹配
6、(按年龄、性别、住址、民族(按年龄、性别、住址、民族等)一个或多个对照,组成对照组。可采等)一个或多个对照,组成对照组。可采用传统的用传统的个体匹配病例对照研究的分析方个体匹配病例对照研究的分析方法法或比例危险回归模型分析方法。或比例危险回归模型分析方法。分类:分类:前瞻性巢式病例对照研究(前瞻性巢式病例对照研究(prospective nested case-control study)回顾性巢式病例对照研究(回顾性巢式病例对照研究(retrospective nested case-control study)特点特点:在疾病诊断前收集暴露资料,选择偏倚和信息偏在疾病诊断前收集暴露资料,选
7、择偏倚和信息偏倚较小;倚较小;研究中的病例与对照来自于同一队列,其可比性研究中的病例与对照来自于同一队列,其可比性好;好;统计效率和检验效率高于传统病例对照研究,可统计效率和检验效率高于传统病例对照研究,可以计算发病率以计算发病率 研究样本量小于经典的队列研究,节约人力和物研究样本量小于经典的队列研究,节约人力和物力;力;符合因果推断的时间顺序,论证强度高符合因果推断的时间顺序,论证强度高 实施方法实施方法确定某一人群作为研究队列确定某一人群作为研究队列收集队列内每个成员的收集队列内每个成员的相关信息和生物标本相关信息和生物标本随访一段预定的时间随访一段预定的时间确定随访期内发生研究确定随访期
8、内发生研究病例的全部病例组成病例组病例的全部病例组成病例组用用危险集危险集抽样为每个病例抽取抽样为每个病例抽取一定数量的对照组组成对照组一定数量的对照组组成对照组抽取已收集到的两组成员相关信息抽取已收集到的两组成员相关信息和生物标本做必要的化验作统计分析,和生物标本做必要的化验作统计分析,计算发病率、相对危险度,进行分层分析,计算发病率、相对危险度,进行分层分析,条件条件Logistic回归分析或回归分析或Cox回归分析回归分析得出结论得出结论实例实例 一项环氧化酶选择性和非选择性抑制非甾体类一项环氧化酶选择性和非选择性抑制非甾体类抗炎药与严重性冠心病关系的巢式病例对照研抗炎药与严重性冠心病关
9、系的巢式病例对照研究,其过程和结果如下:究,其过程和结果如下:(Lancet 2005;365:47581)1、研究目的、研究目的 比较罗非昔布和传统非甾体类抗炎药或塞来昔布发比较罗非昔布和传统非甾体类抗炎药或塞来昔布发生严重冠心病的危险性生严重冠心病的危险性2、研究队列、研究队列 美国加利福尼亚州美国加利福尼亚州Kaiser Permanente医疗保健机构医疗保健机构1999年年1月月1日到日到2001年年12月月31日期间使用日期间使用环氧化酶环氧化酶2选择性抑选择性抑制剂(塞来昔布或罗非昔布)或非选择性(其他)非甾体制剂(塞来昔布或罗非昔布)或非选择性(其他)非甾体类抗炎药物治疗一次以
10、上的所有类抗炎药物治疗一次以上的所有18-84岁的人为观察队列岁的人为观察队列。非肿瘤、肾衰、肝衰、严重呼吸疾病、器官移植或。非肿瘤、肾衰、肝衰、严重呼吸疾病、器官移植或HIV/AIDS。队列成员共队列成员共39639人(男人(男24430,女,女15209)。)。3、资料收集、资料收集 收集所有观察对象人口学、治疗、死亡状态、死因等资料收集所有观察对象人口学、治疗、死亡状态、死因等资料4、随访发现严重心脏疾病患者、随访发现严重心脏疾病患者 研究队列中研究队列中1999年至年至2001年期间共随访年期间共随访2302029人年,出人年,出现现8199个病例。个病例。5、用匹配的方法选择对照、用
11、匹配的方法选择对照 每个病例从队列中选择每个病例从队列中选择4个对照,按照年龄、性别和个对照,按照年龄、性别和卫生计划区域(南或北)配对,排除卫生计划区域(南或北)配对,排除了了56例病例和例病例和1300例对照,最后符合标准的有例对照,最后符合标准的有8143例病例和例病例和31496个对照。个对照。v联系强度的计算联系强度的计算 比值比(比值比(odds ratio OR)OR=1 无关联无关联 95%CL 比值比与关联强度的关系比值比与关联强度的关系 比值比关联的强度比值比关联的强度 0.91.0 1.01.1 无无 0.70.8 1.21.4 弱弱 0.40.6 1.52.9 中等中等
12、 0.10.3 3.09.9 强强 10 很强很强 6、资料分析、资料分析 应用条件应用条件logistic回归进行分析结果:回归进行分析结果:罗非昔布罗非昔布(各剂量各剂量)对塞来昔布调整对塞来昔布调整OR为为1.59(95%CI:1.10-2.32,P=0.015)25mg/day以下以下罗非昔布的罗非昔布的OR=1.47(95%CI:0.99-2.17,P=0.054)25mg/day以上以上OR=3.58(95%CI:1.27-10.11,P=0.016)萘普生对既往萘普生对既往NSAID的的OR为为1.14(95%CI:1.00-1.30,P=0.05)7、结论、结论 结果表明罗非昔
13、布较塞来昔布增加严重性冠心结果表明罗非昔布较塞来昔布增加严重性冠心病的危险,病的危险,而萘普生没有而萘普生没有抗严重性冠心病的作抗严重性冠心病的作用。用。表表1 服用选择性非甾体类抗炎药和传统非甾体类抗炎药服用选择性非甾体类抗炎药和传统非甾体类抗炎药或塞来昔布发生急性心肌梗塞危险性或塞来昔布发生急性心肌梗塞危险性病例病例对照对照粗粗OR(95%CI)调整调整OR(95%CI)P与既往使用与既往使用NSAIDs比较比较既往使用既往使用4658187201.001.00近来使用近来使用172062581.12(1.05-1.20)1.11(1.03-1.19)0.004当前使用当前使用塞来昔布塞来
14、昔布126 4911.05(0.86-1.28)0.84(0.67-1.04)0.12布洛芬布洛芬67025731.07(0.98-1.18)1.06(0.96-1.17)0.27萘普生萘普生36714091.07(0.95-1.21)1.14(1.00-1.30)0.05罗非昔布(所有剂量)罗非昔布(所有剂量)681961.39(1.05-1.83)1.34(0.98-1.82)0.066罗非昔布罗非昔布25mg/day 581881.23(0.92-1.66)1.23(0.89-1.7100.21罗非昔布罗非昔布25mg/day 10 85.03(1.98-12.76)3.00(1.09-
15、8.31)0.03其他其他NSAIDs53418491.19(1.07-1.32)1.13(1.01-1.27)0.03调整年龄、性别、卫生计划区,心血管危险性得分,非心脏不良的治疗,给药时间,非心脏原因的急诊观察,激素替代治疗,高剂量的强的松治疗。调整年龄、性别、卫生计划区,心血管危险性得分,非心脏不良的治疗,给药时间,非心脏原因的急诊观察,激素替代治疗,高剂量的强的松治疗。病例病例对照对照粗粗OR(95%CI)调整调整OR(95%CI)P与塞来昔布比较与塞来昔布比较塞来昔布塞来昔布1264911.001.00既往使用既往使用4658187200.95(0.78-1.16)1.11(0.96
16、-1.48)0.12近来使用近来使用172062581.07(0.87-1.31)1.32(1.06-1.65)0.015当前使用当前使用布洛芬布洛芬67025731.02(0.82-1.27)1.26(1.00-1.60)0.054萘普生萘普生36714091.02(0.81-1.28)1.36(1.06-1.75)0.016罗非昔布(所有剂量)罗非昔布(所有剂量)681961.32(0.94-1.85)1.59(1.10-2.32)0.015罗非昔布罗非昔布25mg/day581881.17(0.82-1.67)1.47(0.99-2.17)0.054罗非昔布罗非昔布25mg/day10
17、84.78(1.85-12.38)3.58(1.27-10.11)0.016其他其他NSAIDs53418491.13(0.91-1.41)1.35(1.06-1.72)0.015续表续表1 服用选择性非甾体类抗炎药和传统非甾体类抗炎药服用选择性非甾体类抗炎药和传统非甾体类抗炎药或塞来昔布发生急性心肌梗塞危险性或塞来昔布发生急性心肌梗塞危险性调整年龄、性别、卫生计划区,心血管危险性得分,非心脏不良的治疗,给药时间,非心脏原因的急诊调整年龄、性别、卫生计划区,心血管危险性得分,非心脏不良的治疗,给药时间,非心脏原因的急诊观察,激素替代治疗,高剂量的强的松治疗。观察,激素替代治疗,高剂量的强的松治
18、疗。(二)病例(二)病例-队列研究(队列研究(case-cohort study)也称病例参比式研究是也称病例参比式研究是1986年年Prentice RLPrentice RL提出提出 概念:概念:是以队列中所有随访的病例作为病例组,从整个是以队列中所有随访的病例作为病例组,从整个队列中采用队列中采用随机或分层随机抽样随机或分层随机抽样的办法,选取一定比例的办法,选取一定比例的样本作为对照组,然后用一定的统计方法比较分析两的样本作为对照组,然后用一定的统计方法比较分析两组资料,以探索影响疾病发生、疾病生存时间、预后等组资料,以探索影响疾病发生、疾病生存时间、预后等的因素。的因素。特点特点 是
19、将队列设计和病例对照设计相互交叉,融合两者是将队列设计和病例对照设计相互交叉,融合两者的优点后而形成的一种设计方法。的优点后而形成的一种设计方法。资料统计方法相对危险度的资料统计方法相对危险度的准准似然函数估计法(分似然函数估计法(分析多变量)、改进的相对危险度估计法(分析单变析多变量)、改进的相对危险度估计法(分析单变量)。专用统计软件(量)。专用统计软件(EPICURE)处理。处理。实施方法实施方法选择队列选择队列确定病例组确定病例组确定对照组确定对照组收集资料收集资料统计分析统计分析该病发病率发病人数子队列的人数1用病例对照研究的分析方法估计用病例对照研究的分析方法估计OR值,值,用准似
20、然危险度估计用准似然危险度估计RR值,值,用用Cox模型进行多因素分析等等。模型进行多因素分析等等。估计关联强度估计关联强度 准似然危险度(准似然危险度(pseudo risk):):将子队列按有暴露史和无暴露史将子队列按有暴露史和无暴露史分为两个相对人群,分别计算暴露组和非暴露组的发病比率。分为两个相对人群,分别计算暴露组和非暴露组的发病比率。111BAriskpseudo000BAriskpseudo0011BABARR A1和和A0分别表示子队列暴露组中的发病人数和非暴露组的发病分别表示子队列暴露组中的发病人数和非暴露组的发病人数,人数,B1和和B0有暴露组的人数和无暴露组的人数。有暴露
21、组的人数和无暴露组的人数。病例病例-队列研究与巢式病例对照研究不同奌:队列研究与巢式病例对照研究不同奌:1、对照是随机选取,不与病例进行配比。、对照是随机选取,不与病例进行配比。2、随机对照中的成员如发生被研究的疾病,既作、随机对照中的成员如发生被研究的疾病,既作为对照,又同时作为病例。为对照,又同时作为病例。3、可以同时研究几种疾病,不同疾病有不同的病、可以同时研究几种疾病,不同疾病有不同的病例组,但对照组都是同一组随机样本。例组,但对照组都是同一组随机样本。实例实例 一项关于血液透析病人注射右旋糖酐铁出现副一项关于血液透析病人注射右旋糖酐铁出现副作用的病例作用的病例-队列研究,其过程和结果
22、如下:队列研究,其过程和结果如下:(Am J Kidney Dis 2001;37:743-749,859-861)1、研究目的、研究目的 探讨血液透析病人注射右旋糖酐铁发生副作用探讨血液透析病人注射右旋糖酐铁发生副作用的危险性的危险性 2、确定队列、确定队列 以北美以北美Fresenlus医疗保健中心全体血液透析病医疗保健中心全体血液透析病人作为观察队列。全队列建于人作为观察队列。全队列建于1998年年10月,随月,随访至访至1999年年3月,共月,共6个月,队列成员中共有个月,队列成员中共有86479人注射过右旋糖酐铁。人注射过右旋糖酐铁。3、资料收集、资料收集 人口学特征、药物暴露史、副
23、作用等资料人口学特征、药物暴露史、副作用等资料 4、随访副作用、随访副作用 在在1998年年10月至月至1999年年3月间月间841252次静脉注次静脉注射右旋糖酐铁中共发生副作用射右旋糖酐铁中共发生副作用165次,发生率次,发生率为为0.000196%,约,约20/10万万 5、资料分析、资料分析 用非参数统计比较该病例用非参数统计比较该病例-队列研究中的两组比队列研究中的两组比较,较,Dexferrum和和InFed之间的比较之间的比较t检验和检验和Fish确切概率分析。估计确切概率分析。估计OR值及值及95%CI。应用应用Dexferrum的患者较接受的患者较接受InFed治疗的患者药物
24、治疗的患者药物副作用发生率高副作用发生率高8.1倍。倍。6、随访发生的副作用、随访发生的副作用 临床上由右旋糖酐铁所至的严重副作用罕见,临床上由右旋糖酐铁所至的严重副作用罕见,难以预测,和使用的剂型有关。难以预测,和使用的剂型有关。表表2 与右旋糖酐铁相关的主要的药物不良反应发生率与右旋糖酐铁相关的主要的药物不良反应发生率症状或体征症状或体征病例比例(病例比例(%)呼吸困难呼吸困难43低血压低血压23胸痛胸痛19意识丧失意识丧失*12血管性水肿血管性水肿9可逆性精神症状可逆性精神症状*6癫痫癫痫*5心律失常心律失常3呼吸暂停呼吸暂停2心跳停止心跳停止1死亡死亡0.6*包括神经症状包括神经症状表
25、表3 与右旋糖酐铁相关的轻微不良反应与右旋糖酐铁相关的轻微不良反应症状或体征症状或体征病例比例(病例比例(%)恶心恶心34面色潮红面色潮红27搔痒症搔痒症25呕吐呕吐23皮疹皮疹16腹部痛腹部痛13发汗发汗10颤动颤动6腹泻腹泻5不安不安3肌痛肌痛3疲劳疲劳2其他(视力模糊、抽筋等)其他(视力模糊、抽筋等)6 Dexferrum:P1=110/166228=6.62/10万万 InFed:P0=55/675024=0.815/10万万 OR=(110/166228)/(55/675024)=8.12 OR 95%CI为为8.08-8.17评论评论 奥地利维也纳大学医院的奥地利维也纳大学医院的H
26、orl博士在相关评论中与博士在相关评论中与Chertow博士等博士等讨论说:讨论说:“尽管作者们证实,在临床上由右旋糖酐铁所至的严重尽管作者们证实,在临床上由右旋糖酐铁所至的严重副作用罕见,但用蔗糖铁和葡萄糖铁替代右旋糖酐铁治疗更安全副作用罕见,但用蔗糖铁和葡萄糖铁替代右旋糖酐铁治疗更安全。”Horl博士强调:博士强调:“不能因为价格较贵而对肾性贫血的终末期肾病不能因为价格较贵而对肾性贫血的终末期肾病患者不使用这些替代性铁剂,患者不使用这些替代性铁剂,”他总结说:他总结说:“应该应用应该应用rhuEPO 和更安全的铁剂纠正患者的贫血,而不是右旋糖酐铁。和更安全的铁剂纠正患者的贫血,而不是右旋糖
27、酐铁。”Chertow博士指出:博士指出:“静脉注射无右旋糖酐的铁剂替代右旋糖酐铁静脉注射无右旋糖酐的铁剂替代右旋糖酐铁治疗是一个复杂的问题,需要从危险性、潜在效益和相对价格等治疗是一个复杂的问题,需要从危险性、潜在效益和相对价格等方面进行综合考虑。方面进行综合考虑。”(三)病例三)病例-病例研究病例研究 (case-only study)探讨不同临床类型或具有某方面标志的病例与探讨不同临床类型或具有某方面标志的病例与无标志的病例与危险因素之间的关系和相互作无标志的病例与危险因素之间的关系和相互作用,按病例对照研究的方式处理资料,以探讨用,按病例对照研究的方式处理资料,以探讨不同临床类型的危险
28、因素的差异或者这个生物不同临床类型的危险因素的差异或者这个生物学标志物与该疾病之间的关系和相互作用。学标志物与该疾病之间的关系和相互作用。由于研究对象均为患病的病例,不另外设未患由于研究对象均为患病的病例,不另外设未患病的对照组,直接采取对患病的两个亚组比较病的对照组,直接采取对患病的两个亚组比较,并按病例对照研究的方式处理资料。,并按病例对照研究的方式处理资料。特点特点 可免除选择对照,特别是免除从对照组中可免除选择对照,特别是免除从对照组中去收集资料,尤其是收集生物标本的麻烦去收集资料,尤其是收集生物标本的麻烦。实例实例 陈建锋等人在江苏太仓、如东两地收集脑卒中陈建锋等人在江苏太仓、如东两
29、地收集脑卒中的妇女,采取病例的妇女,采取病例-病例研究方法探讨了复方病例研究方法探讨了复方口服避孕药(口服避孕药(COC)在脑卒中所起的作用在脑卒中所起的作用,旨在进一步澄清复方口服避孕药(旨在进一步澄清复方口服避孕药(COC)对对脑卒中可能的危险。脑卒中可能的危险。发现110例女性脑卒中患者中,用病例-病例研究方法,分析复方口服避孕药(口服避孕药(COC)在脑在脑卒中所起的作用中所起的作用,在影响出血性和缺血性脑卒中之间的差异。表表4 单因素分析结果单因素分析结果 因素因素OROR95%CI2PCOC3.031.27-7.297.670.01BMI242.991.17-7.756.470.0
30、1与高血压、高血脂、糖尿病、生活紧张程度、口味偏咸等无关(与高血压、高血脂、糖尿病、生活紧张程度、口味偏咸等无关(P0.05)。Logistic 回归模型显示回归模型显示:最后进入模型的变量是服用最后进入模型的变量是服用COC和和 BMI24,服用服用COC与出血性脑卒中的联系几乎与出血性脑卒中的联系几乎是缺血性脑卒中的是缺血性脑卒中的3倍(倍(OR=2.93,95%CI:1.28-6.72,2=6.48,P=0.01)因素因素脑卒中例数脑卒中例数OROR95%CI趋势检验趋势检验出血性出血性缺血性缺血性2P服用时间(月)服用时间(月)不服用不服用28331.001871.350.38-4.8
31、160952.120.56-8.3711.560.00118018210.612.07-72.65与从未使用与从未使用COC相比,相比,COC发生出血性脑卒中的危险性较发生发生出血性脑卒中的危险性较发生缺血性脑卒中大,随着使用年限的增加有增大的趋势,特别是使缺血性脑卒中大,随着使用年限的增加有增大的趋势,特别是使用用15年以上时,这种趋势就更为明显(年以上时,这种趋势就更为明显(OR=10.61)。)。表表5 使用使用COC时间与妇女脑卒中的关系时间与妇女脑卒中的关系 因素因素脑卒中例数脑卒中例数OROR95%CI趋势检验趋势检验出血性出血性缺血性缺血性2P停药时间(月)停药时间(月)不服用不服用28331.0011126.481.19-46.47601162.160.63-7.644.070.0431801362.550.77-8.78表表6 停用停用COC时间与妇女脑卒中的关系时间与妇女脑卒中的关系 随着随着 停用停用用年限的增加有下降的趋势,但停用用年限的增加有下降的趋势,但停用15年以上时,年以上时,这种不良反应还存在(这种不良反应还存在(OR=2以上)。以上)。Thank you all!