如何在R中將基於不同數據集的圖層與ggplot2合並

Question

我想制作一個顯示時間序列的圖，然后顯示一組預測分布，每個分布都以小提琴圖表示。 下面的示例代碼分別創建和繪制一個時間序列（作為線圖）和兩個小提琴圖。

set.seed(12345)
x <- data.frame(time=1:50, dat=rnorm(50))

y1 <- rnorm(500)
y2 <- rnorm(500, sd=5)
y <- data.frame(time=as.factor(c(rep(51,500),rep(52,500))), dat=c(y1,y2)) 

ggplot(x, aes(x=time, y=dat)) +
  geom_line()

ggplot(y, aes(x=time, y=dat)) +
  geom_violin()

如何將它們組合成一個圖表，並具有從時間點1到50（沿x軸）的線圖，然后分別是在時間點51和52的兩個小提琴圖？

Answer 1

我不確定您是否可以在同一軸上繪制離散變量和連續變量。 因此，您必須妥協。 Markus選擇離散化x變量，而我更喜歡使y變量連續。 注意，我已經更改了y的生成方式（刪除了因子）。

library(ggplot2)

set.seed(12345)
x <- data.frame(time=1:50, dat=rnorm(50))

y1 <- rnorm(500)
y2 <- rnorm(500, sd=5)
y <- data.frame(time=c(rep(51, 500), rep(52, 500)), dat=c(y1,y2)) 

ggplot(x, aes(x = time, y = dat)) +
  theme_bw() +
  scale_x_continuous(limits = c(0, 52)) +
  geom_line() + 
  geom_violin(data = y, aes(group = as.factor(time)))

Answer 2

您需要將y$time因子級別轉換為integer ，添加分組變量，然后將data = ...移動到特定幾何。

# Transform your factor variable to its factor levels
y$time <- as.integer(levels(y$time))[y$time]

# Plot with grouping for the violin plot (as x is not a factor anymore) 
ggplot() + 
    geom_line(data = x, aes(x = time, y = dat)) + 
    geom_violin(data = y, aes(x = time, y = dat, group = time))