Stata如何制作OR值的森林图并调整横坐标

2016-08-02 12:34 来源:丁香园 作者:guiminhuang
字体大小
- | +

Stata 操作因为涉及命令,使用过程中可能经常遇到各种各样的小问题,除了善用help,自己去探索解决,不仅趣味丛生更能学会举一反三,下面的文章抛砖引玉,为大家介绍了利用Stata软件制作OR值的森林图并调整横坐标的具体方法。

Stata软件制作森林图,可以采用命令metan或者点击User---Meta-analysis下面的选项。

数据:数据中输入4个变量分别为factor OR CI1 CI2 (分别代表OR值,CI下限、CI上限),并输入数值。

作图如下:

12060455.jpg
图1采用原始数据OR,CI1,CI2作图的话,同时不勾选Exponentiate(图的横坐标标尺是一致的、等距的)

采用原始数据OR,CI1,CI2作图,可以发现age/B/C这三组数据的CI比较接近,近似一个点,在图上无法清除的显示出来。该森林图的标记线在X=0处,如果要使X=1,可以在Graph Opts增加一个null(1)即可;若不写,此处默认是null(0),即X=0处。下面展示使用2种方法实现横坐标调整:

方法一

1. 采用命令metan,将数据进行对数转换,产生新变量,程序如下:
gene lnor=ln(OR)
gene lnci1=ln(CI1)
gene lnci2=ln(CI2)

2. 使用eform命令:Metan lnor lnci1 lnci2, eformlabel(namevar=factor)。

3. 在X轴上添加坐标数值,常常0.5,1.0,10也要写上,用xlabel语句,如xlabel(0.5,1.0,10),数值可随意增加和修改,用逗号分开即可;X轴的坐标尺采用xtick命令,xtick只显示坐标轴下面的小线条即坐标尺,不出现数值,如xtick(2,3,4,5,6,7,8,9)。

4.总结来说,做OR的森林图程序命令就是metan lnor lnci1 lnci2,eformlable(namevar=factor) xlable(o.5,1.0,10) xtick(2,3,4,5,6,7,8,9) nowt nooverallnobox。

方法二

点击User—Meta-analysis—of binary and continuous进行操作:1. 对数ln转换(用excel或者前面的gene命令语句)。2. 做危险因素OR的森林图,常常不需要最后一行stata直接给出的overall,就用nooverall,即overall那一行的OR不出现;不需要出现最右边加权项的值,则用noowt;图中的OR点上有时候去会看到灰色的大方块,这个灰色方块表示权重的大小,常常不好看,若不需要显示则用nobox。3. 在effect选上Exponentiate。

76514611.jpg

45923515.jpg

20917901.jpg

作图如下,可以看到横坐标标尺不一样:

95009802.jpg
图2 对数转化后的图形,显示更加清晰

这是采用了横坐标对数转换的图。(即使用ln数据,和选上Effect中的Exponentiate的结果)对比一下,第二个更加直观。关于横坐标转换,标尺不一致的问题大家应该了解了,核心就是对原始数据对数转换,写程序的话,加上eform语句;若使用鼠标,也选对数的数据,同时选上effect中的exp选项。

后续的图表编辑,修改Study id和变量、调整位置等,都是鼠标单击双击点选的问题了,比较繁杂,大家可以在不断摸索中进步。

编辑: 任悠悠

版权声明

本网站所有注明“来源:丁香园”的文字、图片和音视频资料,版权均属于丁香园所有,非经授权,任何媒体、网站或个人不得转载,授权转载时须注明“来源:丁香园”。本网所有转载文章系出于传递更多信息之目的,且明确注明来源和作者,不希望被转载的媒体或个人可与我们联系,我们将立即进行删除处理。同时转载内容不代表本站立场。

  • App下载