我有一个转换矩阵如下:
1. A A B
2. B C A
3. A C C
其中每列代表句点,每行代表一个代理,每个字母代表一个状态。我想创建一个像Sankey Diagram这样的图,它显示每个时期从状态到状态的转换。
代理人的身份并不重要。我所以想有一个情节
this。
It seems我可以使用networkD3或googleVis软件包。然而,由于每个节点的位置是由包装内生决定的,我不知道如何将时间方面放在X轴上。
任何帮助或替代可视化建议都将受到高度赞赏, 非常感谢,
您可以通过以下方式重现样本数据:
transitiondata <- data.frame("t1"=c("A","B","A"),
"t2"=c("A","C","C"),
"t3"=c("B","A","C"))
答案 0 :(得分:0)
未来的自我回答:ggalluvial
软件包是在此期间开发的,非常适合此任务。为了使用它,我们需要提供整洁的数据。
让我们加载所需的库:
library(ggplot2)
library(ggalluvial)
library(tidyr)
library(dplyr
然后需要先为数据创建标识符,然后才能将其转换为整洁的格式。因此,新数据如下:
transitiondata$id <- c("id1","id2","id3")
转换为整齐的格式
transitiondata_tidy <- transitiondata %>%
gather(time, state, t1,t2,t3) %>%
mutate(time = as.factor(time), state = as.factor(state))
这是我们的数据的样子:
id time state
1 id1 t1 A
2 id2 t1 B
3 id3 t1 A
4 id1 t2 A
5 id2 t2 C
6 id3 t2 C
7 id1 t3 B
8 id2 t3 A
9 id3 t3 C
ggplot2
和ggalluvial
可以达到目的:
ggplot(transitiondata_tidy,
aes(x = time, stratum = state, alluvium = id, fill = state, label = state)) +
geom_stratum() +
geom_text(stat = "stratum", size = 3) +
geom_flow(fill = "darkgrey", color = "black")
我们的过渡(Sankey)图已准备就绪: