我有一个包含 2 列的 data.frame:节点 A,节点 B。框架中的每个条目都意味着节点 A 和 B 之间的图中的一条边。
必须有一个很好的单行来将此 data.frame 转换为邻接列表。有什么提示吗?
既然你标记了这个igraph,那么使用内置功能怎么样?
> g <- graph.data.frame( edges )
> adjlist <- get.adjedgelist(g)
唯一需要注意的是顶点的索引为零,这将随着 igraph 0.6 的变化而变化。
又快又脏...
> edges <- data.frame(nodea=c(1,2,4,2,1), nodeb=c(1,2,3,4,5))
> adjlist <- by(edges, edges$nodea, function(x) x$nodeb)
> for (i in as.character(unique(edges$nodea))) {
+ cat(i, ' -> ', adjlist[[i]], '\n')
+ }
1 -> 1 5
2 -> 2 4
4 -> 3
> adjlist
edges$nodea: 1
[1] 1 5
------------------------------------------------------------
edges$nodea: 2
[1] 2 4
------------------------------------------------------------
edges$nodea: 4
[1] 3
> edges <- data.frame(nodea=c(1,2,4,2,1), nodeb=c(1,2,3,4,5))
> attach(edges)
> tapply(nodeb,nodea,unique)
$`1`
[1] 1 5
$`2`
[1] 2 4
$`4`
[1] 3
您甚至如何表示 R 中的邻接表?它需要一组相邻节点的可变大小列表;那么你必须使用 list(); 但是在 R 中有什么好处呢?
我可以想到具有类似 sapply 功能的蹩脚技巧,但它们对每个节点进行线性扫描。但是玩了 1 分钟,这里是:pairlists 的列表,其中每对的第二项是邻接列表。输出比数据结构实际更疯狂。
> edgelist=data.frame(A=c(1,1,2,2,2),B=c(1,2,2,3,4))
> library(plyr)
> llply(1:max(edgelist), function(a) list(node=a, adjacents=as.list(edgelist$B[edgelist$A==a])))
[[1]]
[[1]]$node
[1] 1
[[1]]$adjacents
[[1]]$adjacents[[1]]
[1] 1
[[1]]$adjacents[[2]]
[1] 2
[[2]]
[[2]]$node
[1] 2
[[2]]$adjacents
[[2]]$adjacents[[1]]
[1] 2
[[2]]$adjacents[[2]]
[1] 3
[[2]]$adjacents[[3]]
[1] 4
[[3]]
[[3]]$node
[1] 3
[[3]]$adjacents
list()
[[4]]
[[4]]$node
[1] 4
[[4]]$adjacents
list()