Download
1 / 72
870 likes | 1.46k Views
交互作用与效应修饰. 交互作用 ( interaction). 指两个或多个因素相互依赖发生作用而产生的一种效应。生物学上是指两个或多个因素互相依赖发生作用而引起疾病或预防和控制疾病。即生物学交互作用 (biological interaction). 研究资料中各研究因素间存在交互作用时，说明各研究因素的作用不是独立的，一个因素水平改变时与它有交互作用的因素之效应也将发生改变 。 交互作用是描述资料本身所表现出来的关于两个或多个暴露因素联合效应的一个统计学术语，并不表示一定具有生物学意义。有交互作用可能有生物学意义，也可能没有生物学意义。. 效应修饰作用.
E N D
交互作用( interaction) • 指两个或多个因素相互依赖发生作用而产生的一种效应。生物学上是指两个或多个因素互相依赖发生作用而引起疾病或预防和控制疾病。即生物学交互作用(biological interaction) .
研究资料中各研究因素间存在交互作用时，说明各研究因素的作用不是独立的，一个因素水平改变时与它有交互作用的因素之效应也将发生改变 。 • 交互作用是描述资料本身所表现出来的关于两个或多个暴露因素联合效应的一个统计学术语，并不表示一定具有生物学意义。有交互作用可能有生物学意义，也可能没有生物学意义。
效应修饰作用 • 当暴露因素按第三变量分层后估计暴露在每一层中与疾病的联系强度时，效应修正被定义为暴露因素在各层中与疾病的联系强度(测量的效应)因第三变量的存在情况不同而大小不同。 • 该第三变量称为效应修正因素( effect modification factor）。
第三变量在一项研究中是否成为效应修正因素，取决于选用判断暴露和疾病之间联系的指标是用率差( rate difference , RD)还是用率比( rate ratio , RR) 。 • 又称效应变异（effect variation），或效应不一致性或异质性（heterogeneity of effect）。
例如，性别与髋骨骨折的联系受到年龄的修饰；例如，性别与髋骨骨折的联系受到年龄的修饰； • 女性髋骨骨折危险性为男性的2至3倍 • 年轻组：男性髋骨骨折危险性高于女性 • 年老组：女性髋骨骨折危险性明显高于男性 • 说明性别与骨折的联系被年龄修饰。
体重与乳癌的联系受到绝经状态的修饰； • 肥胖与乳癌之间的联系被绝经状态修饰， • 肥胖是绝经后妇女乳癌的危险因素， • 肥胖不是绝经前妇女乳癌的危险因素。
孕妇吸烟与婴儿低出生体重的联系 • 吸烟对出生低体重的有害影响在高龄孕妇较低龄孕妇更为明显。 • 母亲年龄是修饰因素。
EB病毒感染与非洲儿童淋巴瘤的联系 • 生命早期感染EB病毒与非洲儿童Burkitt淋巴瘤的联系多见于幼年生活在疟疾广泛流行区的儿童，少见于疟疾非流行区。 • 当地普遍存在疟疾是一种效应修饰因素。
效应修饰与混杂 • 混杂是一种偏倚，是研究者希望避免的，或在必要时希望从资料中消除的。 • 效应修饰是对效应本身的详细描述，它是一种需要报告的发现，而不是一种需避免的偏倚。 • 流行病学分析总目的是消除混杂，并发现和描述效应修饰。
混杂是否存在取决于研究设计。 • 混杂源于选择研究对象的源人群中混杂变量与研究因素的相关性。因此，在选择研究对象时对研究对象进行适当限制，可以防止一个变量成为混杂因素。 • 效应修饰并不取决于研究设计，是一种与研究设计无关的自然现象，是研究欲尽可能揭示和描述的现象，它超出了恒定的研究设计概念。
效应修饰是一种恒定的自然现象，这是一种相对（非绝对）的概念，故其并不与任何生物学特性相对应。效应修饰是一种恒定的自然现象，这是一种相对（非绝对）的概念，故其并不与任何生物学特性相对应。 • 在研究效应修饰时必须规定欲测量的是何种效应。 • 广义上讲，效应修饰是指对一种效应的修饰，但并未指明修饰何种效应测量。流行病学中有两种常用的危险性效应测量，即率差和率比，以及其它一些不常用的测量。 • 如果不规定修饰何种测量（率差或率比），效应修饰的概念就太模糊，使人无法描述其特性。
暴露组与非暴露组之间的发病率比在不同年龄组是恒定的。暴露组与非暴露组之间的发病率比在不同年龄组是恒定的。 • 但其发病率差则随年龄而增加。以率差为指标，年龄是效应修饰因素。
暴露组与非暴露组之间的发病率差在不同年龄时是恒定的。暴露组与非暴露组之间的发病率差在不同年龄时是恒定的。 • 对率比这个效应来说，年龄修饰暴露效应。
例如 • 如果每天喝5杯咖啡，可使男性心肌梗死发生率增加40%，（率比RR=1.4），但使女性心肌梗死发生率只增加10%（RR=1.1），这种情况称作性别对RR的修饰（或RR的性别变异或性别的RR不一致性），性别即为喝咖啡致心肌梗死RR的修饰因素。
又例如，每天喝5杯咖啡，可使男性心肌梗死发生率增加到400/10万人年，但只使女性心肌梗死发生率增加到40/10万人年。这种情况称为性别对率差（AR）的修饰，性别即为喝咖啡致心肌梗死AR的修饰因素。又例如，每天喝5杯咖啡，可使男性心肌梗死发生率增加到400/10万人年，但只使女性心肌梗死发生率增加到40/10万人年。这种情况称为性别对率差（AR）的修饰，性别即为喝咖啡致心肌梗死AR的修饰因素。
再一例子，如果每天喝5杯咖啡，无论对男性或女性，都可使心肌梗死发生率增加22%（RR=1.22），这种情况称为性别间的RR一致性（homogeneity），也即性别不存在效应修饰。再一例子，如果每天喝5杯咖啡，无论对男性或女性，都可使心肌梗死发生率增加22%（RR=1.22），这种情况称为性别间的RR一致性（homogeneity），也即性别不存在效应修饰。 • 需注意的是，效应修饰和效应不一致性并非效应的固有特性，而仅是效应测量方法的一个特征。
例如，每天喝5杯咖啡，可使男性心肌梗死发生率从1000/10万人年增加到1220/10万人年，而使女性心肌梗死发生率从400/10万人年增加到488/10万人年，男女性的RR均=1.22，呈现一致性。例如，每天喝5杯咖啡，可使男性心肌梗死发生率从1000/10万人年增加到1220/10万人年，而使女性心肌梗死发生率从400/10万人年增加到488/10万人年，男女性的RR均=1.22，呈现一致性。 • 男性的AR=220/10万人年，女性AR=88/10万人年，男女性的AR不一致或AR被性别所修饰。 • 这个例子说明，效应修饰不一定等于生物学的交互作用（interaction）（如协同或拮抗作用）。
效应修饰作用的分析 • 效应修饰可采用分层的方法进行分析，估计每层的效应，并进行比较。 • 但如果每层观察对象的数目太少，使每层的效应不稳定，尤其在调整了混杂因素后。这样，由于随机误差，层间的效应估计可能波动很大。
检验分层资料不一致性的统计学方法，在许多情况下检验效率很低，用常规的显著性水平（如α=0.05），多数的层间不一致性无法检出。用数学模型拟合的方法可以解决其中的一些问题。但由于混杂、测量误差引起的偏倚在各层间不同，对效应修饰无论采用观察或数学模型拟合的方法，仍会有一些偏倚。检验分层资料不一致性的统计学方法，在许多情况下检验效率很低，用常规的显著性水平（如α=0.05），多数的层间不一致性无法检出。用数学模型拟合的方法可以解决其中的一些问题。但由于混杂、测量误差引起的偏倚在各层间不同，对效应修饰无论采用观察或数学模型拟合的方法，仍会有一些偏倚。
交互作用模型 • 相加模型(additive model) • 相加模型假定若交互作用不存在时，两个或两个以上因子共同作用于某一事件时，其效应等于这些因子单单独作用时的和，有时称之具有可加性。
相乘模型( multiplicative model) • 相乘模型假定若交互作用不存在时，两个或两个以上因子共同作用于某一事件时，其效应等于这些因子单独作用时的积。
相加模型：R11－R00＝40-1=39 (R10-R00)+(R01-R00)=8-1+5-1=12 相乘模型：R11/R00=40/1=40 (R10/R00)(R01/R00)=(8/1)(5/1)=40
混杂与交互作用的存在形式 • 四种形式：①有混杂但无交互作用②无混杂但有交互作用③混杂和交互作用同时存在④混杂和交互作用均不存在。 • 如一项研究资料的cRR与aRR不同， 表明其中存在着混杂。将资料按某种特性分层(如可疑的效应修正因子) ，计算各层相对危险度，若层别相对危险度不等，说明资料中存在着交互作用。
例如：某前瞻性研究观察某团体中甲苯磺丁脲治疗组和安慰剂组各种原因死亡数情况。例如：某前瞻性研究观察某团体中甲苯磺丁脲治疗组和安慰剂组各种原因死亡数情况。
交互作用的识别 • 交互作用识别的一般过程 • 在判断交互作用存在与否之前，首先要明确所研究的因素与事件之间是否存在统计学联系。如果有联系，接下来看这一联系是否由偏倚或混杂所致。若存在偏倚或混杂因素，采用适当的方法加以改进和调整后再分析交互作用。
叉生分析 • 叉生分析：x 、y 为两个危险因素, ORx 、ORy 、ORxy分别表示仅有x ,仅有y 和x 、y 都存联合作用模型。若:ORxy > ORx + ORy - 1, 则认为x 、y 间可能存在交互作用。但要作统计学推断,则需作χ2 检验。
叉生分析仅可以给人们提示某两个因素间是否可能存在交互作用,并且叉生分析的统计学检验方法本身尚存在一定的缺陷;另外,在用叉生分析法进行交互作用分析时,未考虑其他不参加叉生的危险因素的作用。因此, 可在多因素logistic 回归模型中加入交互作用项,以平衡其他因素的影响。
分层分析 • 分层分析是比较经典的识别交互作用的方法。可以通过比较按照可疑交互因素分层后层间的效应测量值－相对危险度(RR) 或率差(RD) 来判断是否产生交互作用。如果各层之间的效应测量值RR 或RD 不同，则可能存在交互作用。但是鉴于各层RR 和RD变异可能是机遇所致，因此必须进行统计学检验。
分层分析 • 如果效应一致性的假设不合理，那么按第三变量分层后，各层中效应不一致的情况，不应采用简单合并估计效应的方法。 • 如果第三变量各层间的效应不同时，即存在效应修饰，资料分析和表达的关键是描述分层因素（即第三变量）是如何修饰效应的。
分层分析时的注意问题 • 评价效应修饰常依赖于机械地应用统计学检验方法。 • 建立在统计学效应评价基础上的流行病学交互作用，有时是难以捉摸的，用纯粹机械的方法会使问题变得模糊不清。 • 用统计学方法评价效应修饰时应小心。 • 由于检验的备择假设是非特异性的，一般用于效应修饰检验的普通统计学方法，效率较低。因此，“无显著性的”P值较难正确解释。
由于选择偏倚、分类错误、混杂、以及其它偏倚和致病效应，很少期望用任何尺度测量的效应都会一样精确。由于选择偏倚、分类错误、混杂、以及其它偏倚和致病效应，很少期望用任何尺度测量的效应都会一样精确。 • 一致性效应的无效假设常只是一项统计学设计，这类统计学设计应只能被看作是近似真实情况，仍存有疑问。
分层分析的步骤 • 第一步是按第三变量分层后，检查各层的效应估计。即使基本参数一致，层间估计值还会有某些随机变异，但过大的随机变异或明显的非随机变异可能就是效应修饰的证据。 • 关于效应修饰的判断不应只限于资料的一些表面现象，如果有效应修饰的话，第二步还应在评价过程中结合以往研究的其它知识或一般生物学知识予以考虑。
第三步，用多种统计学检验方法来检验层间一致性。第三步，用多种统计学检验方法来检验层间一致性。 • 不同的危险性指标，如率差和率比，其效应修饰的测量要求有不同的评价方法，比的测量一致性意味差的测量效应修饰，反之，差的测量一致性可能意味比的测量效应修饰。 • 统计学检验效应一致性的无效假设，一般有两类，一是按一致效应的直接加权合并估计，另一是按最大似然估计。
效应修饰只能按某一种效应测量指标来进行描述。效应修饰只能按某一种效应测量指标来进行描述。 • 如果，按率差或按率比，都没有效应修饰，除非未暴露组的疾病率与可能的效应修饰因素之间没有联系，其仍可能在其它测量指标（如SMR）上存在效应修饰。 • 处理效应修饰的基本分析方法是分层分析。 • 将资料按第三变量进行分层分析有两个目的，一个是评价和消除混杂，另一个是评价和描述效应修饰。 • 分层分析是处理这两种问题的好方法。
以色列缺血性心脏病研究 • 心肌梗死与收缩血压关系的病例对照研究──────────────────── 心肌梗死 对照 计──────────────────── • 收缩血压≥140 29 711 740 • 收缩血压< 140 27 1244 1271 • 计 56 1955 2011────────────────────cOR=1.88，χ2=5.56
以色列缺血性心脏病研究（分层） • 心肌梗死与收缩血压关系的病例对照研究────────────────────年龄≥60岁 心肌梗死 对照 计──────────────────── • 收缩血压≥140 9 115 124 • 收缩血压< 140 6 73 79 • 计 15 188 203─────────────────OR1 = 0.95, χ12 = 0.008
以色列缺血性心脏病研究（分层） • ────────────────────年龄< 60岁 心肌梗死 对照 计──────────────────── • 收缩血压≥140 20 596 616 • 收缩血压< 140 21 1171 1192 • 计 41 1767 1808────────────────────OR2 = 1.87, χ22 = 4.04
以色列缺血性心脏病研究 • cOR被看作是不同年龄层ORi（本例中为OR1和OR2）的平均值，如果被平均的各个ORi相同，并包括i个个体的所有样本估计，则这种平均值可以解释为只存在抽样变异。 • 如果血压与心肌梗死的联系在不同年龄层是各异的，我们就称年龄与血压之间存在交互作用，Miettinen认为将其称之为效应修饰更合适。
交互作用的检验 • 上例，两个年龄层的ORi分别为0.95和1.87，也许OR1=0.95反映了年龄≥60岁者的真正OR=1.00，OR2=1.87反映了年龄<60岁者的真正OR=2.00。 • 如果各层的ORi确实不同，是否可计算它们总的平均值呢？ • 首先，应确定各层的ORi是否反映了层间的基本差异。用交互的检验（即血压对心肌梗死危险性的影响在不同年龄组是否各异？）可回答是否需要将各个ORi值合并。 • 为了便于计算可信区间和检验是否存在交互作用，将上述两层的数据总结于表下，作一致性检验。
为了检验交互作用，可用下式计算χ2k-1（一致性检验或同质性检验），据此来确定各层的ln ORi与ln OR之间的差异是否超出随机变异所能解释的范围。 （ln ORi− ln OR）2 χ22-1 = ∑────────── Var（ln ORi） （ −0.0513 − 0.4567）2（0.6259 − 0.4567）2 = ────────── + ───────── 0.300 0.100 = 1.15 P > 0.20。
交互作用的检验 • 上例的资料总结如下： 层 ORi ln ORi Var（ln ORi） 1/ Var（ln ORi）= wi 1 0.95 -0.0513 1/9+1/6+1/115+1/73 = 0.300 3.33 2 1.87 0.6259 1/20+1/21+1/596+1/1171 = 0.100 10.00
