如何在R的forestplot中绘制极高OR值？解决x轴截断问题

问题出在forestplot默认的x轴刻度生成逻辑上——当你扩展clip上限后，它不会自动生成覆盖高值范围的对数刻度。要同时满足展示221的置信区间上限和保留低值区分度，需要手动指定对数刻度的刻度值，结合调整clip参数：

修改后的代码 #

# Cochrane data from the 'rmeta'-package
base_data <- tibble::tibble(mean  = c(19.92 , 41.46, 11.67, 11.69, 25.44, 105.89, 145.45),
                            lower = c(17.09, 34.70, 9.04, 10.92, 19.78, 67.40, 95.64),
                            upper = c(23.22, 49.54, 15.07, 12.51, 32.73, 166.37, 221.22),
                            study = c("Autism", "Conduct Problems", "Tic Disorder", "ADHD",
                                      "OCD", "Schizophrenia", "Manic Bipolar"),
                            OR = c("19.92" , "41.46", "11.67", "11.69", "25.44", "105.89", "145.45"))

base_data |>
  forestplot(labeltext = c(study, OR),
             clip = c(0.1, 221),  # 扩展clip到数据的最大上限
             xlog = TRUE,
             xticks = c(0.1, 1, 10, 50, 100, 200)) |>  # 手动指定对数刻度的刻度值
  fp_set_style(box = "royalblue",
               line = "darkblue",
               summary = "royalblue") |> 
  fp_add_header(study = c("", "Study"),
                OR = c("", "OR")) |>
  fp_append_row(mean  = 60.22,
                lower = 41,
                upper = 83,
                study = "Summary",
                OR = "60.22",
                is.summary = TRUE) |> 
  fp_set_zebra_style("#EFEFEF")

关键调整说明 #

1. clip参数：直接设置为c(0.1, 221)，完全覆盖数据的置信区间范围（从0.1到最高221.22），确保所有数据点和置信区间都能完整显示。
2. xticks参数：手动指定包含低值和高值的刻度点，比如c(0.1, 1, 10, 50, 100, 200)。这些值在对数轴上会均匀分布，既保留了0.1到10区间的刻度区分度，又能展示到200（接近你的221上限）的高值刻度。
3. 若需要更精细的刻度，可根据数据调整xticks的数值，比如添加25、150等，确保刻度符合你的展示需求。

内容的提问来源于stack exchange，提问作者Riley