Stata 操作因为涉及命令,使用过程中可能经常遇到各种各样的小问题,除了善用help,自己去探索解决,不仅趣味丛生更能学会举一反三,下面的文章抛砖引玉,为大家介绍了利用Stata软件制作OR值的森林图并调整横坐标的具体方法。
Stata软件制作森林图,可以采用命令metan或者点击User---Meta-analysis下面的选项。
数据:数据中输入4个变量分别为factor OR CI1 CI2 (分别代表OR值,CI下限、CI上限),并输入数值。
作图如下:
图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。
作图如下,可以看到横坐标标尺不一样:
图2 对数转化后的图形,显示更加清晰
这是采用了横坐标对数转换的图。(即使用ln数据,和选上Effect中的Exponentiate的结果)对比一下,第二个更加直观。关于横坐标转换,标尺不一致的问题大家应该了解了,核心就是对原始数据对数转换,写程序的话,加上eform语句;若使用鼠标,也选对数的数据,同时选上effect中的exp选项。
后续的图表编辑,修改Study id和变量、调整位置等,都是鼠标单击双击点选的问题了,比较繁杂,大家可以在不断摸索中进步。