生物医学网络科学Biomedical Network Science
脑连接组:结构与功能连接07 / 20

阈值与加权

阅读
约 2 分钟

得到 N×NN\times N 的相关矩阵后,很多图论指标要求一个二值网络,于是需要阈值化。这一步看起来是技术细节,实际上对结果影响极大:同一数据,不同阈值,得到的是不同的网络。

绝对阈值:保留 |rij|>τ|r_{ij}| > \tau 的边,τ\tau 固定。不同被试的网络密度会不同。

比例阈值(密度阈值):每个被试保留最强的前 p%p\% 的边,使所有人网络密度相同。

网络密度:D=2EN(N−1)D = \dfrac{2E}{N(N-1)},即实际边数与可能边数之比。

为什么要关心密度?因为几乎所有图论指标都随密度变化:密度越高,聚集系数越高、路径越短、效率越高。如果患者组整体 FC 偏弱,用绝对阈值时患者网络更稀疏,于是“效率降低”可能只是密度降低的副产品(van Wijk 等, 2010, PLoS ONE)。比例阈值消除了密度差异,但也有代价:当某组整体 FC 偏弱时,为了凑够同样的边数,会把更多接近噪声的弱连接纳入网络,反过来让该组网络看起来更“随机”(van den Heuvel 等, 2017, NeuroImage)。两种阈值偏差的量级估计见数学推导与算法。

实际研究中常见的折中方案:

  • 扫描一段密度范围(如 5%–40%),在每个密度上计算指标,报告曲线或曲线下面积(AUC),而不是只挑一个阈值;
  • 直接用加权网络,计算加权聚集系数、加权效率等,避免阈值(见加权脑网络指标与归一化);
  • 保证连通性:先保留最小生成树(MST),再按强度加边;
  • 统计阈值:只保留在组水平显著不为零的边(需做多重比较校正),或对结构网络使用组一致性阈值;
  • 负相关的处理:多数图论指标不直接支持负权,可以取绝对值、只看正相关、或分别分析正负网络,但必须说明选择。

阈值选择不能事后挑。 如果尝试了 10 个阈值、只报告“显著”的那个,就是典型的研究者自由度问题。阈值方案应在分析前确定,或者报告全部结果。本专题的“阈值”交互演示可以直观看到:同一个相关矩阵,阈值从 0.2 提到 0.5,枢纽脑区的排序会明显改变。