社区检测如何从带协变量的噪声数据中恢复脑网络结构?

贝叶斯模型联合多份噪声网络与节点协变量,更稳健地恢复共享脑网络社区结构。

直接答案

这篇论文的做法是:把同一批节点的多份噪声网络与节点级协变量放在同一个模型里联合建模。其贝叶斯框架将共享的潜在社区结构与按社区块结构定义的假阳性边和假阴性边分开,并在网络信号较弱时让协变量补充信息。模拟研究中,它恢复社区的准确率高于现有方法,当网络噪声处于中等至较高水平,或可用的噪声网络数量很少时,提升最大。

1篇文献引用

本文由 WisPaper 驱动的搜索和论文分析生成。

噪声与协变量通常被分开处理

功能连接并非直接观测得到,而是从神经影像数据中估计出来的,因此同一个人的重复测量可能产生同一网络的多个噪声版本。在群体层面,目标往往是从一组个体网络中恢复一个共享的组织模板,而这些个体网络可以被视为该模板的噪声实现。

节点级协变量提供了第二类信息。它们是每个感兴趣脑区的属性,例如空间、解剖或功能特征。当边稀疏或较弱时,这些属性可以补充网络拓扑。论文的前提是,现有方法只解决了这个问题的一部分:它们通常要么假设网络无噪声,要么忽略协变量信息。

多份噪声副本背后是一个潜在网络

模型把真实网络视为度校正随机块模型(DC-SBM)。节点被分配到各个社区,一条边存在的概率取决于两个节点所属的社区,再加上一个节点特有的度校正项。这样,不同脑区可以在整体连接强度上存在差异,而不会仅因此就被划入不同社区。

随后,每个观测到的噪声网络都通过一个块结构噪声模型与潜在边联系起来。对于两端分别属于社区 k 和 l 的节点对,模型估计一个假阳性率和一个假阴性率。由于这些错误率按社区对定义,噪声可以在社区内部与社区之间不同,而不必被强制为单一的全局水平。

协变量通过一个聚类模型进入:每个节点的协变量来自与其所属社区对应的分布。模型同时处理连续和分类协变量,并假设协变量在社区成分内相互独立。正是这一点,使得解剖或功能属性能够在网络信号较弱或稀疏时提供信息。

社区数量由数据选择而非预先固定

该框架不要求事先指定社区数量 K。后验采样由马尔可夫链蒙特卡洛算法完成,K 则通过最小化广泛适用信息准则(WAIC)来选择。该准则把对噪声网络和对协变量的拟合程度,与一个基于后验方差的惩罚项结合起来。

这对比较方法的读者很重要,因为一个需要预先给定 K 的方法,回答的是另一个实际问题。论文报告了一个开源 R 包 rcdc 来实现该方法,并使用标签重映射和初始化策略来处理标签切换问题并加快收敛。

噪声中等至较高时提升最大

模拟研究在 200 个节点、8 个等大社区上生成真实网络,并设置强、弱两种社区结构。噪声网络按社区内部与社区之间的假阳性率和假阴性率生成,噪声网络数量在 4、10、20 之间变化。协变量包括 8 个连续变量和 2 个分类变量,并设置强、弱两种协变量信号。

聚类准确率用调整兰德指数(ARI)衡量,它比较估计的社区分配与真实分配。论文报告,其方法 RCDC 在所有设置下都优于所有对比方法,包括不同的网络噪声水平、噪声网络数量以及结构和协变量信号强度。摘要指出,当网络噪声处于中等至较高水平,或只有少量噪声网络可用时,提升更大。

比较结果也显示了其他方法在哪些情况下失效。当只有 4 个噪声网络时,随着假阳性率上升,CASCORE 和 SPC 的准确率大幅下降;随着网络数量增加,这一效应变得不那么明显。GBM 虽然对网络噪声建模,但在 4 个噪声网络时表现不佳,在 10 和 20 个时呈现非单调模式,在噪声同质附近准确率最低。作者将此解释为 GBM 需要相对较多的噪声网络和足够强的结构信号。

删边比加边造成的损害更小

模拟把假阳性与假阴性分开考察。当假阴性率上升时,所有方法的表现都相对稳定。作者给出的解释是:假阴性增加会同时删除社区内部和社区之间的边,而真实网络相对稀疏,因此社区结构的强度基本不变。

这是作者对一种模拟模式的解释,而不是关于噪声网络的普遍规律。它受限于研究中使用的稀疏真实网络。这种不对称性对读者仍有参考价值:在这些模拟中,虚假边是更具破坏性的噪声,尤其是在可用噪声网络很少时。

应用识别出可解释且与疾病相关的结构

该方法被应用于阿尔茨海默病神经影像学计划(ADNI)和人类连接组计划(HCP)的功能脑网络。在 ADNI 应用中,论文报告了具有生物学可解释性的结构,并捕捉到与疾病相关的重组。在 HCP 应用中,论文报告了个体层面的变异,与健康年轻成人的群体层面功能组织一致。

HCP 分析纳入了 376 个受试者层面的脑网络,比较了受试者之间的成对调整兰德指数,并与一个参考划分进行比较。论文将这些应用描述为展示了可解释的社区结构,以及与疾病严重程度相关的系统性网络重组模式。这些是应用层面的发现;模拟结果仍是对比竞争方法准确率的直接证据。

证据受所用噪声模型限制

噪声模型假设错误率具有块结构,并在社区对之间对称,同时对假阳性率和假阴性率施加 Beta 先验。讨论部分列出了作者认为有用的扩展方向:加权噪声网络、具有依赖关系或非结构化异质性的噪声,以及用非参数先验或变分算法降低计算成本。

模拟比较也仅限于所测试的设置:200 个节点、8 个社区、强和弱结构信号、强和弱协变量信号,以及 4、10 或 20 个噪声网络。所报告的优势应放在这些条件内理解,而不是当作对任意网络规模或噪声模式的保证。

网络少且噪声大时优先联合建模

论文中最具决策相关性的比较,是只使用一种信息来源的方法与同时使用两种信息的该框架之间的比较。报告的最大提升出现在网络噪声中等至较高,或只有少量噪声网络可用时,而这正是单个噪声网络携带可靠信号很少的情形。

对决定从这项研究中获取什么的读者来说,优先事项是这种组合:在多个实现上显式地对噪声建模,并在边较弱时让节点协变量提供信息。模拟显示,这种组合在所测试的设置中优于替代方法,脑网络应用则显示所得社区可以作生物学解释。边界同样具体:优势是在块结构、对称噪声以及所报告的模拟条件下展示的,而不是在任意噪声依赖或加权网络下展示的。

本文引用的文献

1

Robust Community Detection for Noisy Networks with Covariates: Application to Functional Brain Networks

一个贝叶斯框架联合建模多份噪声网络实现与节点协变量,比竞争方法更准确地恢复共享潜在社区,当网络噪声中等至较高或可用噪声网络很少时提升更大。