审稿人：影响因素森林图太难看了，你能不能画一个NEJM风格？

妙趣横生统计学

已于 2024-08-14 14:01:41 修改

阅读量1k

点赞数 15

文章标签：人工智能回归数据挖掘算法

于 2024-08-09 11:29:22 首次发布

本文链接：https://blog.csdn.net/weixin_44693403/article/details/141074577

版权

详情请点击下方：

零基础两天学会R语言！！！

森林图目前在相关文献中可以说是非常常见了，不只是亚组分析中，普通的回归分析结果也可以用森林图进行可视化展示，不仅可以帮助我们更好地理解不同变量间的一致性和差异性，新颖的图片还能为文章增色不少。

当下主流绘制森林图的方式还是R语言，其他也有使用GraphPad Prism或者excel，但就R语言而言，一行行的代码参数着实让人头大，不仅需要一定的代码基础，想要画出好看的森林图，参数也要调上好久，效率着实是低了一些......

# 定义一个简单的主题，大家可以随意发挥自己的审美！
tm <- forest_theme(base_size = 10,           # 设置基础字体大小
                   refline_col = "red4",     # 设置参考线颜色为红色
                   arrow_type = "closed",    # 设置箭头类型为闭合箭头
                   footnote_col = "blue4")   # 设置脚注文字颜色为蓝色


# 绘制森林图
p <- forest(result2[,c(1, 2, 9, 6)],   # 选择要在森林图中使用的数据列，这里包括变量名列、患者数量列、绘图要用的空白列和HR（95%CI）列
            est = result2$HR,          # 效应值，也就是HR列
            lower = result2$lower.95,  # 置信区间下限
            upper = result2$upper.95,  # 置信区间上限
            sizes = result2$se,        # 黑框框的大小
            ci_column = 3,             # 在第3列（可信区间列）绘制森林图
            ref_line = 1,              # 添加参考线
            arrow_lab = c("Low risk", "High Risk"),  # 箭头标签，用来表示效应方向，如何设置取决于你的样本情况
            xlim = c(-1, 10),          # 设置x轴范围
            ticks_at = c(-0.5, 1, 3, 5, 7),  # 在指定位置添加刻度
            theme = tm,                # 添加自定义主题
            footnote = "This is the demo data. Please feel free to change\nanything you want.")  # 添加脚注信息
p
# 加列P值
p <- forest(result2[,c(1, 2, 9, 6, 7)],  # 加了P值列哟！选择要在森林图中使用的数据列，这里包括变量名列、患者数量列、绘图要用的空白列和HR（95%CI）列
            est = result2$HR,            # 效应值，也就是HR列
            lower = result2$lower.95,    # 置信区间下限
            upper = result2$upper.95,    # 置信区间上限
            sizes = result2$se,          # 黑框框的大小
            ci_column = 3,               # 在第3列（可信区间列）绘制森林图
            ref_line = 1,                # 添加参考线
            arrow_lab = c("Low risk", "High Risk"),   # 箭头标签，用来表示效应方向，如何设置取决于你的样本情况
            xlim = c(-1, 10),            # 设置x轴范围
            ticks_at = c(-0.5, 1, 3, 5, 7),   # 在指定位置添加刻度
            theme = tm,                  # 添加自定义主题
            footnote = "This is the demo data. Please feel free to change\nanything you want.")   # 添加脚注信息
p

现在郑老师风暴统计平台，可以画出NEJM森林图了！！！！只要三步走，自动绘制出森林图！！！不需要写代码！！！

NEJM: