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

任意度分布网络的生成:配置模型与隐变量模型

阅读
约 4 分钟

很多时候我们并不关心度分布“从哪里来”,只想问:如果网络除了度分布之外完全随机,会是什么样子?这需要能按指定度分布或度序列生成网络的“静态”模型。它们既是研究度异质性效应的理论工具,也是检验真实网络其他结构(聚类、社团、度相关、基序)是否“超出度序列所能解释”的零模型。

配置模型

配置模型(configuration model):给定度序列 {k1,…,kN}\{k_1,\dots,k_N\}(∑iki=2L\sum_i k_i=2L 为偶数),给节点 ii 装上 kik_i 个“半边”(stub),把全部 2L2L 个半边随机两两配对,每一对构成一条边。所有配对方式等概率。

一个可手算的小例子:度序列 (2,1,1)(2,1,1) 有 4 个半边 a1,a2a_1,a_2(节点 1)、bb(节点 2)、cc(节点 3)。4 个半边共有 3 种配对方式:{a1a2,bc}\{a_1a_2, bc\} 产生节点 1 的自环加边 2–3;{a1b,a2c}\{a_1b, a_2c\} 与 {a1c,a2b}\{a_1c, a_2b\} 都产生路径 2–1–3。所以按配置模型抽样,有 1/3 的概率得到“带自环”的图——小网络、度大的节点特别容易出现自环与重边。

配置模型的主要性质(稀疏、NN 大):

  • 连边概率:节点 ii、jj 之间的期望边数为 kikj2L−1≈kikj2L\dfrac{k_ik_j}{2L-1}\approx\dfrac{k_ik_j}{2L}。这正是模块度定义中的期望项(见社团结构与社区发现)。
  • 余度与朋友悖论:沿一条随机边到达的节点,其度为 kk 的概率正比于 kpkkp_k;去掉来时那条边后的“余度”(excess degree)分布为 qk=(k+1)pk+1⟨k⟩q_k=\dfrac{(k+1)p_{k+1}}{\langle k\rangle}。随机邻居的平均度是 ⟨k2⟩/⟨k⟩\langle k^2\rangle/\langle k\rangle,总不小于 ⟨k⟩\langle k\rangle——“你的朋友平均比你朋友多”。在 2<γ<32<\gamma<3 的网络中这个比值随规模发散。
  • 巨连通分支:Molloy 与 Reed(1995)证明,当 ⟨k2⟩−2⟨k⟩>0⇔⟨k2⟩⟨k⟩>2\langle k^2\rangle-2\langle k\rangle>0 \quad\Longleftrightarrow\quad \frac{\langle k^2\rangle}{\langle k\rangle}>2 时出现巨连通分支(ER 网络中即 ⟨k⟩>1\langle k\rangle>1)。直观上:顺着一条边走到的节点,平均还能再引出 ⟨k2⟩/⟨k⟩−1\langle k^2\rangle/\langle k\rangle-1 条新边,大于 1 就能“繁衍”下去。⟨k2⟩\langle k^2\rangle 发散的无标度网络几乎总满足这个条件,这也是它们对随机删点稳健的原因(见网络鲁棒性、渗流与级联失效)。
  • 路径长度:第一层邻居平均 z1=⟨k⟩z_1=\langle k\rangle 个,第二层约 z2=⟨k2⟩−⟨k⟩z_2=\langle k^2\rangle-\langle k\rangle 个,平均距离近似为 ℓ≈1+ln⁡(N/z1)ln⁡(z2/z1)\ell\approx 1+\dfrac{\ln(N/z_1)}{\ln(z_2/z_1)}(Newman、Strogatz 与 Watts 2001)。
  • 聚集系数:C=(⟨k2⟩−⟨k⟩)2⟨k⟩3NC=\dfrac{(\langle k^2\rangle-\langle k\rangle)^2}{\langle k\rangle^3N},随 NN 趋于 0;但当 ⟨k2⟩\langle k^2\rangle 随 NN 发散时下降得很慢,度异质性本身就会“制造”一部分三角形。这就是为什么判断小世界性时要用保度随机化而不是 ER 网络作参照。

自环、重边与结构截断

配置模型生成的是多重图。度有限方差时,自环与重边的期望数是 O(1)O(1),可以忽略或删除;但在 γ<3\gamma<3 的重尾序列中,枢纽之间的期望边数 kikj/2Lk_ik_j/2L 可能大于 1,重边大量出现。删除它们会改变度序列,并人为制造“枢纽彼此回避”的负度相关。Boguñá、Pastor-Satorras 与 Vespignani(2004)指出,不想产生这种相关的简单图,最大度应不超过结构截断 ks∼⟨k⟩Nk_s\sim\sqrt{\langle k\rangle N};而 2<γ<32<\gamma<3 时自然截断 N1/(γ−1)N^{1/(\gamma-1)} 大于 N\sqrt N,两者冲突。因此,在重尾网络中“只由度序列决定、没有度相关”的简单图在数学上本来就难以实现——实际网络中观察到的枢纽间连接偏少,一部分可能只是这种结构约束(见度相关性与网络结构的高阶统计)。

此外,并不是每个整数序列都能实现为简单图:判断一个度序列是否“可图”(graphical)可用 Erdős–Gallai 条件。

保度随机化

如果手头已经有一个真实网络,最常用的做法不是从头配对,而是反复做双边交换:随机选两条边 (a,b)(a,b)、(c,d)(c,d),改连为 (a,d)(a,d)、(c,b)(c,b),在不产生自环和重边时接受。每次交换都保持所有节点的度不变,交换足够多次后得到“与原网络度序列相同、其余随机”的简单图(算法见数学推导与算法,零模型选择的系统讨论见网络构建与统计验证)。

保度零模型揭示的蛋白质网络结构:Maslov 与 Sneppen(2002, Science)把酵母蛋白质相互作用网络与转录调控网络分别和保度随机网络比较,发现度高的蛋白之间直接相连的比例低于随机预期,枢纽倾向于连到低度蛋白。正因为比较对象保持了完全相同的度序列,这一差异不能归因于“枢纽多”本身。上一小节的结构截断提醒我们:解读这类结果时,还要区分生物学上的“枢纽回避”与简单图约束造成的部分。

隐变量模型

配置模型固定每个节点的确切度。另一类静态模型只固定期望度:给每个节点一个隐变量(权重、适应度)ηi\eta_i,节点对 (i,j)(i,j) 独立地以概率 pij=min⁡(1,ηiηjN⟨η⟩)p_{ij}=\min\!\left(1,\ \frac{\eta_i\eta_j}{N\langle\eta\rangle}\right) 连边。这样节点 ii 的期望度就是 ηi\eta_i,给定 ηi\eta_i 时实际度近似服从均值为 ηi\eta_i 的泊松分布。只要隐变量分布 ρ(η)\rho(\eta) 有幂律尾 ρ(η)∼η−γ\rho(\eta)\sim\eta^{-\gamma},度分布在大 kk 处就有同样的尾部 pk∼k−γp_k\sim k^{-\gamma}。这一类模型有多个名字:Chung–Lu 模型(期望度模型)、Caldarelli 等(2002)的“内禀适应度”模型、Goh、Kahng 与 Kim(2001)的“静态模型”(取 ηi∝i−α\eta_i\propto i^{-\alpha},得到 γ=1+1/α\gamma=1+1/\alpha)。

隐变量模型的优点是各边独立,期望值可以解析计算,也天然没有重边;缺点与配置模型相同——当 ηiηj/(N⟨η⟩)>1\eta_i\eta_j/(N\langle\eta\rangle)>1 需要截断时,就会引入度相关。它的思想在网络推断中也很常见:“每个节点有一个潜在的活跃度或可结合性,连边概率由两端的潜变量决定”。

静态模型只“复制”度分布,不解释度分布的来源。用配置模型重现了某个网络的某项性质,只能说明这项性质可以由度序列解释,不能说明这个网络是“随机配对”生成的。反过来,真实网络与配置模型的差异(更多三角形、社团、度相关)正是值得研究的“非随机结构”。