生物医学网络科学Biomedical Network Science
网络结构指标与随机图07 / 19

加权网络的指标

阅读
约 2 分钟

生物医学网络大多带权:功能连接的相关系数、纤维束的条数或各向异性分数、共表达强度、STRING 的置信分数。二值化会丢失信息,并引入阈值选择问题;直接在加权网络上定义指标是另一条路线。

强度与强度–度关系

强度(strength):si=∑jwijs_i=\sum_jw_{ij}。有向网络区分入强度与出强度。

若边权与端点度无关,则 ⟨s|k⟩≈⟨w⟩k\langle s\,|\,k\rangle\approx\langle w\rangle k。Barrat 等(2004)在若干真实网络中观察到 s(k)∼kβs(k)\sim k^{\beta},β>1\beta>1 表示高度节点的边也更“重”。空手道网络中度与强度的相关系数约 0.98,强度排序与度排序的前五名完全相同;但在功能连接这类几乎全连通的加权网络中,度失去意义,强度成为主要的“度”类指标。

加权聚集系数

把三角形的“强弱”计入聚集系数有多种方式,常用两种:

CiO=1ki(ki−1)∑j≠h(ŵijŵihŵjh)1/3,CiB=1si(ki−1)∑j≠hwij+wih2aijaihajh,C_i^{\text{O}}=\frac{1}{k_i(k_i-1)}\sum_{j\ne h}\left(\hat w_{ij}\hat w_{ih}\hat w_{jh}\right)^{1/3},\qquad C_i^{\text{B}}=\frac{1}{s_i(k_i-1)}\sum_{j\ne h}\frac{w_{ij}+w_{ih}}{2}\,a_{ij}a_{ih}a_{jh},

求和遍历 ii 的有序邻居对。COC^{\text{O}} 由 Onnela 等(2005)提出,ŵ=w/max⁡w\hat w=w/\max w,用三角形三条边权重的几何平均衡量“三角形强度”,是 NetworkX 与 BCT(clustering_coef_wu)的默认加权定义;CBC^{\text{B}} 由 Barrat 等(2004)提出,只看 ii 自己的两条边,衡量“ii 的强连接是否落在三角形里”。所有权重相等时,两者都退化为无权 CiC_i。

算例:节点 ii 连向 jj、hh、ll,权重分别为 4、2、2;jj–hh 之间有一条权重为 1 的边,最大权重为 4。无权 Ci=2⋅1/(3⋅2)=0.33C_i=2\cdot1/(3\cdot2)=0.33。Onnela:唯一三角形的强度 (1×0.5×0.25)1/3=0.5(1\times0.5\times0.25)^{1/3}=0.5,有序对计两次,CiO=2×0.5/6≈0.17C_i^{\text{O}}=2\times0.5/6\approx0.17,因为第三条边很弱而低于无权值。Barrat:si=8s_i=8,CiB=2×4+22/(8×2)=0.375C_i^{\text{B}}=2\times\frac{4+2}{2}/(8\times2)=0.375,因为 ii 的最强连接参与了三角形而高于无权值。

两者回答不同的问题,数值可能朝相反方向偏离无权值。Onnela 定义对最大权重敏感:加入一条极端大的边会压低所有节点的 COC^{\text{O}}。

加权路径:先把权重变成长度

最短路径算法把边属性当作长度(越大越远),而多数生物医学网络的权重是强度(越大越近)。必须先转换:ℓij=1/wij\ell_{ij}=1/w_{ij}、ℓij=−log⁡wij\ell_{ij}=-\log w_{ij}(w∈(0,1]w\in(0,1],如概率或置信度,路径长度对应概率乘积的负对数)或 ℓij=1−wij\ell_{ij}=1-w_{ij}。转换方式没有唯一正确答案,结果会随之变化,应作为敏感性分析的一部分。代码示例中,空手道网络节点 0 到 33 的正确最短路径(ℓ=1/w\ell=1/w)经过强度为 5、5、4 的三条边;若误把强度当长度,算法会选择强度为 2、1 的两条最弱的边。

Opsahl 等(2010)用一个参数 α\alpha 在“只看边数”和“只看权重”之间插值:广义度 ki1−αsiαk_i^{1-\alpha}s_i^{\alpha},广义长度 ℓij=wij−α\ell_{ij}=w_{ij}^{-\alpha}。α=0\alpha=0 回到无权指标,α=1\alpha=1 为纯加权指标。

负权重(如偏相关网络中的负边)不能直接用于最短路径(Dijkstra 要求非负长度),也会让强度正负抵消。常见做法是分别分析正、负子网络,或使用符号专用的定义(BCT 的 clustering_coef_wu_sign;症状网络中的预期影响,见症状中心性及其解释)。直接取绝对值会把“抑制”与“促进”混为一谈。