在数据分析和可视化中,森林图(Forest Plot)是展示多组数据比较结果的有效工具,特别是在医学研究中评估风险比或危险比(Hazard Ratio, HR)时非常常见。本文将探讨如何在R语言中使用ggplot2包来创建一个精美的森林图,并解决一些常见的问题。
问题背景
假设我们有一个数据集final_cox_dat,包含了不同组别(Group A, Group B, Group C)在加权和非加权模型下的估计值及其置信区间。以下是数据样本:
library(ggplot2) final_cox_dat <- data.frame( comparator = rep(c("Group A", "Group B", "Group C"), each = 2), estimate = c(1.2, 1.5, 0.8, 0.9, 1.4, 1.7), conf.low = c(1.0, 1.3, 0.7, 0.8, 1.3, 1.5), conf.high = c(1.5, 1.7, 0.9, 1.0, 1.6, 1.9), model = rep(c("Weighted", "Unweighted"), times = 3) )创建森林图
首先,我们使用ggplot2来绘制基本的森林图:
forest_plot <- ggplot(final_cox_dat, aes(x = estimate, y