使用 ggalluvial 创建冲积图时出现问题

Question

我有这个数据集：

type <- c(1, 2, NA, 1, 2, NA)
freq <- c(75, 12, 25, 69, 22, 32)
time <- c("before", "before", "before", "after", "after", "after")

df <- data.frame(type , freq, time )

我需要制作一个像这样的图表（显然有不同的数据）：

我尝试按照此处的指南进行操作： https://cran.r-project.org/web/ packages/ggalluvial/vignettes/ggalluvial.html

使用此代码：

ggplot(modechoice,
       aes(x = time, stratum = type, alluvium = time,
           y = freq,
           fill = type, label = type)) +
  scale_x_discrete(expand = c(.1, .1)) +
  geom_flow() +
  geom_stratum(alpha = .5) +
  geom_text(stat = "stratum", size = 3) +
  theme(legend.position = "none") +
  ggtitle("x")

但我收到错误，我的数据未被识别为冲积层。我做错了什么？

Answer 1

您应该有第四个变量，它与您的 time 变量链接，这就是我创建一个额外变量作为示例的原因。您可以使用以下代码：

type <- c("1", "2", "NA", "1", "2", "NA")
freq <- c(75, 12, 25, 69, 22, 32)
time1 <- c("before", "before", "before", "after", "after", "after")
time2 <- c("after", "before", "after", "after", "before", "after") 


df <- data.frame(type , freq, time1, time2)

library(tidyverse)
library(ggalluvial)

ggplot(df,
       aes(y = freq, axis1 = time1, axis2 = time2)) +
  geom_alluvium(aes(fill = type), width = 1/12) +
  geom_stratum(width = 1/12, fill = "black", color = "grey") +
  geom_label(stat = "stratum", aes(label = after_stat(stratum))) +
  scale_x_discrete(limits = c("time1", "time2"), expand = c(.05, .05)) +
  scale_fill_brewer(type = "qual", palette = "Set1") +
  ggtitle("x")

输出：

在 Victor Nielson 评论后编辑：

    type <- c("1", "2", "NA", "1", "2", "NA")
freq <- c(75, 12, 25, 69, 22, 32)
time1 <- c("before", "before", "before", "before", "before", "before")
time2 <- c("after", "after", "after", "after", "after", "after") 


df <- data.frame(type , freq, time1, time2)

library(tidyverse)
library(ggalluvial)

ggplot(df,
       aes(y = freq, axis1 = time1, axis2 = time2)) +
  geom_alluvium(aes(fill = type), width = 1/12) +
  geom_stratum(width = 1/12, fill = "black", color = "grey") +
  geom_label(stat = "stratum", aes(label = after_stat(stratum))) +
  scale_x_discrete(limits = c("time1", "time2"), expand = c(.05, .05)) +
  scale_fill_brewer(type = "qual", palette = "Set1") +
  ggtitle("x")

输出：