spss:交互作用指数及其计算 -凯发k8国际首页登录

　　在研究数据的分析过程中，我们经常会遇到交互作用的问题，而且这是一个不太好解释和描述的问题。本文就跟大家一起来讨论一下rothman k.j.等人提出的用来描述交互作用的指标及其计算。rothman k.j.是《modern epidemiology》一书的主编，该书在流行病学领域地位非凡。

　　当两个或多个因素同时作用于一个结局时，就可能产生交互作用，又称为效应修饰作用(effect modification)。当两个因素同时存在时，所导致的效应(a)不等于它们单独效应相加(b c)时，则称因素之间存在交互作用。当a=b c时称不存在交互效应;当a>b c时称存在正交互作用，又称协同作用(synergy);当a

rothman等人在1972发表文章讨论吸烟与饮酒联合暴露对口腔癌的效应，其文章主要数据如下：

　　rothman等人提出的基于相加模型的交互作用指标有：

　　1、交互作用指数(the synergyindex, s)：s = (a-1)/[(b-1) (c-1)]，当s = 1时，说明吸烟饮酒联合效应等于各自单独效应之和，故无交互;当s > 1时，说明正交互;当s < 1时，说明负交互。

　　2、交互作用归因比(attributableproportion of interaction, api)：api = [a-(b c) 1]/a，体现的是总效应中吸烟和饮酒交互作用效应所占的比例。

　　3、纯交互作用归因因比：api’=[a-(b c) 1]/(a-1)，跟api的差别是分母多减去个1，体现的是吸烟和饮酒引起的效应中，两者交互效应所占的比例。

　　4、交互作用超额相对危险度(relativeexcess risk of interaction, reri)：reri = a-(b c) 1，体现的是吸烟与饮酒联合作用与单独作用之和的差值。

　　对于上述交互作用指标的计算有多种方法，包括直接计算法和回归模型估计法。直接法就是手工直接计算，根据上述公式计算，s指数=(9.04-1)/(2.96 3.33-2)=1.87。回归模型法则是将数据放入回归模型进行计算，本文列举如下：

　　1、回归模型法a

　　将上述数据整理成如下格式：

　　在中用case进行加权，然后以cancer为因变量，以smoking、alcohol和smoking*alcohol为自变量计算逻辑回归结果如下：

　　此时s指数=[exp(1.086 1.204-0.089)-1]/[exp(1.086) exp(1.204)-2]=1.87。

　　其他指数类似。

　　2、回归模型法b

　　将上述数据整理成如下格式，也就是设置哑变量的形式，包括四类不饮酒不吸烟，饮酒不吸烟：不饮酒吸烟和饮酒吸烟，以不饮酒不吸烟为参照。

　　在中用count进行加权，然后以y为因变量，以al、sm和add为自变量计算逻辑回归，同时保存协方差阵，结果如下：

　　此时s指数=[exp(2.201)-1]/[exp(1.086) exp(1.204)-2]=1.87。

　　其他指数类似。

　　同样可对上述指标进行可信区间的估计，估算过程需要使用上述我们保存出来的协方差矩阵，具体算法请参阅参考文献中hosmer, 向惠云等人的文章。