R语言中使用activate子集化igraph图报错，如何保留布局筛选节点？

想要基于顶点属性V(graph)$period筛选出特定节点绘图，同时保留原布局，但运行以下代码时出现错误：

plot.igraph(graph, layout=lo, vertex.label=NA, rescale=T, vertex.size = 4) %>% 
  tidygraph::activate(nodes) %>% 
  filter(period == 1)

错误信息：

Error in UseMethod("activate") : 
no applicable method for 'activate' applied to an object of class "NULL"

解决方案 #

错误原因 #

plot.igraph()的返回值不是可被tidygraph处理的图对象（实际返回的是绘图设备相关内容，甚至可能是NULL），所以不能直接用activate()和filter()对绘图结果进行筛选操作。正确思路是先子集化图对象，再用原布局绘制子图。

方法1：igraph原生实现 #

先筛选出符合条件的顶点生成子图，再提取原布局中对应节点的坐标进行绘图：

# 筛选period为1的顶点，生成子图
sub_graph <- induced_subgraph(graph, vids = V(graph)[period == 1])
# 提取原布局中对应子图节点的坐标
sub_layout <- lo[V(graph)[period == 1], ]
# 绘制子图，使用提取后的布局
plot.igraph(sub_graph, layout = sub_layout, vertex.label = NA, rescale = T, vertex.size = 4)

方法2：tidygraph实现 #

先将igraph对象转为tidygraph进行筛选，再转回igraph，最后用原布局对应坐标绘图：

library(tidygraph)

# 转换为tidygraph并筛选节点
sub_tidy_graph <- graph %>%
  as_tbl_graph() %>%
  activate(nodes) %>%
  filter(period == 1) %>%
  as.igraph()

# 提取原布局中对应节点的坐标
sub_layout <- lo[V(graph)[period == 1], ]
# 绘制子图
plot.igraph(sub_tidy_graph, layout = sub_layout, vertex.label = NA, rescale = T, vertex.size = 4)

关键要点 #

布局lo是原全图的坐标矩阵，子集化后必须提取对应节点的行作为子图的布局，这样才能保证子图节点位置和原图一致。

内容的提问来源于stack exchange，提问作者flâneur