交互作用的检验

1、交互作用与效应修饰交互作用与效应修饰 交互作用( interaction) 指两个或多个因素相互依赖发生作用而产生的一种效应。生物学上是指两个或多个因素互相依赖发生作用而引起疾病或预防和控制疾病。即生物学交互作用(biological interaction) . 研究资料中各研究因素间存在交互作用时，说明各研究因素的作用不是独立的，一个因素水平改变时与它有交互作用的因素之效应也将发生改变 。 交互作用是描述资料本身所表现出来的关于两个或多个暴露因素联合效应的一个统计学术语，并不表示一定具有生物学意义。有交互作用可能有生物学意义，也可能没有生物学意义。效应修饰作用 当暴露因素按第三变量分层后估

2、计暴露在每一层中与疾病的联系强度时，效应修正被定义为暴露因素在各层中与疾病的联系强度(测量的效应)因第三变量的存在情况不同而大小不同。 该第三变量称为效应修正因素( effect modification factor）。 第三变量在一项研究中是否成为效应修正因素，取决于选用判断暴露和疾病之间联系的指标是用率差( rate difference , rd)还是用率比( rate ratio , rr) 。 又称效应变异（effect variation），或效应不一致性或异质性（heterogeneity of effect）。 例如，性别与髋骨骨折的联系受到年龄的修饰； 女性髋骨骨折危险性为

3、男性的2至3倍 年轻组：男性髋骨骨折危险性高于女性 年老组：女性髋骨骨折危险性明显高于男性 说明性别与骨折的联系被年龄修饰。 体重与乳癌的联系受到绝经状态的修饰； 肥胖与乳癌之间的联系被绝经状态修饰， 肥胖是绝经后妇女乳癌的危险因素， 肥胖不是绝经前妇女乳癌的危险因素。 孕妇吸烟与婴儿低出生体重的联系 吸烟对出生低体重的有害影响在高龄孕妇较低龄孕妇更为明显。 母亲年龄是修饰因素。 eb病毒感染与非洲儿童淋巴瘤的联系 生命早期感染eb病毒与非洲儿童burkitt淋巴瘤的联系多见于幼年生活在疟疾广泛流行区的儿童，少见于疟疾非流行区。 当地普遍存在疟疾是一种效应修饰因素。效应修饰与混杂n 混杂是一种

4、偏倚，是研究者希望避免的，或在必要时希望从资料中消除的。n 效应修饰是对效应本身的详细描述，它是一种需要报告的发现，而不是一种需避免的偏倚。n 流行病学分析总目的是消除混杂，并发现和描述效应修饰。n 混杂是否存在取决于研究设计。n 混杂源于选择研究对象的源人群中混杂变量与研究因素的相关性。因此，在选择研究对象时对研究对象进行适当限制，可以防止一个变量成为混杂因素。n 效应修饰并不取决于研究设计，是一种与研究设计无关的自然现象，是研究欲尽可能揭示和描述的现象，它超出了恒定的研究设计概念。n 效应修饰是一种恒定的自然现象，这是一种相对（非绝对）的概念，故其并不与任何生物学特性相对应。n 在研究效应

5、修饰时必须规定欲测量的是何种效应。n 广义上讲，效应修饰是指对一种效应的修饰，但并未指明修饰何种效应测量。流行病学中有两种常用的危险性效应测量，即率差和率比，以及其它一些不常用的测量。n 如果不规定修饰何种测量（率差或率比），效应修饰的概念就太模糊，使人无法描述其特性。暴露组与非暴露组之间的发病率比发病率比在不同年龄组是恒定的。但其发病率差发病率差则随年龄而增加。以率差为指标，年龄是效应修饰因素。暴露组与非暴露组之间的发病率差发病率差在不同年龄时是恒定的。对率比率比这个效应来说，年龄修饰暴露效应。n 如果每天喝5杯咖啡，可使男性心肌梗死发生率增加40%，（率比rr=1.4），但使女性心肌梗死发

6、生率只增加10%（rr=1.1），这种情况称作性别对rr的修饰（或rr的性别变异或性别的rr不一致性），性别即为喝咖啡致心肌梗死rr的修饰因素。例如n 又例如，每天喝5杯咖啡，可使男性心肌梗死发生率增加到400/10万人年，但只使女性心肌梗死发生率增加到40/10万人年。这种情况称为性别对率差（ar）的修饰，性别即为喝咖啡致心肌梗死ar的修饰因素。 再一例子，如果每天喝5杯咖啡，无论对男性或女性，都可使心肌梗死发生率增加22%（rr=1.22），这种情况称为性别间的rr一致性（homogeneity），也即性别不存在效应修饰。 需注意的是，效应修饰和效应不一致性并非效应的固有特性，而仅是效应测

7、量方法的一个特征。 例如，每天喝5杯咖啡，可使男性心肌梗死发生率从1000/10万人年增加到1220/10万人年，而使女性心肌梗死发生率从400/10万人年增加到488/10万人年，男女性的rr均=1.22，呈现一致性。 男性的ar=220/10万人年，女性ar=88/10万人年，男女性的ar不一致或ar被性别所修饰。 这个例子说明，效应修饰不一定等于生物学的交互作用（interaction）（如协同或拮抗作用）。效应修饰作用的分析效应修饰作用的分析 效应修饰可采用分层的方法进行分析，估计每层的效应，并进行比较。 但如果每层观察对象的数目太少，使每层的效应不稳定，尤其在调整了混杂因素后。这样，

8、由于随机误差，层间的效应估计可能波动很大。 检验分层资料不一致性的统计学方法，在许多情况下检验效率很低，用常规的显著性水平（如=0.05），多数的层间不一致性无法检出。用数学模型拟合的方法可以解决其中的一些问题。但由于混杂、测量误差引起的偏倚在各层间不同，对效应修饰无论采用观察或数学模型拟合的方法，仍会有一些偏倚。交互作用模型 相加模型(additive model) 相加模型假定若交互作用不存在时，两个或两个以上因子共同作用于某一事件时，其效应等于这些因子单单独作用时的和，有时称之具有可加性。 r11-r00=(r10-r00)+(r01-r00) 相乘模型( multiplicative

9、model) 相乘模型假定若交互作用不存在时，两个或两个以上因子共同作用于某一事件时，其效应等于这些因子单独作用时的积。 r11/r00=(r10/r00)(r01/r00)相加模型：r11r0040-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不同， 表明其中存在着混杂。将资料按某种特性分层(如可疑的效应修正因子)

10、 ，计算各层相对危险度，若层别相对危险度不等，说明资料中存在着交互作用。 例如：某前瞻性研究观察某团体中甲苯磺丁脲治疗组和安慰剂组各种原因死亡数情况。交互作用的识别 交互作用识别的一般过程 在判断交互作用存在与否之前，首先要明确所研究的因素与事件之间是否存在统计学联系。如果有联系，接下来看这一联系是否由偏倚或混杂所致。若存在偏倚或混杂因素，采用适当的方法加以改进和调整后再分析交互作用。叉生分析 叉生分析：x 、y 为两个危险因素, orx 、ory 、orxy分别表示仅有x ,仅有y 和x 、y 都存联合作用模型。若:orxy orx + ory - 1, 则认为x 、y 间可能存在交互作用。

11、但要作统计学推断,则需作2 检验。 叉生分析仅可以给人们提示某两个因素间是否可能存在交互作用,并且叉生分析的统计学检验方法本身尚存在一定的缺陷;另外,在用叉生分析法进行交互作用分析时,未考虑其他不参加叉生的危险因素的作用。因此, 可在多因素logistic 回归模型中加入交互作用项,以平衡其他因素的影响。 分层分析 分层分析是比较经典的识别交互作用的方法。可以通过比较按照可疑交互因素分层后层间的效应测量值相对危险度(rr) 或率差(rd) 来判断是否产生交互作用。如果各层之间的效应测量值rr 或rd 不同，则可能存在交互作用。但是鉴于各层rr 和rd变异可能是机遇所致，因此必须进行统计学检验。

12、 分层分析 如果效应一致性的假设不合理，那么按第三变量分层后，各层中效应不一致的情况，不应采用简单合并估计效应的方法。 如果第三变量各层间的效应不同时，即存在效应修饰，资料分析和表达的关键是描述分层因素（即第三变量）是如何修饰效应的。 评价效应修饰常依赖于机械地应用统计学检验方法。 建立在统计学效应评价基础上的流行病学交互作用，有时是难以捉摸的，用纯粹机械的方法会使问题变得模糊不清。 用统计学方法评价效应修饰时应小心。 由于检验的备择假设是非特异性的，一般用于效应修饰检验的普通统计学方法，效率较低。因此，“无显著性的”p值较难正确解释。分层分析时的注意问题 由于选择偏倚、分类错误、混杂、以及其

13、它偏倚和致病效应，很少期望用任何尺度测量的效应都会一样精确。 一致性效应的无效假设常只是一项统计学设计，这类统计学设计应只能被看作是近似真实情况，仍存有疑问。分层分析的步骤 第一步第一步是按第三变量分层后，检查各层的效应估计。即使基本参数一致，层间估计值还会有某些随机变异，但过大的随机变异或明显的非随机变异可能就是效应修饰的证据。 关于效应修饰的判断不应只限于资料的一些表面现象，如果有效应修饰的话，第二步第二步还应在评价过程中结合以往研究的其它知识或一般生物学知识予以考虑。 第三步第三步，用多种统计学检验方法来检验层间一致性。 不同的危险性指标，如率差和率比，其效应修饰的测量要求有不同的评价方

14、法，比的测量一致性意味差的测量效应修饰，反之，差的测量一致性可能意味比的测量效应修饰。 统计学检验效应一致性的无效假设，一般有两类，一是按一致效应的直接加权合并估计，另一是按最大似然估计。 效应修饰只能按某一种效应测量指标来进行描述。 如果，按率差或按率比，都没有效应修饰，除非未暴露组的疾病率与可能的效应修饰因素之间没有联系，其仍可能在其它测量指标（如smr）上存在效应修饰。 处理效应修饰的基本分析方法是分层分析分层分析。 将资料按第三变量进行分层分析有两个目的，一个是评价和消除混杂，另一个是评价和描述效应修饰。 分层分析是处理这两种问题的好方法。以色列缺血性心脏病研究以色列缺血性心脏病研究

15、心肌梗死与收缩血压关系的病例对照研究 心肌梗死 对照 计 收缩血压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 61

16、6 收缩血压 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

17、.00，or2=1.87反映了年龄 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据此，尚缺乏证据说明表2中两层的ori不是来自同一总体，因此，可以认为他们来自同一总体，总估计orw值1.58（见下式）就可以被解释为收缩血压升高与发生心肌梗死危险性之间的总的关系。 3.33（ 0.0513）+ 10.0

18、0（0.6259）ln orw = = 0.4567 3.33 + 10.00orw = e0.4567 = 1.58（hauck加权法）p上例只有15例心肌梗死病人大于60岁，各层ori中，所用的检验方法检出其差异大于机遇可能性（抽样误差）的能力（检验效率）十分有限。p检出样本资料差异的能力称检验效率。p其它因素相同时，样本越小，检出差异的检验效率（1 ）也越小。上例，由于缺乏交互证据，利用前述计算的总估计值orw是合理的，并可按下式计算orw的的95%可信区间：可信区间：ln orw = 0.4567se ln orw = 1/wi = 1/（3.33+10.00） = 0.273995%

19、的ln orw的可信区间 = ln orw 1.96（1/wi） = 0.4567 1.96（0.2739） = 0.08010.9935orw = e-0.0801e0.9935 = 0.922.70。直接加权合并估计直接加权合并估计 直接加权合并估计检验的基本原则是将各层的估计与合并的估计进行比较，用其差值平方除以层估计效应的方差。 对效应的差测量对效应的比测量 需进行对数转换。 ln(ri) ln(rw)22n-1 = var（ln ri） 上式中ri = 第i层的比效应估计，rw = 直接合并的比效应。 吸烟和不吸烟英国男医生冠心病死亡专率吸烟和不吸烟英国男医生冠心病死亡专率 英国男医

20、生的年龄别冠心病死亡专率 吸烟者（1） 不吸烟者（0） 年龄 rr 死亡数 人年数 率 死亡数 人年数 率 a1i n1i r1i b0i n0i r0i3544 32 52407 .00061 2 18790 .00011 5.744554 104 43248 .00240 12 10673 .00112 2.145564 206 28612 .00720 28 5710 .00490 1.476574 186 12663 .01469 28 2585 .01083 1.367584 102 5317 .01918 31 1462 .02120 0.90 计 630 142247 .0044

21、3 101 38220 .00264 吸烟和不吸烟英国男医生冠心病死亡专率吸烟和不吸烟英国男医生冠心病死亡专率按上表资料计算的层发病率差估计、以及资料合并的方差和权数 发病率差的 方差 权数年龄 估计(1/万) var(i) (wi ) （1）（3） (1) (2) (3) （4）3544 5.04 1.73 57.8 291.3124554 12.8 16.1 6.21 79.4885564 23.0 111 0.90 20.76574 38.6 585 0.19 7.3447584 20.2 1810 0.06 1.212 计 6.10* 65.16 397.632 * 合并的发病率差估计

22、 = （1）（3） / （3），wi = 1 / var（i）同质性检验2 = 8.38，=4，p=0.08，说明各层率差无明显差异（在=0.05水平上）。 down氏综合征与母亲孕前使用杀精药氏综合征与母亲孕前使用杀精药 down氏综合征与母亲孕前使用杀精药病例对照研究母亲年龄（岁） 35 35 计 杀精药 + 计 + 计 + 计down氏 3 9 12 1 3 4 4 12 16综合征对照 104 1059 1163 5 86 91 109 1145 1254计 107 1068 1175 6 89 95 113 1157 1270down氏综合征与母亲孕前使用杀精药氏综合征与母亲孕前使用杀精药 各层比值比及其对数转换、方差、权数母亲年龄 比值比 比值比对数 比值比对数的方差 权数 （岁） （ori） （ln ori） var（ln ori） wi 35 3.4 1.22 0.46 2.20 35 5.7 1.75 1.54 0.65 * wi = 1 /（1/a + 1/b + 1/c + 1/d），同质性检验2 = 0.14， = 1，p = 0.7，说明各层的发病率比无明显差异。down氏综合征与母亲孕前使用杀精药氏综合征与母亲孕前使用杀精药 直接加权合并估计直接加权合并估计的假设是各层内观察频数