R - ggraph 网络图:指定节点坐标



尝试使用 ggraph 包绘制网络并指定每个节点的坐标。虽然我可以使用igraph包实现这一点 - 但我无法使用ggraph包执行此操作。

# reproducible example to generate a random graph
library(igraph)
g1 <- erdos.renyi.game(20, 1/2)
plot(g1)
# function to produce coordinates for each node in order of the node
# degree (number of links per node)
coord <- function(g){   
    n.nod <- length(V(g))
    mat.c <- matrix(0, nrow = n.nod, ncol = 2)
    deg <- degree(g)
    uniq.deg <- unique(deg)
    min.d <- min(deg); max.d<- max(deg)
    spa <- 10/(max.d - min.d) # how much to increment 0 to 10 space along y-axis
    divi.y.1 <- seq(from = min.d, to=max.d, by = 1)
    divi.y.2 <- seq(from = 0, to=10, by = spa)  # both have same length
    ind.x <- c(); ind.x[1] = 0
    for(x in 2:n.nod){
      ind.x[x] <- ind.x[x-1] + 0.1
      if(ind.x[x] >= 10){
        ind.x[x] = ind.x[x-1] - 0.1
      }  
    } 
 d1 <- data.frame(divi.y.1, divi.y.2)
 # plotting space of grid is (0,0), (0, 10), (10, 10), (10, 0)
 for(i in 1:n.nod){
    # y-axis in order
    inD <-  which(d1$divi.y.1 == deg[i])
    mat.c[i, 2] <- d1[inD,2]
 } 
 mat.c[, 1] <- ind.x
 return(data.frame(mat.c))
}

以下是绘制igraph对象的"老式"方法:

# plot igraph object - the old fashion way
x11()
plot(g1, layout = coord(g1), rescale = T, frame=T,
 vertex.frame.color = "black", edge.color = "lightsteelblue",
 edge.width = 1, vertex.label.cex = 1,
 vertex.label.color = "black", 
 main = "Nodes in order of degree: top have more links, bottom fewer links")

此处提供了ggraph文档的链接。另请参阅 GitHub 存储库以获取ggraph包安装说明(需要>R.3.3 版本(。下面是有效的 ggraph 图(但我没有指定每个节点的坐标(:

library(ggraph)
V(g1)$NaMe <- seq(1:20)
x11()
   ggraph(g1, 'igraph', algorithm = 'kk') +  
   geom_edge_link(colour = "black", alpha = 0.8, show.legend = F) + 
   geom_node_label(aes(label = NaMe)) + ggtitle("ggraph plot: How to allocate coordinate for each node?") + 
   ggforce::theme_no_axes()

这是我尝试为每个节点制作具有指定坐标的ggraph图。按照类似的例子和早期的尝试,ggraph()将节点的坐标传递到图上,我尝试了以下内容:

g <- make_ring(10) + make_full_graph(5)
coords <- layout_(g, as_star())
plot(g, layout = coords)
# try to replicate this example:
coords2 <- add_layout_(g1, coord(g1))

还尝试使用此功能。这很困难,因为文档中没有示例。

Lay <- layout_igraph_manual(g1, coord(g1))
Lay <- layout_igraph_igraph(g1, coord(g1))
x11()
ggraph(g1, 'igraph', algorithm = 'kk' ) + add_layout_(Lay) +
# layout_igraph_circlepack() + 
geom_edge_link(colour = "black", alpha = 0.8, show.legend = F) + 
geom_node_label(aes(label = NaMe)) + ggtitle("ggraph plot: Cannot allocate    coordinate for each node") + 
ggforce::theme_no_axes()

我不相信你直接打电话给layout_igraph_manual是作者的意图。如果将布局命名为"手动",则可以传入手动布局顶点所需的参数。如果我理解正确的话,这让我得到了你所追求的:

coords <- coord(g1)
names(coords) <- c("x","y")
coords <- as.data.frame(coords)
ggraph(g1, layout = "manual", circular = FALSE, node.positions = coords)+
  geom_edge_link(colour = "black", alpha = 0.8, show.legend = F) + 
  geom_node_label(aes(label = NaMe)) + ggtitle("ggraph plot: How to allocate coordinate for each node?") + 
  ggforce::theme_no_axes()

根据文档,手动布局必须是数据帧并包含列xy,这一点很重要。我已经对您的coord()函数的输出执行此操作。您可能只想更新它以返回该形式的数据。

使用 ggraph 时,您通常不必考虑节点坐标 - 它们会根据所选布局自动提供。如果您想事先自己计算布局,只需使用手动布局 - 查看文档以了解layout_igraph_manual

谢谢@ThomasP85(哇!ggraph的开发者(。在本例中,我想指定每个节点的坐标。我已经尝试了layout_igraph_manuallayout_add_layout_和许多其他方法,但似乎不起作用。

这到底是怎么回事

co2<- data.frame(coord(g1))  # my mistake above, the output of coord() should be a matrix
colnames(co2)<- c("x", "y")
Lay<- createLayout(g1, layout = "nicely") # try any layout
Lay$x<- co2$x  # <== now overwrite their coordinates
Lay$y<- co2$y

x11()  # <==== this is ggraph version of the igraph plot above 
ggraph(graph=g1, data=Lay) + 
geom_edge_link(colour = "black", alpha = 0.5, show.legend = F) + 
geom_node_label(aes(label = NaMe)) + ggtitle("Nodes with higher degree at top, lower degree nodes at the bottom") + 
ggforce::theme_no_axes()

@ThomasP85如果您可以使用layout_igraph_manual提供上述代码的工作示例,将不胜感激。

最新更新