生物医学网络科学Biomedical Network Science
Methods & Research PracticePresenting results and research frontiers

结果呈现与研究前沿

怎样用图表和报告可信地呈现网络分析结果,以及高阶、时间、多层网络与图表示学习等前沿方向

知识页面
15 页
阅读
约 47 分钟
讲解视频
1 段
代码示例
13 个
结果呈现科学报告图表设计网络可视化同行审阅零模型高阶网络超图单纯复形拓扑数据分析时间网络多层网络相互依赖网络网络嵌入

一项网络分析最终要以图表和文字的形式被他人阅读、检验和使用。本专题的前半部分讨论结果呈现:怎样让每张图回答一个问题并带上对照与不确定性,网络图本身怎样画才不误导,一篇网络分析报告应包含哪些部分、方法部分怎样写才能被复现,口头报告怎样组织,以及怎样审阅他人的网络分析研究。

后半部分把各专题串联起来:先概括贯穿网络科学的五条主线——度与枢纽、模块、网络上的过程、网络是构建出来的、与什么比较;再分页介绍超越“静态、单层、两两相互作用”这一经典网络表示的几个前沿方向:

  • 高阶相互作用:用超图和单纯复形描述三个及以上元素的群体作用,并以拓扑数据分析(持续同调)刻画加权网络的“形状”;
  • 时间网络:边只在特定时刻存在时,传播必须沿时间一致的路径进行;
  • 多层与相互依赖网络:多种关系、多种组学或多个模态放在同一个数学对象里;
  • 网络表示学习与图神经网络:把节点、边或整张图映射为向量,用于预测药物–靶点关系、蛋白质功能和脑网络分类;
  • 可解释性与因果:预测性能好的模型是否告诉我们机制,统计依赖的网络能否回答干预问题。

每个方向都给出定义、可手算的小例子、代表性文献与开放问题,研究前沿与开放问题一页把它们汇总成表。几个方向的共同教训是:表示方法越复杂,“节点和边怎样定义、与什么比较、结论对分析选择是否稳健”这些老问题就越重要。

本专题与网络分析流程与 NetworkX 实践、研究问题、可复现分析与敏感性分析共同构成“计算方法与研究规范”领域:前两者讨论怎样做出可信的分析,本专题讨论怎样把它可信地呈现出来,以及方法还可以向哪里发展。

Contents知识页面

01研究结果的图表呈现基本原则 - 一图一结论:每张图回答一个问题,标题写结论,而不是写“度分布图”。例如“AD 基因的最大连通分支显著大于度匹配随机期望”。 - 总有对照:网络指标单独给出意义不大,要…02网络可视化:原则与常见错误网络图(节点–链接图)是网络研究中最醒目也最容易误导的图。它的价值在于展示局部结构——一条通路、一个模块、一个候选基因的邻居、一组枢纽之间的连接——而不是证明全局性质。全局性质(度…03科学报告的结构与要点报告的组成 一篇网络分析研究报告通常包含以下部分: 部分 要回答的问题 --- --- 摘要 问题、数据、方法、主要发现、意义,各一两句 引言 生物医学背景是什么?为什么要用网络?…04贯穿网络科学的五条主线主线一:度与枢纽。 网络结构指标与随机图中的度分布与小世界与无标度网络模型中的无标度模型 → 网络可控性与稳定性中蓄意攻击枢纽使网络崩解 → 脑连接组中的枢纽脑区与富人俱乐部 → …05高阶相互作用:超图与单纯复形为什么需要高阶表示 普通网络只表示两两关系,而许多生物过程涉及三个以上元素同时作用:多个亚基组成的蛋白复合物、多底物多产物的代谢反应、多个转录因子共同调控一个基因、多个脑区的协同激…06拓扑数据分析与持续同调入门思想 很多网络分析需要先选一个阈值把加权网络变成二值网络,而结论常常随阈值改变(见阈值、加权与稀疏化)。拓扑数据分析(topological data analysis, TDA;…07时间网络从静态网络到时间网络 许多生物医学网络中,边只在特定时刻存在:医院里医护人员与患者的接触、性接触网络、动物个体间的接触、神经元之间的放电序列、症状随日常生活的波动。把所有接触累加成…08多层与相互依赖网络统一框架 当同一组实体之间存在多种关系(PPI、共表达、遗传相互作用),或多组实体之间存在跨类型的关系(基因–蛋白–代谢物–疾病),或网络随时间改变时,可以把它们放在一个多层网络中…09网络表示学习:谱嵌入与随机游走嵌入编码器–解码器框架 网络嵌入(network embedding,图表示学习)把每个节点 映射为一个低维向量 ,使得向量之间的几何关系(内积、距离)反映网络中的关系(邻接、多步可达…10图神经网络消息传递框架 图神经网络(graph neural network, GNN)通过在邻居之间反复“传递消息”学习节点、边或整张图的表示。与浅层嵌入不同,GNN 的参数在所有节点间共…11网络科学中的可解释性与因果两种“解释” 讨论网络模型的可解释性时,要区分两个问题: - 模型解释:模型为什么做出这个预测?例如 GNN 的预测依赖于哪些节点、边和特征。 - 机制解释:生物系统为什么表现出这…12研究前沿与开放问题前沿方向一览 方向 核心问题 代表性文献 开放问题 --- --- --- --- 高阶相互作用 群体相互作用如何改变结构与动力学 Battiston 等 2020;Iacopin…13数学推导与算法一、贝蒂数的秩公式 维链群 是以 -单纯形为基的向量空间(维数 ),边界算子 的矩阵即 。由 ,有 ,同调群 ,其维数 霍奇分解: ,三部分两两正交。由于 是两个半正定矩阵之和, …14生物医学应用复杂网络分析在生物医学工程中的意义,在于它提供了一种介于还原与整体之间的描述层次:既不把系统看成孤立元件的总和,也不停留在“一切都相互关联”的泛泛之谈,而是把“关联的模式”本身作为…15适用范围与常见误用1. 误区:报告越长、图越多越好。→ 正解:每张图回答一个问题;冗余图表会稀释主要结论。 2. 误区:口头报告时把所有分析步骤按时间顺序讲一遍。→ 正解:按研究问题组织,先讲结论,…

VIDEO讲解视频

把结果讲清楚:图表、报告与前沿方向6 分 25 秒
文字稿

大家好,这一期讨论两件事,怎样可信地呈现网络分析,以及方法还能向哪里发展。前半部分讲图表、网络图、报告结构和审阅。后半部分讲高阶相互作用、时间网络、图神经网络和因果这些前沿方向。

每张图回答一个问题,标题写结论,而不是写度分布图。网络指标单独给出意义不大,要和零模型、另一组或另一个阈值放在一起。给出误差线、置信区间、随机基线分布或敏感性曲线。上千个节点的毛线团图不传达信息,只画模块、枢纽或子网络。连续变量用 viridis 这样感知均匀的渐变,避免彩虹色,也不要只靠红绿对比。一张好的结果图至少包含三样东西,观测值、参照和不确定性。

比如这张图,右边把两种策略的曲线画在同一坐标系里对照。蓝线是随机失效,红线是蓄意攻击,差别一眼就能看出来。另一种常用图式,是把零模型分布画成直方图,真实值画成竖线,并标出z分数。度分布最好用对数坐标下的互补累积分布,线性分箱的直方图尾部噪声很大。

全网的节点链接图无法支持存在模块或存在枢纽这样的结论,结论应来自定量分析。只画最强的5%的边,却把图上的结构当作发现,展示用的阈值要在图注中说明。患者组和对照组各自做一次力导向布局,位置不同,读者就无法对照。布局依赖随机初值,要固定种子或保存节点坐标。布局中两个模块靠得近,并不意味着它们之间连接多。稠密网络可以改用按社团排序的邻接矩阵,对角块一目了然。

数据要写出数据库与版本号,比如 STRING 的版本和置信度阈值。节点代表什么,边是物理相互作用、相关、偏相关还是因果,有向还是无向。阈值的类型和取值,以及是否在一个范围内重复分析。零模型的类型和随机网络数目,多重比较的校正方法和效应量。改变阈值、图谱、数据库版本或算法之后,结论是否改变。公开代码、数据位置和软件版本,让别人只看方法部分就能算出同样的数值。

以12分钟的报告为例,开场1分钟内就让听众知道你要证明什么。幻灯片标题写成一句完整的结论,讲图时告诉听众看哪里。在被问到之前,主动说出最重要的一两个局限,以及你怎样检验了它们。审阅别人的研究时,先用自己的话一句话复述它的问题和主要结论。最常见的有效意见不是你的结论错了,而是现有分析还不足以支持这个结论。最值得追问的是,两组网络密度是否相同,零模型是否保留了该保留的性质。

第一条主线是度与枢纽,从无标度模型,到攻击枢纽,再到邻近度的度匹配。第二条是模块,从模块度,到脑功能模块、疾病模块和共表达模块。第三条是网络上的过程,传播、同步、症状激活和网络扩散。第四条是网络是构建出来的,阈值、偏相关和网络推断决定了边的含义。第五条是与什么比较,任何网络指标只有在合适的零模型下才有意义。这五条主线在前沿方向中依然成立,只是换了形式。

普通网络只表示两两关系,一个三元复合物和三对两两结合,投影后是同一个三角形。Iacopini 等人发现,三体感染率足够大时,流行病转变会从连续变成爆发式。时间网络里,A在时刻1接触B,B在时刻2接触C,A能到达C,C却不能到达A。拓扑数据分析不选单一阈值,而是记录环和空腔在全部阈值上的出现与消失。所以把接触聚合成静态网络,通常会高估传播的可达范围。还要注意,用团复形把普通网络升级成单纯复形,并不会带来新的群体相互作用信息。

图神经网络在邻居之间反复传递消息,每多一层,节点就多汇集一步邻域的信息。Xu 等人证明了它的区分能力有上限,比如分不清两个三角形和一个六元环。图卷积本质上是一种平滑,层数多了节点会变得难以区分,所以一般只用2到3层。链接预测按边随机划分时,测试药物的相似药物已在训练集中,性能会被高估。还要与简单方法比较,代码示例中图卷积网络在空手道网络上并不比标签传播更准。Stokes 等人发现抗生素候选分子 halicin,关键在于模型预测之后的湿实验验证。

找出模型依赖的子图,只说明模型学到了什么,它可能来自数据偏倚而不是生物机制。功能连接、偏相关和共表达网络的边,都不能直接回答干预某个节点会怎样。最直接的因果证据来自干预,比如基因敲除、CRISPR 扰动、药物处理和脑刺激。可信的研究会与基线比较,在独立数据中验证,再用干预实验检验关键预测。

小结一下,每张结果图都要同时给出观测值、参照和不确定性。方法部分要写到别人能从原始数据重建同一个网络。方法再复杂,节点和边怎样定义、与什么比较,这些老问题依然存在。建议打开代码示例,试着把零模型分布和真实值画在同一张图里。

INTERACTIVE交互可视化

CODE计算实现

13 个完整可运行的 Python 示例,每个都给出代码、输出图和解读;可在浏览器中直接运行,或下载 Jupyter Notebook · Python 脚本。

References核心文献

Bassett & Sporns 2017:Network neuroscience

  • 背景:网络科学进入神经科学后,需要一个统一的视角来整合从分子到行为的多尺度网络。Nature Neuroscience 20:353–364,DOI 10.1038/nn.4502。
  • 核心内容:综述网络神经科学的核心概念(模块、枢纽、动态、控制)及其跨尺度应用,讨论多层网络和网络控制等方向。
  • 阅读提示:读的时候对照贯穿网络科学的五条主线,标出文中对应的段落;这篇文章几乎覆盖了网络结构指标、网络模型、社团结构、传播、同步、可控性与脑网络等专题的概念。
  • 思考题:文中提到的哪一个方向与你的研究最接近?若把你的研究推进一步,会用到哪个专题的方法?

Battiston 等 2020:Networks beyond pairwise interactions: structure and dynamics

  • 背景:两两相互作用的网络框架能否描述群体相互作用?Physics Reports 874:1–92,DOI 10.1016/j.physrep.2020.05.004。
  • 核心内容:系统介绍超图、单纯复形等高阶表示,高阶网络上的随机游走、传播、同步等动力学,以及在脑、生态、社会系统中的应用。
  • 阅读提示:篇幅很长,建议只读引言和与网络上的传播与免疫、同步与协同动力学相关的传播、同步章节;数学细节可以跳过。对照高阶相互作用:超图与单纯复形中的定义与算例。
  • 思考题:在你的研究对象中,有哪些“三个以上元素同时作用”的现象被两两网络遗漏了?

Holme & Saramäki 2012:Temporal networks

  • 背景:接触、通信与神经活动都发生在特定时刻,静态网络无法表示时间顺序。Physics Reports 519(3):97–125,DOI 10.1016/j.physrep.2012.03.001。
  • 核心内容:时间网络的表示(接触序列、区间图)、时间一致路径与可达性、时间网络上的中心性与基序、随机化参照模型,以及传播等动力学过程。
  • 阅读提示:重点读时间一致路径和随机化参照模型两部分,对照时间网络中的算例与代码。
  • 思考题:你的数据如果有时间戳,聚合成静态网络时丢失了什么?用哪种随机化可以检验这些信息是否重要?

Kivelä 等 2014:Multilayer networks

  • 背景:多重网络、时间网络、相互依赖网络等概念在不同文献中名称和定义各异。Journal of Complex Networks 2(3):203–271,DOI 10.1093/comnet/cnu016。
  • 核心内容:提出以节点–层二元组为基础的统一框架与术语,整理了多层网络的结构指标、动力学与各类特例之间的对应关系,是多组学网络的理论基础。
  • 阅读提示:先读一般定义与特例分类部分,再对照多层与相互依赖网络与多层网络。
  • 思考题:你的多组学或多模态数据中,哪些层共享同一组节点(多重网络),哪些层是不同类型的实体(需要层间依赖边)?

Grover & Leskovec 2016:node2vec: Scalable feature learning for networks

  • 背景:怎样从网络结构中自动学习节点特征,同时兼顾社团(同质性)与结构角色?Proceedings of the 22nd ACM SIGKDD,855–864,DOI 10.1145/2939672.2939754。
  • 核心内容:用参数 pp、qq 控制的二阶有偏随机游走生成节点序列,再用 skip-gram 学习嵌入,在节点分类与链接预测上与 DeepWalk 等方法比较。
  • 阅读提示:对照网络表示学习中的转移概率算例;注意链接预测实验中测试边如何从网络中移除、负样本如何选取。
  • 思考题:在 PPI 网络上,你希望嵌入捕捉的是“同一复合物”(社团)还是“同为枢纽”(结构角色)?qq 应如何选择?

Zitnik、Agrawal 与 Leskovec 2018:Modeling polypharmacy side effects with graph convolutional networks

  • 背景:多药联用的副作用难以通过临床试验穷尽。Bioinformatics 34(13):i457–i466,DOI 10.1093/bioinformatics/bty294。
  • 核心内容:构建蛋白–蛋白、药物–蛋白、药物–药物(按副作用类型区分)的多关系图,用图卷积网络预测药物组合的特定副作用。
  • 阅读提示:对照疾病模块、网络邻近度与药物重定位中的异质网络与网络药理学;注意作者怎样划分训练和测试数据,思考是否存在泄漏风险。
  • 思考题:与网络邻近度方法相比,这类深度学习方法的优势和代价各是什么?

Krzywinski 等 2012:Hive plots—rational approach to visualizing networks

  • 背景:力导向布局的“毛线团”图难以解读、难以比较。Briefings in Bioinformatics 13(5):627–644,DOI 10.1093/bib/bbr069。
  • 核心内容:提出把节点按规则放到放射状轴上的布局方法,使网络图的形状由节点属性决定,可重复、可比较。
  • 阅读提示:对照网络可视化:原则与常见错误,思考“规则布局”与“力导向布局”各自适合回答什么问题。
  • 思考题:如果要比较患者组与对照组的脑网络,你会用什么规则安排节点的位置?

Data & tools相关数据与软件

Further延伸资源

类型 资源 说明
教材 Network Science(Barabási) 免费在线教材,可作为系统学习网络科学的参考教材;Ch. 8 Network Robustness(级联失效)与 Ch. 10 Spreading Phenomena(接触网络上的传播)可作为相互依赖网络与时间网络的背景
教材 Hamilton 2020,Graph Representation Learning(Morgan & Claypool / Springer),doi.org/10.1007/978-3-031-01588-5;作者主页提供预印本 cs.mcgill.ca/~wlh/grl_book 网络嵌入、知识图谱嵌入、GNN 的理论动机(WL 检验)与图生成模型;本专题网络表示学习与 GNN 两页按其覆盖面组织
教材 Bianconi 2018,Multilayer Networks: Structure and Function(Oxford University Press),doi.org/10.1093/oso/9780198753919.001.0001 多层网络的结构指标、渗流与扩散理论
教材 Bianconi 2021,Higher-Order Networks(Cambridge Elements),doi.org/10.1017/9781108770996 单纯复形、霍奇拉普拉斯与高阶动力学的简明导论
综述 Boccaletti 等 2014,doi.org/10.1016/j.physrep.2014.07.001 多层网络的结构与动力学综述
综述 Bick 等 2023,What are higher-order networks?,doi.org/10.1137/21M1414024 何时需要高阶表示、各类高阶表示之间的关系
综述 Torres 等 2021,doi.org/10.1137/20M1355896 复杂系统的表示选择:两两网络、超图、时间网络、多层网络各适合什么问题
综述 Holme 2015,doi.org/10.1140/epjb/e2015-60657-4 时间网络理论的较新进展
综述 Carlsson 2009,Topology and data,doi.org/10.1090/S0273-0979-09-01249-X 拓扑数据分析的经典导论
综述 Zitnik 等 2019,doi.org/10.1016/j.inffus.2018.09.012 生物医学多源数据整合的机器学习原则、实践与机会
综述 Li、Huang 与 Zitnik 2022,doi.org/10.1038/s41551-022-00942-x 图表示学习在生物医学与医疗中的应用
综述 Gaudelet 等 2021,doi.org/10.1093/bib/bbab159 图机器学习在药物发现与开发中的应用
综述 Gehlenborg 等 2010,doi.org/10.1038/nmeth.1436;O’Donoghue 等 2018,doi.org/10.1146/annurev-biodatasci-080917-013424 组学与生物医学数据可视化的工具与原则
论文 Merico 等 2009,doi.org/10.1038/nbt.1567 用网络解读生物数据时的布局与视觉映射要点
论文 Crameri 等 2020,doi.org/10.1038/s41467-020-19160-7;Wong 2011,doi.org/10.1038/nmeth.1618 科学图表中颜色的误用与色觉障碍友好配色
论文 Weissgerber 等 2015,doi.org/10.1371/journal.pbio.1002128 用展示原始数据分布的图代替均值柱状图
论文 Barabási, Gulbahce & Loscalzo 2011,doi.org/10.1038/nrg2918 网络医学综述,对应分子网络与网络医学专题;论文精读见 Barabási 2011
论文 Broido & Clauset 2019,doi.org/10.1038/s41467-019-08746-5 “无标度网络是否普遍”之争,对应小世界与无标度网络模型,体会“与什么比较”;论文精读见 Broido 2019
论文 Buldyrev 等 2010,doi.org/10.1038/nature08932 相互依赖网络的级联失效;论文精读见 Buldyrev 2010
论文 Sandve 等 2013,doi.org/10.1371/journal.pcbi.1003285 提交报告前对照的可复算十条规则
工具 NetworkX 文档 本站代码示例使用的 Python 库,网络分析中最常查阅

本专题正文引用的其他文献:Belkin 与 Niyogi 2003(10.1162/089976603321780317);Perozzi 等 2014,DeepWalk(10.1145/2623330.2623732);Qiu 等 2018(10.1145/3159652.3159706);Yue 等 2020(10.1093/bioinformatics/btz718);Kipf 与 Welling 2017(ICLR,arXiv:1609.02907);Hamilton、Ying 与 Leskovec 2017,GraphSAGE(NeurIPS);Veličković 等 2018,GAT(ICLR);Xu 等 2019,GIN(ICLR);Gilmer 等 2017,MPNN(ICML);Ying 等 2019,GNNExplainer(NeurIPS);Li、Han 与 Wu 2018(10.1609/aaai.v32i1.11604);Stokes 等 2020(10.1016/j.cell.2020.01.021);Gligorijević 等 2021,DeepFRI(10.1038/s41467-021-23303-9);Jumper 等 2021,AlphaFold(10.1038/s41586-021-03819-2);Kawahara 等 2017,BrainNetCNN(10.1016/j.neuroimage.2016.09.046);Li 等 2021,BrainGNN(10.1016/j.media.2021.102233);Rudin 2019(10.1038/s42256-019-0048-x);Schölkopf 等 2021(10.1109/JPROC.2021.3058954);Klamt 等 2009(10.1371/journal.pcbi.1000385);Iacopini 等 2019(10.1038/s41467-019-10431-6);Skardal 与 Arenas 2019(10.1103/PhysRevLett.122.248301);Benson、Gleich 与 Leskovec 2016(10.1126/science.aad9029);Giusti 等 2015(10.1073/pnas.1506407112);Petri 等 2014(10.1098/rsif.2014.0873);Sizemore 等 2018(10.1007/s10827-017-0672-6);Goh 与 Barabási 2008(10.1209/0295-5075/81/48002);Scholtes 等 2014(10.1038/ncomms6024);Lambiotte、Rosvall 与 Scholtes 2019(10.1038/s41567-019-0459-y);Vanhems 等 2013(10.1371/journal.pone.0073970);De Domenico 等 2013(10.1103/PhysRevX.3.041022);Gómez 等 2013(10.1103/PhysRevLett.110.028701);Battiston、Nicosia 与 Latora 2014(10.1103/PhysRevE.89.032804);Mucha 等 2010(10.1126/science.1184819);Halu 等 2019(10.1038/s41540-019-0092-5);Valdeolivas 等 2019(10.1093/bioinformatics/bty637);De Domenico 2017(10.1093/gigascience/gix004);Pratapa 等 2020(10.1038/s41592-019-0690-6);Gysi 等 2021(10.1073/pnas.2025581118);Lex 等 2014,UpSet(10.1109/TVCG.2014.2346248);Holten 2006(10.1109/TVCG.2006.147);Behrisch 等 2016(10.1111/cgf.12935);Krzywinski 等 2009,Circos(10.1101/gr.092759.109)。

Open questions开放问题

  1. 用五条主线中的任意两条,为你正在进行的一项网络研究画一张“知识来源图”:每一步分别用到了哪些专题的概念?哪一步还需要补充其他领域的方法?
  2. 在高阶网络、图神经网络、多组学网络三个方向中,选一个与你的研究最相关的,设想一个具体的生物医学问题,并说明它为什么不能用普通的两两网络很好地回答。
  3. 找一张你或他人论文中的网络图,按网络可视化:原则与常见错误中的清单逐条检查:它支持的结论是什么?换成矩阵图或统计图是否更清楚?
  4. 对一个带时间戳的生物医学数据集(如医院接触、电子病历中的诊断序列),聚合网络与时间网络的结论在什么条件下会不同?你会用哪种随机化来检验?
  5. 设计一个药物–靶点链接预测的评估方案:训练/测试怎样划分才能模拟“为一个全新药物预测靶点”的真实场景?应该与哪些简单基线比较?
  6. 统计依赖网络中的“枢纽”在什么条件下可以被解释为干预靶点?需要补充哪些数据或实验?