数据可视化——R语言forestplot包绘制组间差异对比图(森林图forest plot)
概述:使用R语言中的forestplot包绘制组间差异对比图。forestplot包本来用于绘制森林图,此处笔者将此用于绘制组间差异对比图,异曲同工,为另一篇博文:数据可视化——R语言ggplot2包绘制组别间指标差异对比图(箱形图及误差条图)提供了另一种实现方案。
森林图(forest plot)常用于Meta分析,可用于表达统计指标的效应量和置信区间。本文中,笔者将森林图中的效应量用于表示各组间的属性指标差异的中值,森林图中的置信区间用于表示各组间的属性指标差异的四分位距的范围。
使用工具:R语言中的forestplot包
如果没有安装forestplot包,输入install.packages(“forestplot”)进行安装。
以下示例中的数据模拟产生:共5个模型,每个模型都有一个属性指标,且每个模型的每个属性指标都进行了100次重复测量。绘制差异对比图时需要额外的计算,即每两两模型之间指标属性相减求出差异,然后寻找差异值的上、下四分位数和中值即可。
示例代码如下:
rm(list=ls()) #清除工作区
library(forestplot)
#生成模拟数据
# 共5个模型,每个模型都有一个属性指标,且每个模型的每个属性指标都进行了100次重复测量
group1 <- rnorm(100)+2
group2 <- rnorm(100)+1
group3 <- rnorm(100)+0.5
group4 <- rnorm(100)
group5 <- rnorm(100)-0.5
Data <- data.frame(group1=group1,group2=group2,group3=group3,group4=group4