【发布时间】:2020-05-26 12:49:55
【问题描述】:
考虑一个将每个变量相互关联的数据框:
iris_cor <- structure(list(x = c("Sepal.Length", "Sepal.Length", "Sepal.Length",
"Sepal.Length", "Sepal.Width", "Sepal.Width", "Sepal.Width",
"Sepal.Width", "Petal.Length", "Petal.Length", "Petal.Length",
"Petal.Length", "Petal.Width", "Petal.Width", "Petal.Width",
"Petal.Width"), y = c("Sepal.Length", "Sepal.Width", "Petal.Length",
"Petal.Width", "Sepal.Length", "Sepal.Width", "Petal.Length",
"Petal.Width", "Sepal.Length", "Sepal.Width", "Petal.Length",
"Petal.Width", "Sepal.Length", "Sepal.Width", "Petal.Length",
"Petal.Width"), r = c(NA, -0.117569784133002, 0.871753775886583,
0.817941126271576, -0.117569784133002, NA, -0.42844010433054,
-0.366125932536439, 0.871753775886583, -0.42844010433054, NA,
0.962865431402796, 0.817941126271576, -0.366125932536439, 0.962865431402796,
NA)), row.names = c(NA, -16L), class = c("tbl_df", "tbl", "data.frame"
), .Names = c("x", "y", "r"))
如果我想对此进行网络图,我会执行以下操作:
library(igraph)
iris_igraph <- graph_from_data_frame(iris_cor)
现在,如果我只对与 Sepal.Length 相关的相关性感兴趣,我不清楚如何有效地提取这些信息。尝试这两种方法都不成功。第一个是因为错误,第二个是因为我只得到Sepal.Length 的第一个实例而不是全部:
> subgraph.edges(iris_igraph, E(iris_igraph)[name == "Sepal.Length"])
Error in eval(x$expr, data, x$env) : object 'name' not found
>
> subgraph.edges(iris_igraph, V(iris_igraph)[name == "Sepal.Length"])
IGRAPH dc7408b DN-- 1 1 --
+ attr: name (v/c), r (e/n)
+ edge from dc7408b (vertex names):
[1] Sepal.Length->Sepal.Length
这会产生所需的输出:
subgraph.edges(iris_igraph, E(iris_igraph)[1:4])
但是,当 igraph 对象的所需子集不按顺序或有很多节点/边时,编号索引不实用。
谁能推荐一种将 igraph 对象子集作为边缘名称的方法?
【问题讨论】: