带江河绘图包的Sankey地块

Question

附魔。

编辑:解决方案

正如MartineJ和emilliman5所指出的，节点应该有唯一的标签(如下所示)。

library("riverplot")
nodes<-structure(list(ID = c("2011+", "2011-", "2016+", "2016-"), x = c(20, 
20, 30, 30), y = c(50, 40, 50, 40)), class = "data.frame", row.names = c(NA, 
-4L))
edges<-structure(list(N1 = c("2011+", "2011-", "2011+", "2011-"), N2 = 
c("2016+", "2016-", "2016-", "2016+"), Value = c(461, 7, 0, 46)), class = 
"data.frame", row.names = c(NA, -4L))

river <- makeRiver(nodes,edges)
riverplot(river)

我一直在尝试绘制一个Sankey图/江河图(使用河岸包)来描述癌症注册是如何随时间发展的，不过到目前为止，这段代码并没有给我带来什么成功。有人能告诉我这段代码的错误吗？

Warning message: In checkedges(x2$edges, names(x2)) : duplicated edge information, removing 1 edges

这是可疑代码：

library(“riverplot”)

edges<-structure(list(N1 = c("+", "-", "+", "-"), N2 = c("+", "-", "-", "+"), Value = c(664L, 50L, 0L, 46L)), .Names = c("N1", "N2", "Value"), class = "data.frame", row.names = c(NA, -4L))

nodes = data.frame(ID = unique(c(edges$N1, edges$N2)), stringsAsFactors = FALSE)

nodes$x = c(1,2)
rownames(nodes) = nodes$ID

rp <- list(nodes=nodes, edges=edges)

class(rp) <- c(class(rp), "riverplot")

plot(rp)

以及包含在代码中的数据：

N1    N2    Value
 +     +     664
 -     -      50
 +     -       0
 -     +      46

永远感激。

Answer 1

您的节点需要唯一命名，然后使用nodes$labels将其更改为：

library(riverplot)

edges<-structure(list(N1 = c("+", "-", "+", "-"), N2 = c("+", "-", "-", "+"), Value = c(664L, 50L, 0L, 46L)), .Names = c("N1", "N2", "Value"), class = "data.frame", row.names = c(NA, -4L))
edges$N1 <- paste0(edges$N1, "a")
edges$N2 <- paste0(edges$N2, "b")
nodes = data.frame(ID = unique(c(edges$N1, edges$N2)), stringsAsFactors = FALSE)

nodes$x = c(1,1,2,2)
nodes$labels <- as.character(substr(nodes$ID, 1, 1))
rownames(nodes) = nodes$ID

rp <- list(nodes=nodes, edges=edges)

class(rp) <- c(class(rp), "riverplot")

plot(rp)

​

​

Answer 2

看起来，您在N1 (和N2)中多次使用相同的值。试着使它们都不同(每一栏)，然后再试一次，f.i.：

N1: plus1 minus1 plus2 minus2

如果您希望只显示+和-：在makeRiver中，有一个选项**node_labels **