1、流行病学数据统计分析策略流行病学数据统计分析策略 中南大学湘雅公共卫生学院中南大学湘雅公共卫生学院 谭红专谭红专 教授教授 2023/2/26 2 主要内容主要内容 概述概述 专题研究的数据分析专题研究的数据分析 常规收集的流行病学数据的分析常规收集的流行病学数据的分析 2023/2/26 3 第一节第一节 概述概述 流行病学研究与流行病学数据流行病学研究与流行病学数据 流行病学数据分析的目的流行病学数据分析的目的 流行病学数据中的变量分类流行病学数据中的变量分类 流行病学数据分析的内容流行病学数据分析的内容 2023/2/26 4 一、流行病学研究与流行病学数据一、流行病学研究与流行病学数据
2、 狭义的流行病学数据流行病学数据（epidemiological data）指任何来自针对明确研究目的而开展特定流行病学研究所收集的资料。广义的流行病学数据还包括出于其他目的和用途而收集的、可以用来定量地探索疾病健康或医疗卫生服务其他问题的资料。即一切可为流行病学所用的数据。目前目前的临床大数据是个宝库。的临床大数据是个宝库。2023/2/26 5 二、流行病学数据分析的目的二、流行病学数据分析的目的 流行病学数据分析是针对明确研究问题、具有针对明确研究问题、具有明确目的明确目的、采用特定方法对流行病学数据进行统计整理、统计描述、统计推断和总结。流行病学研究的问题涉及疾病分布、病因与危险因素、
3、诊断、预防和治疗效果评价等。数据分析的根本是正确地回答有意义的问题，而不是统计学数字游戏。2023/2/26 6 流行病学数据分析的目的流行病学数据分析的目的 估计有关统计学指标，如相对危险度；估计该统计学指标的可信区间；控制可能的混杂因素；分析剂量反应关系；分析可能的效应修饰因子；分析可能存在的偏倚。2023/2/26 7 三、流行病学数据中的变量分类三、流行病学数据中的变量分类 典型的流行病学研究中变量可根据其用途分为五类：暴露变量：因，自变量，可以有多个 结局变量：果，因变量，是数据分析的核心 混杂因素：需要控制 效应修饰因素：需要描述 其他变量：2023/2/26 8 按统计学分类按统
4、计学分类 定性：二分类和多分类（无序多分类）定量：在多数流行病学研究中被转化为分类变量 半定量（分级，有序多分类）同一个研究中，因变量和自变量的关系应该是固定的、不能互换。但在不同的研究中，一个变量在这个研究中可能是果，在另一个研究中也可能是因。例如，高血压可能是遗传的结果，血压也可能是心脑血管事件的原因。2023/2/26 9 四、流行病学数据分析的内容四、流行病学数据分析的内容 描述研究对象的数量变动；变量分类和数据整理；描述和比较组间基线资料；估计结局事件发生频率；估计效应大小及其可信区间；识别和控制混杂；识别和测量效应修饰作用；识别和测量剂量反应关系；其他分析。2023/2/26 10
5、 第二节第二节 流行病学专题研究的数据分析流行病学专题研究的数据分析 一、描述研究对象的数量变动一、描述研究对象的数量变动 数据分析的第一步须对研究对象数数据分析的第一步须对研究对象数量随时间的变动进行描述，这是判断选量随时间的变动进行描述，这是判断选择偏倚是否存在及其大小的主要依据；择偏倚是否存在及其大小的主要依据；对选样、抽样和筛选都应交代原则和方对选样、抽样和筛选都应交代原则和方法；对退出和失访都应记录数量和原因法；对退出和失访都应记录数量和原因。2023/2/26 11 2023/2/26 12 二、流行病学变量的分类与数据转换二、流行病学变量的分类与数据转换 转换目的：习惯、专业考量
6、、模型限制 转换类型：定量转定性或分级，定性转定量 转换常用原则：借鉴既往同类研究的分级方法；按照通用的分级标准，如年龄可按国际通用的婴儿、幼儿、青少年、青年、中年、老年的年龄切点来分级；把研究人群分成每组人数均等的35组等。2023/2/26 13 三、描述基线资料三、描述基线资料 基线资料就是有关研究对象代表性和混杂因素的信息资料。基线资料描述的基本目的：一是交待研究对象的特征，提供其代表人群的信息，据此判断结果的外推人群和外推性；二是评估暴露组和非暴露组之间混杂因素的可比性，检查存在混杂的可能性。描述基线资料时，通常在总体描述基础上，对暴露组和非暴露组还要分别进行描述。2023/2/26
7、 14 四、结局事件发生的频率四、结局事件发生的频率 一般要求用率表示，如发病率。发病率由三个因素决定：观察时间、可能发病的总人数和实际发病人数。一般情况下，发病率的分子是在一定时间内发病的总人数。分母则因研究设计不同而异。当发病率较低，且每个人的观察时间长短基本一致，可用研究开始时该组人数作分母计算发病率，叫做累积发病率（cumulative incidence）。当结局事件发生率较高（如晚期癌症治疗试验中的死亡）或每个人的观察时间相差较大时，这时可用人时数（如人年数，person-years）做分母计算发病率，称作人时发病率（person-time incidence rate或发病密度（
8、incidence density），2023/2/26 15 五、估计效应值五、估计效应值 效应就是暴露或治疗对结局作用或影响的大小，多用暴露组与非暴露组或治疗组与对照组间发病率的差别来表达。流行病学研究中常见的效应测量指标有：相对危险度RR=I1/I0；归因危险度百分数ARP=(I1-I0)/I1；相对危险减少度RRR=(I1-I0)/I0；比值比OR=ad/bc；率差RD=I1-I0；需治疗人数NNT=1/RD。2023/2/26 16 例如，某项吸烟和肺癌关系的前瞻性研究，纳入无肺癌的吸烟者 10,000 人，无肺癌的非吸烟者 20,000 人，平均随访 10 年，不同研究对象观察时间
9、的差别不超过 3 个月。10 年内吸烟组中新发肺癌 50 例，非吸烟组新发肺癌 20 例。吸烟对肺癌发病的效应测量指标计算如下：吸烟组 10 年累积发病率 I50/10000=0.005 非吸烟组 10 年累积发病率 I020/20000=0.001 相对危险度 RR=I/I0=0.005/0.001=5.0 率差 RD=I-I0=0.005-0.001=0.004 归因危险度百分数 ARP=(I-I0)/I1=(0.005-0.001)/0.005=0.004/0.005=80%比值比 OR=I(1-I)/I(1-I)=0.005(1-0.001)/0.001(1-0.005)=5.02 需
10、治疗人数 NNT=1/0.004=250 2023/2/26 17 六、估计效应可信区间六、估计效应可信区间 由于随机误差，各效应指标的点估计不能代表效应的真实值，可信区间（Confidence Interval，CI）可用来表达由随机误差引起的效应估计值的不确定性。实践中一般采用95%可信区间，可以将其理解为真实效应值有95%的可能在这个区间之内。若效应指标的点估计值为X，标准误为SE，该效应指标的95可信区间：95CI X1.96SE，其中，（X1.96SE）为其下限，（X1.96SE）为其上限。2023/2/26 18 七、识别和控制混杂七、识别和控制混杂 混杂是由于同一个研究里暴露对某
11、疾病的作用与其他病因的作用交织在一起导致暴露效应估计上的误差。混杂必须同时具备三个条件：其一，是疾病确定的病因或危险因素；其二，不是暴露和疾病间的中间因素；其三，可疑的混杂因素在暴露组和非暴露组的分布不均衡。如性别可能就是吸烟与肺癌关系中的一个混杂因素。2023/2/26 19 混杂的控制就是控制混杂的第三个条件，办法有：设计阶段有三种方法：限制、匹配和随机分组；数据分析阶段也有三种方法：直接标化法、分层分析和多元回归分析。2023/2/26 20 八、交互作用八、交互作用 交互作用不同于混杂。混杂是粗效应值估计时的一种偏倚，需要控制和消除；交互作用是效应大小随第三因素暴露强度或剂量变化而变化
12、的现象。区别交互作用和混杂的方法是一致性检验。如果层间差异无显著性（P0.05），说明不存在交互作用，这时可按混杂处理，可以用一个加权平均效应值来表达，如果平均效应值与粗效应值一致，说明无混杂。2023/2/26 21 交互作用模型及其相对性交互作用模型及其相对性 决策应基于相加模型决策应基于相加模型 表 20-3 肺癌发病率交互作用与效应指标的关系（1/10 万）石棉接触史 吸烟史 有 无 吸烟者 50 10 非吸烟者 5 1 率比 10=50/5 10=10/1 率差 45=50-5 9=10-1 2023/2/26 22 九、识别和测量剂量反应关系九、识别和测量剂量反应关系 识别和测量剂
13、量反应关系的分析方法主要有分层分析和回归分析。表 20-5 吸烟量与肺癌的剂量反应关系 每日吸烟支数 0 1-10-20-30-40-合计 RR 1.0 5.0 8.5 10.2 13.5 15.7 9.8 2023/2/26 23 第三节第三节 常规收集的常规收集的 流行病学数据的分析流行病学数据的分析 常规数据分析的特殊性有三点：常规数据分析的特殊性有三点：第一，对研究问题的设定和对分析结果真实性的评估；第二，由于没有预先设定研究问题，可研究的问题取决于资料内容和性质；第三，设计框架不清楚或不完善，存在多种偏倚，结果真实性较差或不确定。本节将以医院的病人资料为例，说明常规数据分析的目的、方
14、法、步骤和注意事项。2023/2/26 24 一、常规数据的特征一、常规数据的特征（一）数据的时间框架和分类（一）数据的时间框架和分类 时间框架指数据的时间特征，主要分为随访数据（队列数据）和非随访数据（现况数据）。现况数据主要探索病人现况、服务状况和诊断准确性，偶尔用于初步探索病因、疗效和不良反应。队列数据可用于大多队列数据可用于大多数临床问题的探索数临床问题的探索。2023/2/26 25（二）变量的特征和分类（二）变量的特征和分类 按照临床特征可将常规数据中的变量分为治疗措施、诊断结果、病因/危险因素、预后因素、结局、服务相关因素（如药价、医疗保险等）。变量决定可能的研究问题变量决定可能
15、的研究问题。例如，数据中无治疗信息，则无法研究疗效和副作用；无诊断结果，则无从研究诊断准确性。无结局信息，则不能研究疗效、预后因素、危险因素。2023/2/26 26 二、形成研究问题二、形成研究问题 根据数据的变量特征，可以初步形成可研究问题的范围。欲进行诊断准确性研究，数据库里必须有检查结果和疾病诊断信息。研究副作用，基线资料中必须有治疗信息，随访资料中须有副作用的信息。评估疗效，基线资料中须有治疗的信息，随访资料中须有结局信息。研究危险因素，基线资料中须有可疑病因的信息，随访资料中须有有关疾病的信息。对研究者最大的挑战不是如何进行分析，而是如何形成研究问题？2023/2/26 27 三、
16、常规数据三、常规数据“研究设计研究设计”的缺陷的缺陷 表 20-7 常见临床问题的最佳研究设计 研究目的 流行病学最佳设计 治疗的效果 随机对照试验 诊断的准确性 现况研究 预后和预后因素 前瞻性研究 病因和危险因素 前瞻性研究 罕见的病因 病例对照研究 罕见的药物副作用 病例对照研究 病人和服务现状 现况研究 总结现有的研究证据 系统综述 2023/2/26 28 研究设计最重要的三个因素是时间、人群和变研究设计最重要的三个因素是时间、人群和变量。量。时间因素包括时点还是时间段，如何定义起始时间、时间走向、变量间的时间关系；人群因素包括人群特征，选择过程，及比较组的形成；变量因素包括测量指标、何时测量、测量的准确性。通过比较常规数据收集的“设计框架”和最佳研究设计，可以发现常规数据的设计缺陷，判断可能引起的偏倚。在时间、人群和变量特征在时间、人群和变量特征上与最佳设计的差异，就是常规数据的缺陷所上与最佳设计的差异，就是常规数据的缺陷所在，在，也是偏倚可能出现的地方。主要缺陷是非随机主要缺陷是非随机,失访或变量缺失。失访或变量缺失。2023/2/26 29 四、估计相应的指标四、估计相
