在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
有更优雅的方法吗?
答案 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]
如果您有任何疑问,请告诉我。