RERI相加交互作用-2.excel篇深度评测

上一章，笔者介绍了相加交互的理论，并手动DIY了诸多关系图，小伙伴们以后遇到相加相乘方向相反也不用百思不得其解了。临床流行病常见的结局都是二分类01，那么探讨二分类因素a和二分类因素b交互作用的时候仅仅将a*b纳入方程可能是远远不够，且结论甚至会错误。那实战中如何快速计算出相加交互作用呢？本公众号从浅入深开始介绍RERI的软件实现，顺序是excel → SAS → R。本章节从手动的excel计算开始。

PS: 请回复“相加交互spss代码”, 即可获得方法2代码。

1.excel工具测评

① PMID: 22253321，knol2012年IJE发表。利用原始变量product term(a*b)，或哑变量化后dummy计算，且提供了RERI、AP、S的95%CI (delta法)及p值。强烈推荐！

② PMID: 16119429，andersson2005年EJE发表。利用哑变量化后dummy计算，且提供了RERI、AP、S的95%CI  (delta法)。附表中还提供了spss、SAS等计算代码。

③ PMID: 17726040， knol2007年IJE发表。利用原始变量product term(a*b)仅提供了RERI、AP、S的点估计值。不推荐。

2.数据库变量赋值

假设多风险因素中研究两暴露因素AB交互作用，数据可以按照两种思路进行演化。1.原始变量直接入选，生成a*b变量；2.生成哑变量dum01 dum10 dum11入选。

注：OR00表示AB均无暴露即OR00=1，OR10表示A暴露B无暴露，OR01表示A无暴露B暴露，OR11表示AB均暴露。交互是AB同等重要，效应修饰则是探索A主B辅，常简写为EM，exposuremodifier。本文暂不深入讨论二者区别。

3. excel实战

实战数据为虚拟数据jiaohu.sav。变量a b不分主次，用来演示logistic、cox的相乘相加交互作用。考虑到SPSS与excel的互动性与普及性，二者结合十分适合手动计算的小伙伴。注意，下图中的β1β2β3β4请返回上一章方法学内容。

3.1 方法1，基于哑变量dum：

运行Andersson2005附表提供的SPSS代码，生成新哑变量ind01 ind10 ind11，带入Andersson2005的excel工具。菜单操作是多元logistic内置的协方差矩阵。缺点是不提供P值(当然也可以自己计算)。

3.2 方法2(推荐！)， 基于原始或哑变量：

本例通过上述2个excel工具，可以计算得到OR10、OR01、OR11，INTM点估计值，也顺利得到了delta法计算的RERI、AP、S 95%CI及P值。

4 发散思维

学会了knol与andersson工具，似乎很强，但是似乎又很弱。对于刚入门的小伙伴，excel+spss是最容易实现你目的的方法。久而久之之，你写论文时会发现一会要提供OR01~OR11的95%CI，一会需要提供P值，一会需要进行78个亚组异质性分析。单纯用excel似乎有点费劲，并且可能需要进行cox交互作用分析。是的，讨论到现在所有的excel工具，都只能实现“非条件logistic”!

是不是想起那句话，你知道的越多发现越多的不知。

5 小结

1.推荐笔者改编的SPSS代码+ knol2012 excel表手动计算RERI、AP、S的95%CI与P值。

2.Knol2012或Andersson2005 excel工具，适合非条件logistic的相加交互作用分析。

原创不易，转载请说明来自本公众号。