使用igraph在R中按边属性对邻居进行排序和子集

时间:2015-03-19 20:36:34

标签: r igraph network-analysis

在R中使用igraph,对于某个节点x,我想根据x与该邻居节点之间的边缘属性列出前三个相邻节点

创建定向的,称重的样本图:

set.seed(42)
library(igraph)
n <- 10
adjm <- matrix(sample(0:100, n*10), nc=n)
colnames(adjm) <- rownames(adjm) <- letters[1:n]
g <- graph.adjacency(adjm, weighted=TRUE)

基于输入邻接矩阵的edge属性(此处为权重),x的前三个外围边缘:

x <- 'e'
adjm[x,][order(adjm[x,], decreasing = TRUE)][1:3]

输出:

 i  a  b 
86 62 40 

当前的方法相当麻烦:为邻居选择邻居和边缘,添加到数据帧,对数据帧进行排序并选择前三:

x <- 'e'
tab <- data.frame(cbind(
  name=V(g)[neighbors(g,x, mode='out')]$name,
  weight=E(g)[x %->% neighbors(g,x, mode='out')]$weight)) # or %--%, %<-%
tab <- tab[order(tab$weight, decreasing=TRUE),]
head(tab,3)

输出:

  name weight
8    i     86
1    a     62
3    c      6

有更优雅的方法吗?

1 个答案:

答案 0 :(得分:4)

我不知道这是否优雅,但肯定更短:

e_edges <- E(g)[from(x)]
e_top_weights <- order(e_edges$weight, decreasing=TRUE)[1:3]
E(g)[ as.vector(e_edges)[e_top_weights] ]

边缘序列:

#> [48] e -> i
#> [41] e -> a
#> [42] e -> b

最后一步很麻烦,因为e_edges为它定义了一个奇怪的索引操作符,你需要将它转换为一个向量。这将在igraph的下一个版本中发生变化,您将能够编写自然的

e_edges[e_top_weights]

如果您有任何疑问,请告诉我。