兄弟我太懂你这种要统一多图刻度的痛点了！之前做系列meta分析要输出统一风格的森林图，在meta包的forest函数上踩了好几个坑，给你把解决方案捋得明明白白：

一、用meta包原生解决（最省心，不用转格式） #

meta包的forest()函数因为是自己封装的绘图逻辑，不支持基础绘图的xlim/xtick这类参数直接硬传，得用它专门设计的参数来搞：

1. xlim：指定X轴的范围，传一个长度为2的数值向量就行
2. at：自定义刻度的位置，传一个数值向量
3. atlab：对应at的刻度标签（可选，默认就是at的数值）

给你把你的代码修改好直接能用：

rm(list = ls())
library(meta)
study_data <-   data.frame(
  study = c("Study 1","Study 2"),                           
  log_effect= c(-0.2242,0.27110),                           
  n = c(858,736),                           
  sei = c(0.124652,0.1508)
)
meta_results <- metagen(TE = log_effect, seTE = sei,
                        studlab = study,
                        data = study_data,
                        sm = "RR",
                        common = FALSE,
                        random = TRUE,
                        overall = TRUE,
                        overall.hetstat = FALSE,
                        text.random = "Pooled effect")

# 核心修改就在这几个参数！
meta::forest(meta_results,
             leftcols = "study",
             rightcols = c("effect", "ci"),
             col.random = "black",
             col.diamond = "black",
             col.square = "black",
             xlab = "Relative risk ratio",
             # 自定义X轴范围，你可以根据自己的效应值调整
             xlim = c(0.5, 2),
             # 自定义刻度位置，和xlim对应上
             at = c(0.5, 1, 1.5, 2),
             # 刻度标签（可选，不想改的话可以删掉这个参数）
             atlab = c("0.5", "1", "1.5", "2")
)

为啥之前会报错？大概率是你把xlim写成了par(xlim=...)或者参数格式错了（比如只传了一个值），按上面的方式直接传进forest函数就没问题。要是还想调刻度字体大小，加个cex.at=0.8（数值越小字越小）就行。

二、如果一定要用forestploter的解决方案 #

forestploter本身不识别meta包的metagen对象，得自己从meta_results里提取核心数据，转成它能读的数据框格式：

library(forestploter)
# 1. 提取单个研究的效应量（把对数RR转成原始RR）
study_df <- data.frame(
  study = meta_results$studlab,
  rr = exp(meta_results$TE),
  lower = exp(meta_results$lower),
  upper = exp(meta_results$upper)
)
# 2. 提取合并效应
pooled_df <- data.frame(
  study = "Pooled effect",
  rr = exp(meta_results$TE.random),
  lower = exp(meta_results$lower.random),
  upper = exp(meta_results$upper.random)
)
# 3. 合并成forestploter需要的格式
plot_df <- rbind(study_df, pooled_df)
# 4. 绘图
forest(plot_df, 
       est = plot_df$rr, 
       lower = plot_df$lower, 
       upper = plot_df$upper,
       xlab = "Relative risk ratio",
       xlim = c(0.5, 2),
       xticks = c(0.5, 1, 1.5, 2))

最后补一句 #

要是按上面的方法还报错，把具体的错误信息贴出来（比如是“参数不存在”还是“长度不匹配”），我再帮你精准排查！