如何在同一 X 轴的同一面板上绘制具有两个不同 y 轴范围的点?

mgh*_*haz 2 r axes ggplot2 multiple-axes

我正在尝试使用 ggplot 创建一个散点图,该散点图共享一个 X 轴,但具有一个具有两个不同比例的 Y 轴。

Y轴的底部有三个刻度,从0%到0.1%,然后是0.1%到1%,最后是10%的规则间隔。

这里的一个例子:

在此处输入图片说明

有没有办法使用 ggplot 在 R 中产生这样的东西?我会修改轴吗?在同一个面板上叠加多个图?或者是其他东西?

Dav*_*son 5

通常不推荐在 ggplot2 中使用不连续的轴,原因在这个问题 中探讨。Hadley Wickham(ggplot2 的创建者)在这里解释道:

我不是这种类型的显示器的忠实粉丝,因为我认为它在视觉上会扭曲。我认为显示两个图更合适 - 一个是所有数据,一个是小值。这样你就可以看到大值在多大程度上支配了小值。

不过,确实有可能!您必须创建自定义轴转换,这不适合胆小的人。

这是一个例子。假设我们有y对数正态分布的数据。

set.seed(20)
dat <- data.frame(x = c(0, rnorm(50)), y = c(0, exp(rnorm(50, -2, 1.5))))
ggplot(dat, aes(x, y)) + geom_point()
Run Code Online (Sandbox Code Playgroud)

原来的

许多点靠近底部:假设我想将所有小于 1 的值放在对数刻度上,将大于 1 的值放在线性刻度上。因此,我创建了一个名为 的自定义转换,combine_trans它结合了对数刻度和线性刻度(这与上面的示例图不完全相同,因为它似乎变为 0,但它可能已经足够接近了)。

combine_trans <- function(breakpoint, tr1, tr2,
                          br1 = tr1$breaks,
                          br2 = tr2$breaks) {
    # combine two transformations.
    # can also be given functions to determine tick marks
    trans_breakpoint <- tr1$transform(breakpoint)
    trans <- function(x) {
        # choose which transformation to apply
        ifelse(x < breakpoint,
               tr1$transform(x),
               tr2$transform(x - breakpoint) - trans_breakpoint)
    }
    inv <- function(x) {
        # inverse of both transformations
        ifelse(x < trans_breakpoint,
               tr1$inverse(x),
               breakpoint + tr2$inverse(x + trans_breakpoint))
    }
    br <- function(x) {
        # combine break choices from both scales
        br1 <- br1(c(x[x < breakpoint], breakpoint))
        br2 <- br2(c(breakpoint, max(x[x > breakpoint])))
        br <- c(br1, br2)
        br[br > 0]
    }
    # combine domains
    dom <- c(max(tr1$domain[1], tr2$domain[1]), min(tr1$domain[2], tr2$domain[2]))
    trans_new("combined", trans, inv, breaks = br, domain = dom)
}

# combine log10 transformation and identity transformation
combined <- combine_trans(1, log10_trans(), identity_trans())

ggplot(dat, aes(x, y)) +
    geom_point() +
    scale_y_continuous(trans = combined) +
    geom_hline(yintercept = 1, lty = 2)
Run Code Online (Sandbox Code Playgroud)

在此处输入图片说明

请注意,我手动添加了一条带 的水平虚线geom_hline,这至少有助于引起对不连续性的注意。

您可以在此处看到不连续变换的另一个示例,它在轴上引入了一个简单的中断。

请注意,此代码很复杂,可能需要自定义。更重要的是,它仍然有些误导:一般来说,您可能应该将数据分成单独的图,或者将所有内容都放在对数刻度上。尽管如此,值得知道的是 ggplot2 可以让您访问这些工具,即使您用它们“用脚射击”!