R- Meta-Analysis -如何从clubSandwich函数创建具有稳健估计的森林图

# 安装和加载必要的包
install.packages("clubSandwich")
install.packages("metafor")
install.packages("ggplot2")

library(clubSandwich)
library(metafor)
library(ggplot2)

# 模拟数据
set.seed(123)
n_studies <- 20
true_effect <- 0.5
study_effects <- rnorm(n_studies, true_effect, 0.2)
study_variances <- runif(n_studies, 0.04, 0.16)
data <- data.frame(
  study = 1:n_studies,
  effect = study_effects,
  variance = study_variances
)

# 随机效应模型
model <- rma(effect, vi = variance, data = data, method = "REML")

# 使用clubSandwich计算稳健标准误
robust_se <- clubSandwich::robust_se(model)

# 创建森林图
forest(model, sei = robust_se$SE)