igraph包中的子图出错

Question

我想得到图中所有周期的子图。我尝试了下面的代码

 for (i in 1:length(cycles)){
    vids<-as.numeric(unlist(cycles[[i]]))
    subgraph<- induced.subgraph(graph, vids)
 }

但它如下所示抛出错误：

Error in .Call("R_igraph_induced_subgraph", graph, vids - 1, impl, PACKAGE = "igraph") : 
  At iterators.c:759 : Cannot create iterator, invalid vertex id, Invalid vertex id

我发现代码适用于cycles列表中的第二个元素，该元素较短但不是第一个元素。 因此，如果我尝试这将有效，

subgraph<- induced.subgraph(g, c(3,4))

但不是

subgraph<- induced.subgraph(g, c(26, 2, 30, 29, 25, 9, 27, 13, 14, 8, 23, 20, 19, 17, 12, 11, 24, 21, 6, 28, 15,3,4))

此外，欢迎任何替代for循环的建议。

可重现的例子：

    library(igraph)
    graph<-graph(c(1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13,14,14,15,15,16,
           16,17,17,18,18,19,19,20,20,21,21,1,22,23,23,22),directed=T)
    V(graph)$name<-c(26, 2, 30, 29, 25, 9, 27, 13, 14, 8, 23, 20, 19, 17, 12, 11, 
             24, 21, 6, 28, 15,3,4)

    cycles<-list(list(26L, 2L, 30L, 29L, 25L, 9L, 27L, 13L, 14L, 8L, 23L, 
        20L, 19L, 17L, 12L, 11L, 24L, 21L, 6L, 28L, 15L), list(4L, 
        3L))

Answer 1

如果将一个数字作为顶点id传递给induced.subgraph函数，它将按编号访问顶点（在图形编号1到23中），导致“无效顶点ID”错误，原因是像26,29和30这样​​的索引。您希望实际按名称引用顶点，您可以通过传递字符串而不是数字来完成：

for (i in 1:length(cycles)){
  vids <- as.character(unlist(cycles[[i]]))
  subgraph<- induced.subgraph(graph, vids)
  plot(subgraph)
}

现在您已成功提取了两个子图：

没有办法避免循环遍历列表来生成子图，尽管你可以用lapply函数隐藏循环：

subgraphs <- lapply(cycles, function(x) induced.subgraph(graph, as.character(unlist(x))))

现在subgraphs是一个图表列表，每个周期一个。在您的情况下，可以使用subgraphs[[1]]和subgraphs[[2]]访问它们。