R语言绘制森林图如何计算OR值

在R语言中，计算OR值可以使用以下代码：

# 假设有两组样本，分别有10个和15个观测值
# 其中第一组中有6个成功事件，第二组中有8个成功事件

# 计算OR值
odds_ratio <- (6/4) / (8/7)  # 或者使用 exp(log(6/4) - log(8/7))

# 输出结果
odds_ratio

在绘制森林图时，可以将每个研究的OR值和其置信区间绘制在图中，例如：

# 假设有3个研究，每个研究的OR值和置信区间如下
studies <- data.frame(
  study = c("Study 1", "Study 2", "Study 3"),
  odds_ratio = c(1.5, 2.0, 0.8),
  lower_ci = c(1.1, 1.3, 0.5),
  upper_ci = c(2.0, 2.7, 1.2)
)

# 绘制森林图
library(forestplot)
forestplot(studies$odds_ratio, 
           lower = studies$lower_ci, 
           upper = studies$upper_ci,
           xlog = TRUE, 
           is.summary = c(FALSE, FALSE, TRUE),
           clip = c(0.1, 10), 
           xlab = "Odds Ratio (95% CI)", 
           col = fpColors(box = "royalblue", line = "darkblue", summary = "royalblue"),
           boxsize = 0.2, 
           txt_gp = fpTxtGp(label = gpar(cex = 0.8), ticks = gpar(cex = 0.7)))

以上代码将绘制出如下的森林图：

![](https://i.imgur.com/xvZvzPb.png