二分类变量的meta分析，为什么OR值是负数？

近段时间，有不少小伙伴给我留言，为什么二分类变量的meta分析，得到的OR值是负数，这个结果怎么解释？ 我们都知道，OR是“Odd ratio”的缩写，中文名称是“比值比”，是两组计数资料频数分布的比值差异，其取值范围是0到+∞。可是，为什么软件会出现OR<0的结果呢？

 
这个事情，看上去是软件的问题，但仔细一想，问题出在我们自己身上。 先说软件的问题，现在很多meta分析教程都是基于低版本软件的，例如Stata 11/12，对于这两个版本，二分类变量的meta分析，若使用OR值为效应值，计算得到的是OR值。 以下面这些数据为例，用OR (95%CI)评估两组死亡风险的差异

 
通过菜单或命令行操作，得到随机效应模型的合并结果，森林图以OR (95%CI)展示结果。

 
数据窗口，软件也生成了两个临时变量_ES和_selogES，也就是OR和selogOR。

 
如果接下来我们要做敏感性分析或发表偏倚检验，需要计算logOR和selogOR，那就要利用以下代码将_ES转换为logOR：gen logOR=log(_ES) 可是，当我们用Stata14.0操作的时候，出来的数据是这样的！两个临时变量，_selogOR变成了_seES，但数据完全一样；而_ES的变量名没变，可数据截然不同，甚至出现了负数。

 
这个时候再进行对数转换，就会出现部分数据缺失的情况。为什么有缺失？负数没法进行自然对数转换，这是中学数学的知识！

 
那是软件出错了吗？软件也没错，它只是在版本升级的同时，做了一些更新，将二分类变量的meta分析结果以对数值的形式展示，也就是此时的_ES实际上就是11/12版本的logOR。
我们可以验证一下，下面的logOR是否与上图的_ES (保留小数点后7位)一模一样？

 
OK，第一个情况解决了。如果你是使用Stata14.0完成了二分类变量的meta分析，得到_ES后就不要再进行对数转换了，因为它就是OR的自然对数。 还有一种更有趣的情况，使用Stata16的新菜单“meta”实现二分类变量的森林图，如果菜单选择有误，得到的森林图都是以对数值显示结果。