在 R 中使用 iGraph 获取路径的权重
Obtain weight of a path with iGraph in R
我有一个从数据框创建的图表,采用 from, to, cost
列的形式。我还有一条有效的路径(作为一系列顶点,采用 igraph
的 vpath
格式)(我的图是有向的)。
是否有任何函数可以根据我的图形和路径给出该路径的成本?
我不是在求最短路径,因为我从all_shortest_paths
获得了路径。但是,我只得到了节点的继承,而没有得到我也需要的路径的权重。
编辑:数据
这是我转换成图表的数据框:http://www.sharecsv.com/s/47209742f0052a37e17db37ea3af63ac/arcsWithCost.csv
而我的路径是15 4 50 212 183 112 114 37 228 119
您将路径作为顶点序列。如果您有边序列,这将很容易 - 只需将每条边的权重相加即可。因此,您需要做的是将顶点序列转换为边序列。这就是 get.edge.ids
所做的,尽管您需要将数据转换为正确的格式。
既然你没有提供任何数据,我就随便举个例子说明。
library(igraph)
set.seed(1234)
g = erdos.renyi.game(10,0.15, directed=TRUE)
E(g)$weight = sample(5, length(E(g)), replace=TRUE)
plot(g)
好的,现在假设我们要对从节点 1-4-10-3 开始的路径上的权重求和。
我假设你有一个列表 c(1,4,10,3)
VP = c(1,4,10,3)
EP = rep(VP, each=2)[-1]
EP = EP[-length(EP)]
E(g)$weight[get.edge.ids(g, EP)]
[1] 1 5 4
sum(E(g)$weight[get.edge.ids(g, EP)])
[1] 10
加法:
添加到问题的数据中,有110个节点,但编号最多为281。这些数字是节点的标签,而不是节点ID。您可以使用标签来访问节点,但必须将它们转换为字符串才能将它们视为标签。此代码适用于您的示例。
VP = c(15, 4, 50, 212, 183, 112, 114, 37, 228, 119)
EP = rep(VP, each=2)[-1]
EP = EP[-length(EP)]
E(g)$cost[get.edge.ids(g, as.character(EP))]
sum(E(g)$cost[get.edge.ids(g, as.character(EP))])
我有一个从数据框创建的图表,采用 from, to, cost
列的形式。我还有一条有效的路径(作为一系列顶点,采用 igraph
的 vpath
格式)(我的图是有向的)。
是否有任何函数可以根据我的图形和路径给出该路径的成本?
我不是在求最短路径,因为我从all_shortest_paths
获得了路径。但是,我只得到了节点的继承,而没有得到我也需要的路径的权重。
编辑:数据
这是我转换成图表的数据框:http://www.sharecsv.com/s/47209742f0052a37e17db37ea3af63ac/arcsWithCost.csv
而我的路径是15 4 50 212 183 112 114 37 228 119
您将路径作为顶点序列。如果您有边序列,这将很容易 - 只需将每条边的权重相加即可。因此,您需要做的是将顶点序列转换为边序列。这就是 get.edge.ids
所做的,尽管您需要将数据转换为正确的格式。
既然你没有提供任何数据,我就随便举个例子说明。
library(igraph)
set.seed(1234)
g = erdos.renyi.game(10,0.15, directed=TRUE)
E(g)$weight = sample(5, length(E(g)), replace=TRUE)
plot(g)
好的,现在假设我们要对从节点 1-4-10-3 开始的路径上的权重求和。
我假设你有一个列表 c(1,4,10,3)
VP = c(1,4,10,3)
EP = rep(VP, each=2)[-1]
EP = EP[-length(EP)]
E(g)$weight[get.edge.ids(g, EP)]
[1] 1 5 4
sum(E(g)$weight[get.edge.ids(g, EP)])
[1] 10
加法:
添加到问题的数据中,有110个节点,但编号最多为281。这些数字是节点的标签,而不是节点ID。您可以使用标签来访问节点,但必须将它们转换为字符串才能将它们视为标签。此代码适用于您的示例。
VP = c(15, 4, 50, 212, 183, 112, 114, 37, 228, 119)
EP = rep(VP, each=2)[-1]
EP = EP[-length(EP)]
E(g)$cost[get.edge.ids(g, as.character(EP))]
sum(E(g)$cost[get.edge.ids(g, as.character(EP))])