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。
作图如下,可以看到横坐标标尺不一样:
这是采用了横坐标对数转换的图。(即使用 ln 数据,和选上 Effect 中的 Exponentiate 的结果)对比一下,第二个更加直观。关于横坐标转换,标尺不一致的问题大家应该了解了,核心就是对原始数据对数转换,写程序的话,加上 eform 语句;若使用鼠标,也选对数的数据,同时选上 effect 中的 exp 选项。
后续的图表编辑,修改 Study id 和变量、调整位置等,都是鼠标单击双击点选的问题了,比较繁杂,大家可以在不断摸索中进步。