如何修复R语言igraph中closeness计算结果的NaN值？

问题根源 #

你的图里节点6和8是完全孤立的——没有任何入边或出边。计算mode="in"的closeness（入近中心性）时，这类节点没有其他节点能到达它的路径，导致无法算出有效数值，返回NaN。

处理方法 #

根据储罐脆弱性分析的业务逻辑，你可以选以下几种常用方式：

1. 直接替换为0 #

如果认为孤立节点的入近中心性为0（比如意味着没有节点依赖它，脆弱性极低），用这段代码替换：

closenessmygraph[is.na(closenessmygraph)] <- 0

2. 贴合现有数据分布替换 #

要是想让孤立节点的取值和已有数据保持一致，可替换成现有非NaN值的最小或最大值：

# 替换为非NaN值的最小值
min_closeness <- min(closenessmygraph, na.rm = TRUE)
closenessmygraph[is.na(closenessmygraph)] <- min_closeness

# 或者替换为最大值
max_closeness <- max(closenessmygraph, na.rm = TRUE)
closenessmygraph[is.na(closenessmygraph)] <- max_closeness

3. 按业务规则自定义值 #

如果你的分析有特定阈值（比如孤立储罐的脆弱性取固定值），直接赋值即可：

closenessmygraph[is.na(closenessmygraph)] <- 0.1  # 这里的0.1可根据实际业务调整

可选预处理方案 #

如果你不想计算时出现NaN，也可以先给孤立节点加自环（注意这会改变图结构，需谨慎）：

# 给节点6、8添加自环，权重可按需设置
mygraph <- add_edges(mygraph, c(6,6,8,8), attr = list(weight = 0))
# 重新计算closeness
closenessmygraph <- closeness(mygraph, mode = "in")

内容的提问来源于stack exchange，提问作者زهرا خدابخش