在R语言中获取社交网络传递率组分，排查igraph计算结果不符问题

需求背景 #

处理社交网络数据时，需要分别获取三角形（闭合三元组）数量和连通三元组（三个节点至少由两条边连接的三元组）数量，用于二项式模型统计分析，而非直接使用transitivity函数返回的聚类系数结果。

问题重现 #

尝试用以下代码计算全局传递性，但结果与transitivity(exampleNet, type="global")不一致：

(length(triangles(exampleNet))/3)/graph.motifs.no(exampleNet, size=3)

igraph官方文档说明 #

"global" 无向图的全局传递性（有向图也会被当作无向图处理），即图中三角形数量与连通三元组数量的比值。对于有向图，边的方向会被忽略。

问题原因 #

你对graph.motifs.no的返回值理解有误：

• triangles(exampleNet)返回每个三角形的三个节点ID，每个三角形被重复计数3次，除以3得到正确的三角形数量，这部分是对的。
• graph.motifs.no(exampleNet, size=3)返回的是所有3节点子图的类型统计向量（无向图下共4种类型）：
  1. 0条边的孤立三元组
  2. 1条边的三元组
  3. 2条边的连通三元组（开三元组）
  4. 3条边的三角形（闭合三元组）
• transitivity函数的分母是所有至少包含两条边的三元组数量（即类型3+类型4的和），但你的代码直接用了整个向量（或其求和值），包含了不符合要求的0边、1边三元组，导致结果偏差。

正确的计算方式 #

以下代码可准确获取所需数值，并验证与transitivity函数的一致性：

# 计算三角形（闭合三元组）数量
tri_count <- length(triangles(exampleNet)) / 3

# 统计所有3节点子图的类型数量
motif_3_counts <- graph.motifs.no(exampleNet, size = 3)

# 若需要transitivity分母定义的"连通三元组"（至少两条边的三元组）
connected_triads_total <- motif_3_counts[3] + motif_3_counts[4]

# 若需要"恰好两条边"的开三元组数量
open_triads_count <- motif_3_counts[3]

# 手动计算全局传递性，与内置函数对比
manual_trans <- tri_count / connected_triads_total
all.equal(manual_trans, transitivity(exampleNet, type = "global"))

内容的提问来源于stack exchange，提问作者Pierre Raynerd