文档建议反馈控制台
首页
学习
活动
专区
圈层
工具
MCP广场
文章/答案/技术大牛
发布
社区首页 >专栏 >R语言绘制森林图

R语言绘制森林图

作者头像
一粒沙
发布于 2019-07-31 09:46:33
发布于 2019-07-31 09:46:33
8.7K0
举报
文章被收录于专栏：R语言交流中心R语言交流中心

在Meta分析中森林图比较常见，其主要是是以统计指标和统计分析方法为基础，用数值运算结果绘制出的图型。它在平面直角坐标系中，以一条垂直的无效线(横坐标刻度为1或0)为中心，用平行于横轴的多条线段描述了每个被纳入研究的效应量和可信区间，用一个棱形(或其它图形)描述了多个研究合并的效应量及可信区间。它非常简单和直观地描述了Meta分析的统计结果，是Meta分析中最常用的结果表达形式。

森林图的类型主要包括以下两种：

1. 二值变量的森林图

当研究对象为二值变量（如发生与不发生）时，采用RR和OR作为统计学指标。

RR（risk ratio,危险比率）=发生某一事件的人数/观察的总人数。

OR（odds ratio，比值比）=发生某一事件的人数/未发生某事件的人数。

2. 连续性变量的森林图

分析指标是连续变量，也称数值变量，可选择加权均数差(WMD)或标准化均数差(SMD)为合并统计量。

加权均数差（weighted mean difference）是对每个研究的均数进行加权处理。

大样本研究，由于效应估计值更为准确，权重较大；

小样本研究，由于估计值不够准确，权重较小。

标准化均数差（standardise mean difference）在每一试验中以不同的测量单位对同一结局描述时，需要进行标准化处理。

森林图中横短线与中线相交表示无统计学意义；横线在左侧说明结局弱于竖线结局；在右侧说明强于竖线结局。最后以菱形所在位置得到总体的评价结果。

以上就是森林图的理论知识。接下来我们介绍其是如何在R语言中实现的。

首先我们需要导入R包forestplot。具体的安装载入不再赘述。

下面我们看下其主要的函数forestplot。

其中主要的参数：

tabletext：主要是以矩阵或者list形式将数据导入函数，最好以矩阵，因为数据一般都是矩阵的。

Mean,lower,uper一般会以list形式直接一次性导入，当然如果数据是一个一个向量还是分开的方便。

is.summary主要的功能是让表格的每一行字体出现差异，从而区分表头。其值主要用TRUE/FALSE进行差异化分配。

Graph.pos主要是定位森林图所在的位置。通过数字来确定为第几列。

hrzl_lines以list形式设置表中线的类型以及长度。Eg：“3”=gar(lwd=1,columns=1:4,col=’red’)意思就是宽度为1，线段延伸至第四列。Col指的颜色。

clip主要x轴的最大最小值。

Col指的其中横线以及点的颜色。eg:col=fpcolors(box=’royblue’,line=’darkblue’, summary=’royblue’, hrz_lines=’red’)设置其中这些部分的的颜色。

txt_gap主要是设置表格中文本的字体

最后我们看两个样例程序：

代码语言：javascript
复制
library(forestplot)
代码语言：javascript
复制
# Cochrane data from the 'rmeta'-package
代码语言：javascript
复制
cochrane_from_rmeta <-
代码语言：javascript
复制
  structure(list(
代码语言：javascript
复制
    mean  =c(NA, NA, 0.578, 0.165, 0.246, 0.700, 0.348, 0.139, 1.017, NA, 0.531),
代码语言：javascript
复制
    lower =c(NA, NA, 0.372, 0.018, 0.072, 0.333, 0.083, 0.016, 0.365, NA, 0.386),
代码语言：javascript
复制
    upper =c(NA, NA, 0.898, 1.517, 0.833, 1.474, 1.455, 1.209, 2.831, NA, 0.731)),
代码语言：javascript
复制
    .Names =c("mean", "lower", "upper"),
代码语言：javascript
复制
    row.names =c(NA, -11L),
代码语言：javascript
复制
    class ="data.frame")
代码语言：javascript
复制
代码语言：javascript
复制
tabletext<-cbind(
代码语言：javascript
复制
  c("", "Study", "Auckland", "Block",
代码语言：javascript
复制
    "Doran", "Gamsu", "Morrison", "Papageorgiou",
代码语言：javascript
复制
    "Tauesch", NA, "Summary"),
代码语言：javascript
复制
  c("Deaths", "(steroid)", "36", "1",
代码语言：javascript
复制
    "4", "14", "3", "1",
代码语言：javascript
复制
    "8", NA, NA),
代码语言：javascript
复制
  c("Deaths", "(placebo)", "60", "5",
代码语言：javascript
复制
    "11", "20", "7", "7",
代码语言：javascript
复制
    "10", NA, NA),
代码语言：javascript
复制
  c("", "OR", "0.58", "0.16",
代码语言：javascript
复制
    "0.25", "0.70", "0.35", "0.14",
代码语言：javascript
复制
    "1.02", NA, "0.53"))

forestplot(tabletext,

graph.pos =4,

hrzl_lines =list("3"=gpar(lty=2),

"11" =gpar(lwd=1, columns=c(1:3,5), col ="#000044"),

"12" =gpar(lwd=1, lty=2, columns=c(1:3,5), col ="#000044")),

cochrane_from_rmeta,new_page =TRUE,

is.summary=c(TRUE,TRUE,rep(FALSE,8),TRUE),

clip=c(0.1,2.5),

xlog=TRUE,

col=fpColors(box="royalblue",line="darkblue",summary="royalblue",hrz_lines ="#444444"))

代码语言：javascript
复制
data(HRQoL)
代码语言：javascript
复制
clrs <-fpColors(box="royalblue",line="darkblue", summary="royalblue")
代码语言：javascript
复制
tabletext <-
代码语言：javascript
复制
  list(c(NA, rownames(HRQoL$Sweden)),
代码语言：javascript
复制
       append(list(expression(beta)), sprintf("%.2f", HRQoL$Sweden[,"coef"])))
代码语言：javascript
复制
tabletext <-tabletext[,1]
代码语言：javascript
复制
forestplot(tabletext,
代码语言：javascript
复制
            fn.ci_norm =c(fpDrawNormalCI, fpDrawCircleCI),
代码语言：javascript
复制
            boxsize = .25, # We set the box size to better visualize the type
代码语言：javascript
复制
            line.margin = .1, # We need to add this to avoid crowding
代码语言：javascript
复制
            mean =cbind(HRQoL$Sweden[, "coef"], HRQoL$Denmark[, "coef"]),
代码语言：javascript
复制
            lower =cbind(HRQoL$Sweden[, "lower"], HRQoL$Denmark[, "lower"]),
代码语言：javascript
复制
            upper =cbind(HRQoL$Sweden[, "upper"], HRQoL$Denmark[, "upper"]),
代码语言：javascript
复制
            clip =c(-.125, 0.075),
代码语言：javascript
复制
            col=fpColors(box=c("blue", "darkred")),
代码语言：javascript
复制
            xlab="EQ-5D index")
代码语言：javascript
复制
forestplot(tabletext,
代码语言：javascript
复制
            fn.ci_norm =c(fpDrawNormalCI, fpDrawCircleCI),
代码语言：javascript
复制
            boxsize = .25, # We set the box size to better visualize the type
代码语言：javascript
复制
            line.margin = .1, # We need to add this to avoid crowding
代码语言：javascript
复制
            mean =cbind(HRQoL$Sweden[, "coef"], HRQoL$Denmark[, "coef"]),
代码语言：javascript
复制
            lower =cbind(HRQoL$Sweden[, "lower"], HRQoL$Denmark[, "lower"]),
代码语言：javascript
复制
            upper =cbind(HRQoL$Sweden[, "upper"], HRQoL$Denmark[, "upper"]),
代码语言：javascript
复制
            clip =c(-.125, 0.075),
代码语言：javascript
复制
            lty.ci =c(1, 2),
代码语言：javascript
复制
            col=fpColors(box=c("blue", "darkred")),
代码语言：javascript
复制
            xlab="EQ-5D index")
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2018-07-03，如有侵权请联系 cloudcommunity@tencent.com 删除
html
编程算法

本文分享自 R语言交流中心 微信公众号，前往查看

如有侵权，请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ，欢迎热爱写作的你一起参与！

html
编程算法
评论
登录后参与评论
0 条评论
热度
最新
登录 后参与评论
推荐阅读
领券

腾讯云开发者

扫码关注腾讯云开发者

扫码关注腾讯云开发者

领取腾讯云代金券

热门产品

热门推荐

更多推荐

Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有 

深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 粤公网安备44030502008569号

腾讯云计算（北京）有限责任公司 京ICP证150476号 |  京ICP备11018762号

问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档

Copyright © 2013 - 2026 Tencent Cloud.

All Rights Reserved. 腾讯云 版权所有

登录 后参与评论