在R中使用modularity()函数计算传球网络子组模块度报错排查

需求背景 #

有一个包含6名球员（A-F）的传球网络数据Passes，其中ABC属于X队，DEF属于Y队，Number列代表球员间的传球次数，需要计算该网络基于队伍划分的模块度。

原始数据：

Passes <- data.frame(
  id = c("A", "B", "C", "C", "C", "D", "E", "E", "F", "F"), 
  tid = c("B", "C", "A", "D", "B", "F", "A", "B", "D", "E"), 
  Number = c(14,5,1,10,23,4,1,3,5,10)
)

错误分析 #

当前代码存在三个核心问题：

• 顶点属性未同步到图对象：给Passes数据框添加了Team列，但未将该属性赋值给igraph图对象的顶点，导致V(g)$Team不存在，触发"Membership is not a numerical vector"错误。
• membership参数长度不匹配：modularity要求membership向量长度与图的顶点数完全一致（此处为6个顶点），但用Passes$id生成的向量长度为10（边的数量），因此出现"Membership vector size differs from number of vertices"错误。
• 模块度概念误解：模块度是针对整个社区划分方案的指标，不能单独计算某一个队伍的模块度，需先给所有顶点分配对应社区编号，再计算整体模块度。

正确实现代码 #

library(igraph)

# 1. 构建带权重的有向图
g <- graph.data.frame(Passes, directed = TRUE)
E(g)$weight <- Passes$Number  # 直接给边赋值权重

# 2. 给每个顶点绑定队伍属性
team_mapping <- c("A"="X", "B"="X", "C"="X", "D"="Y", "E"="Y", "F"="Y")
V(g)$team <- team_mapping[V(g)$name]

# 3. 生成符合要求的membership向量：X队标记为1，Y队标记为2
membership_vec <- as.integer(factor(V(g)$team, levels = c("X", "Y")))

# 4. 计算有向图的模块度
total_modularity <- modularity(g, membership = membership_vec, directed = TRUE)
cat("整体模块度：", total_modularity, "\n")

补充说明 #

• 模块度的核心逻辑是对比社区内部边的权重占比与随机网络的预期占比，因此不存在单独的"某队伍模块度"概念。如果需要拆分各队伍对整体模块度的贡献，可基于模块度公式手动拆解计算。
• 针对有向图，必须在modularity函数中显式设置directed=TRUE，否则会按无向图规则计算，结果不符合传球网络的实际逻辑。

内容的提问来源于stack exchange，提问作者Emmy W