ggplot2中trans_new函数分面设置x轴范围时的异常问题

问题场景 #

尝试创建分面密度图，其中value1变量通过自定义转换函数log10p_trans处理。未设置x轴范围时图形显示正常，但添加limits = c(0, 0.5)参数后，x轴范围未正确适配，图形出现异常截断。

正常运行的代码 #

library(ggplot2)
library(tidyr)
library(reshape)
library(forcats)
library(scales)

log10p_trans <- trans_new(name = "log10p_trans",
                          transform = function(x) log10(x + 1),
                          inverse = function(x) 10^x - 1)

rex <- data.frame("value1" = rexp(100, rate = 1),
                  "value2" = rexp(100, rate = 2))
rex_pivot <- data.frame(pivot_longer(rex, cols = 1:2))

ggplot(rex_pivot, aes(x = value, y = ..density..)) + 
  geom_density() + 
  facet_wrap(facets = ~name, scales = "free_x") +
  ggh4x::scale_x_facet(
    name == "value1",
    trans = log10p_trans)

出现异常的代码 #

ggplot(rex_pivot, aes(x = value, y = ..density..)) + 
  geom_density() + 
  facet_wrap(facets = ~name, scales = "free_x") +
  ggh4x::scale_x_facet(
    name == "value1",
    trans = log10p_trans,
    limits = c(0, 0.5))

问题原因 #

ggh4x::scale_x_facet的limits参数作用于转换后的数据空间，而非原始数据空间。自定义转换log10p_trans会先对原始x值执行log10(x+1)转换，直接传入原始数据的c(0,0.5)作为limits时，程序会将其视为转换后的值，反向转换回原始数据范围是10^0 -1=0到10^0.5 -1≈1.236，这与预期的0-0.5范围不符，导致图形截断异常。

解决方法 #

方法1：传入转换后的对应范围 #

将原始数据的目标范围通过自定义转换函数计算出转换后的值，作为limits参数传入：

ggplot(rex_pivot, aes(x = value, y = ..density..)) + 
  geom_density() + 
  facet_wrap(facets = ~name, scales = "free_x") +
  ggh4x::scale_x_facet(
    name == "value1",
    trans = log10p_trans,
    limits = c(log10(0 + 1), log10(0.5 + 1))  # 转换后的对应值
  )

方法2：使用oob参数处理超出范围值 #

通过scales包的oob函数（如squish或oob_keep），将超出目标范围的数据压缩到边界或保留，避免图形截断：

ggplot(rex_pivot, aes(x = value, y = ..density..)) + 
  geom_density() + 
  facet_wrap(facets = ~name, scales = "free_x") +
  ggh4x::scale_x_facet(
    name == "value1",
    trans = log10p_trans,
    limits = c(0, 0.5),
    oob = scales::squish  # 把超出范围的值压缩到设定的边界
  )

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