0

我有使用该ggalluvial软件包制作冲积地块的经验。但是，我遇到了一个问题，我试图创建一个具有两个不同来源的冲积图，这些来源收敛到 1 个变量。

这是示例数据

library(dplyr)
library(ggplot2)
library(ggalluvial)

data <- data.frame(
  unique_alluvium_entires = seq(1:10),
  label_1 = c("A", "B", "C", "D", "E", rep(NA, 5)),
  label_2 = c(rep(NA, 5), "F", "G", "H", "I", "J"),
  shared_label = c("a", "b", "c", "c", "c", "c", "c", "a", "a", "b")
)

这是我用来制作情节的代码

#prep the data
data <- data %>%
  group_by(shared_label) %>%
  mutate(freq = n())

data <- reshape2::melt(data, id.vars = c("unique_alluvium_entires", "freq"))
data$variable <- factor(data$variable, levels = c("label_1", "shared_label", "label_2"))

#ggplot
ggplot(data,
       aes(x = variable, stratum = value, alluvium = unique_alluvium_entires,
           y = freq, fill = value, label = value)) +
  scale_x_discrete(expand = c(.1, .1)) + 
  geom_flow() +
  geom_stratum(color = "grey", width = 1/4, na.rm = TRUE) +
  geom_text(stat = "stratum", size = 4) +
  theme_void() +
  theme(
   axis.text.x = element_text(size = 12, face = "bold")
  )

结果图（显然我还不能嵌入图像）

如您所见，我可以删除这些NA值，但shared_label不能正确“堆叠”。每个唯一的行应该在列中堆叠在一起shared_label。这也将解决大小问题，使它们沿 y 轴大小相等。

任何想法如何解决这一问题？我已经尝试过ggsankey，但出现了同样的问题，我无法删除NA值。非常感谢任何提示！

4

1 回答 1

0

该图是“流”统计变换的预期结果，是“流”图形对象的默认值。（即geom_flow()= geom_flow(stat = "flow")。）看起来您想要的是指定“冲积层”统计转换。下面我使用了你所有的代码，但只复制和编辑了ggplot()调用。

#ggplot
ggplot(data,
       aes(x = variable, stratum = value, alluvium = unique_alluvium_entires,
           y = freq, fill = value, label = value)) +
  scale_x_discrete(expand = c(.1, .1)) +
  geom_flow(stat = "alluvium") +  # <-- specify alternate stat
  geom_stratum(color = "grey", width = 1/4, na.rm = TRUE) +
  geom_text(stat = "stratum", size = 4) +
  theme_void() +
  theme(
    axis.text.x = element_text(size = 12, face = "bold")
  )
#> Warning: Removed 2 rows containing missing values (geom_text).

由reprex 包于 2021-12-10 创建(v2.0.1)

于 2021-12-10T17:58:11.830 回答