R上带有facet_zoom的ggforce ggplot2误差

Question

我在R4.0.2中有一个data.frame，在一列中有一个连续变量，在另一列中有两个可能值的分类变量(变量'type'：已知或新颖)，我用它来对它们进行不同的颜色(使用一个来自ggsci 2.9包的调色板)。我用stat_bin 3.3.2表示直方图( ggplot2 3.3.2)，我希望使用ggforce 0.3.2的facet_zoom函数只缩放属于其中一种“类型”的数据(使用选项zoom.data，就像在http://cran.univ-paris1.fr/web/packages/ggforce/vignettes/Visual_Guide.html#contextual-zoom上的火山示例中所做的那样)，但是我得到了这个错误：

Error: Aesthetics must be either length 1 or the same as the data (2000): x

可复制的例子：

library(ggplot2)
library(ggsci)
library(ggforce)

testdata <- as.data.frame(sort(rnorm(1000)))
testdata$type <- "known"
testdata[501:1000,2] <- "novel"

# Working code
ggplot(testdata) +
  stat_bin(aes(x=testdata[,1], fill = type), binwidth = 1, color="white") +
  scale_fill_npg() + theme_light() +
  facet_zoom(xlim = c(0, 4), ylim = c(0, 300), horizontal = TRUE, zoom.size = 0.3)

# Desired code
ggplot(testdata) +
  stat_bin(aes(x=testdata[,1], fill = type), data = cbind(testdata, zoom = FALSE), binwidth = 1, color="white") +
  stat_bin(aes(x=testdata[testdata$type == "novel",1]), data = cbind(testdata, zoom = TRUE), binwidth = 0.5) +
  scale_fill_npg() + theme_light() +
  facet_zoom(xlim = c(0, 4), ylim = c(0, 300), horizontal = TRUE, zoom.size = 0.3, zoom.data = zoom)

谢谢!

Answer 1

问题是在第二个stat_bin中将整个数据集作为数据传递。只需传递子设置的df，而不是在aes()中尝试子集。

顺便说一下:我还将数据中的第一个变量重命名为x。

library(ggplot2)
library(ggsci)
library(ggforce)

set.seed(42)

testdata <- data.frame(x = sort(rnorm(1000)))
testdata$type <- "known"
testdata[501:1000,2] <- "novel"

# Desired code
ggplot(testdata) +
  stat_bin(aes(x = x, fill = type), data = cbind(testdata, zoom = FALSE), binwidth = 1, color="white") +
  stat_bin(aes(x = x), data = cbind(testdata[testdata$type == "novel", ], zoom = TRUE), binwidth = 0.5) +
  scale_fill_npg() + theme_light() +
  facet_zoom(xlim = c(0, 4), ylim = c(0, 300), horizontal = TRUE, zoom.size = 0.3, zoom.data = zoom)

​

​

Answer 2

若要仅显示放大图中的type == "novel"数据，请尝试如下：

library(tidyverse)
library(ggsci)
library(ggforce)

testdata <- data.frame(values = sort(rnorm(1000)))
testdata$type <- "known"
testdata[501:1000,2] <- "novel"

# Desired code
ggplot(testdata) +
  stat_bin(aes(x = values, fill = type),
           binwidth = 1, color="white") +
  scale_fill_npg() + theme_light() +
  facet_zoom(zoom.data = ifelse(type == "novel", NA, FALSE), 
             xlim = c(0, 4), ylim = c(0, 300),
             horizontal = TRUE)

​

​