生物医学网络科学Biomedical Network Science
小世界与无标度网络模型07 / 17

测量优先连接与非线性优先连接

阅读
约 2 分钟

BA 模型把优先连接写成一个假设。它在真实网络中是否成立、函数形式是不是线性,可以用网络的时间序列数据直接检验。

从两个时间快照测量 Π(k)\Pi(k)

如果优先连接成立,节点在一小段时间内获得的新边数应当与它当前的度成正比:Δki/Δt∝Π(ki)\Delta k_i/\Delta t\propto\Pi(k_i)。Jeong、Néda 与 Barabási(2003)提出的做法是:

  1. 取网络在 t0t_0 与 t0+Δtt_0+\Delta t 两个时刻的快照(Δt\Delta t 相对网络规模很短);
  2. 对 t0t_0 时度为 kk 的所有节点,统计它们在 Δt\Delta t 内获得的新边数,平均后得到 Π(k)\Pi(k) 的估计;
  3. 由于单个度类的样本很少、噪声大,改看累积函数 π(k)=∑k′<kΠ(k′).\pi(k)=\sum_{k'<k}\Pi(k'). 若 Π(k)∝kα\Pi(k)\propto k^{\alpha},则 π(k)∝kα+1\pi(k)\propto k^{\alpha+1}:无优先连接(α=0\alpha=0)时斜率为 1,线性优先连接时斜率为 2。

用这种方法,引文、互联网等网络得到的 α\alpha 接近 1;科研合作网络与演员合作网络则是亚线性的(α<1\alpha<1)。本专题代码示例在模拟网络上演示了这一测量。

蛋白质相互作用网络中的优先连接:Eisenberg 与 Levanon(2003, Phys. Rev. Lett.)没有时间快照,于是借助跨物种比较给酵母蛋白质估计“演化年龄”(同源基因在多远的物种中仍能找到)。他们发现:越古老的蛋白,在当前相互作用网络中的连接越多;蛋白在演化中获得的相互作用数与其连接度成正比。这是用比较基因组学替代时间序列、在生物网络中检验优先连接的早期范例。注意年龄估计本身依赖序列检测灵敏度,古老且研究充分的蛋白也更可能测到更多相互作用,这两个偏倚会同方向地放大结果。

非线性优先连接

设 Π(k)∝kα\Pi(k)\propto k^{\alpha}。Krapivsky、Redner 与 Leyvraz(2000)的分析表明,只有 α=1\alpha=1 才给出纯幂律:

  • 亚线性(0<α<10<\alpha<1):富者愈富被削弱,度分布变为拉伸指数型 pk∼k−αexp⁡(−ck1−α)p_k\sim k^{-\alpha}\exp\!\big(-c\,k^{1-\alpha}\big);最大度只按 (ln⁡t)1/(1−α)(\ln t)^{1/(1-\alpha)} 增长,枢纽受到明显限制。
  • 线性(α=1\alpha=1):BA 模型,γ=3\gamma=3,kmax∼t1/2k_{\max}\sim t^{1/2}。
  • 超线性(α>1\alpha>1):富者愈富被放大,出现“赢家通吃”:少数(极端情况下一个)节点获得几乎所有新连接,最大度按 tt 线性增长,网络呈“中心–辐条”(hub-and-spoke)结构;α>2\alpha>2 时这种凝聚从一开始就发生。

所以观察到的度分布偏离纯幂律,可能只是优先连接略偏非线性的结果——例如亚线性的合作网络更适合用拉伸指数描述。

测量 Π(k)\Pi(k) 需要节点“出生”与连边的时间信息。只有单一快照时,常见的替代做法(用年龄代理时间、用横断面的度–年龄关系反推)都会混入适应度差异与观测偏倚:一个节点度高,可能因为它早到、可能因为它“更有吸引力”,也可能只是被研究得更多。用快照数据声称“证明了优先连接”需要非常谨慎。