森林图SAS代码

森林图

前面我们在这篇文章里介绍过亚组分析森林图。

我们当时是借助宏来实现森林图的绘制，但是宏没有适配展示生存分析的中位生存时间、有效事件数这些有用的信息，使用起来非常不灵活；另外有时候，我们是用R、python跑出来的亚组分析结果，如果想保持一样的图片出版风格，也不容易统一。如，我们想在文章中加入下图风格的森林图：

这样SAS和R出来的风格很难一致。如果我们还想额外增加一列subgroup人数的count，这个时候就非常需要改动起来非常灵活的SAS代码了。

解决

出于以上的种种考虑，刚好自己最近在做相关的文章分析，我们特别定制了改动非常容易的SAS代码，想着尽可能地保持高灵活高度统一的森林图风格，完整的代码已经上传进我的星球，欢迎感兴趣的或者有实际需求的同学自行下载查看使用。

我自己使用起来的感受，还是非常不错的，可以像上面第二张森林图一样，灵活地额外增加列的信息，如各个亚组的有效事件数count以及percent、交互项p值以及脚注。另外，由于我们不规定森林图之前的分析，所以，可以自由地使用SAS、R、python的回归结果来作图。图片格式可选择有pdf、png。

Did you find this page helpful? Consider sharing it 🙌