开启掘金成长之旅！这是我参与「掘金日新计划 · 12 月更文挑战」的第4天

本文首发于CSDN。

诸神缄默不语-个人CSDN博文目录 cs224w（图机器学习）2021冬季课程学习笔记集合

@[toc]

本章主要内容：本章首先介绍了图中motif / subgraph的概念，以及对motif significance的定义（即定义图中的subgraph要比null model多/少出多少才算显著，以及如何生成null model）。接下来介绍了神经网络模型下的subgraph matching方法（同时也是subgraph的表示方法）。最后介绍如何找到图中出现频率较高的motif / subgraph。

1. Identifying and Counting Motifs in Networks

subgraph subgraph是网络的组成部分，可用于识别和区分不同的网络（可以说是不同种类网络会具有不同特征的subgraph）。使用传统的discrete type matching¹ 方法代价很大，本文会介绍使用神经网络解决subgraph matching问题的方法。
以下图分子式为例：含有羧基（subgraph）的分子（graph）是羧酸（group）

2. Subgraph and Motifs

2.1 Defining Subgraphs and Motifs

subgraph定义对于图 $\mathbf{G}=(\mathbf{V},\mathbf{E})$ 有两种定义其subgraph的方式：
1. node-induced subgraph / induced subgraph：常用定义图中的一个节点子集+原图中两个节点都在该节点子集内的边（即edges induced by the nodes）
  $G'=(V',E')$ 是 node induced subgraph，当且仅当： $V'\subseteq V$ $E'=\{(u,v)\in E\ |\ u,v\in V'\}$ 时
  $G'$ is the subgraph of G induced by $V'$
2. edge-induced subgraph / non-induced subgraph / subgraph 图中的一个边子集+该子集的对应节点
  $G'=(V',E')$ 是 edge induced subgraph，当且仅当： $E'\subseteq E$ $V'=\{v\in V\ |\ (v,u)\in E'\ \text{for some}\ u\}$ 时（我没看懂，意思应该是包含这些边的全部节点。就是说选了边之后节点自然就确定了，所以是non-induced。参考²）
3. 具体使用哪种定义取决于问题领域。如在化学领域中常使用node-induced概念（官能团），在知识图谱中常用edge-induced概念（我们关心的是代表逻辑关系的边）。
前文对subgraph的定义都需要 $V'\subseteq V$ 且 $E'\subseteq E$ ，即 $V'$ 和 $E'$ 都出自原图。如果节点和边出自不同的图但仍有对应关系，如下图所示，我们称这种情况为 $\mathbf{G}_1$ is contained in $\mathbf{G}_2$
graph isomorphism图同构如果 $G_1$ 和 $G_2$ 存在双射关系 $f:V_1\rightarrow V_2$ ，使得当且仅当 $\big(f(u),f(v)\big)\in E_2$ 时， $(u,v)\in E_1$ （即 $G_1$ 中的节点能一一映射到 $G_2$ 中的节点，使节点之间对应的边关系也能同时映射到另一个图所对应的节点之间）。我们称两个图同构。图中ab是uv写错了。如下图左图所示（节点颜色表示映射关系）。因为节点没有固定顺序，所以我们不知道节点之间是怎么映射的，所以我们需要遍历所有可能。检验图是否同构的问题是否NP-hard未知，但至今没有提出polynomial algorithm。
subgraph isomorphism子图同构如果 $G_2$ 的子图与 $G_1$ 同构，我们称 $G_2$ is subgraph-isomorphic to $G_1$ 。我们可以使用node-induced或edge-induced subgraph定义。这一问题是NP-hrad的。节点之间的映射不必唯一。
subgraph举例
1. 所有非同构的、connected、无向的4个节点的图：
2. 所有非同构的、connected、有向的3个节点的图：一般最多也就4-5个节点了
network motifs 定义：recurring, significant patterns of interconnections
1. pattern：小的node-induced subgraph
2. recurring：出现很多次，即出现频率高（以下将介绍如何定义频率）
3. signifant：比预期（如在随机生成的图中）出现的频率高（以下将介绍如何定义随机图）
motif举例：如图所示：左上角就是我们所感兴趣的induced subgraph（motif）。蓝三角内的induced subgraph符合要求，红三角内不是induced subgraph不符合要求。
motif的意义： 1. 帮助我们了解图的工作机制。 2. 帮助我们基于图数据集中某种subgraph的出现和没有出现来做出预测³。
举例：
1. feed-forward loops：神经元网络中用于中和biological noise
2. parallel loops：食物链中（就两种生物以同一种生物为食并是同一种生物的猎物嘛）
3. single-input modules：基因控制网络⁴ 中
subgraph frequency
1. 图级别的subgraph frequency定义设 $G_Q$ 是一个小图， $G_T$ 是目标图数据集。 $G_Q$ 在 $G_T$ 中的频率： $G_T$ 不同的节点子集 $V_T$ 的数目（ $V_T$ induce的 $G_T$ 的subgraph与 $G_Q$ 同构）如下左图中frequency为2（红圈中的两种节点子集），右图中的frequency为 $C_{100}^6$ （图中的排列组合写法可以参考⁵）
2. 节点级别的subgraph frequency定义设 $G_Q$ 是一个小图， $v$ 是其一个节点（anchor）， $G_T$ 是目标图数据集。 $G_Q$ 在 $G_T$ 中的频率： $G_T$ 中节点 $u$ 的数目（ $G_T$ 的subgraph与 $G_Q$ 同构，其同构映射 $u$ 到 $v$ 上） $(G_Q,v)$ 叫node-anchored subgraph 这种定义对异常值比较鲁棒。如在图例中，star subgraph以中心节点为anchor，其在 $G_T$ 中的frequency就是1；若以其外围节点作为anchor，则其frequency就是100
3. 如果数据集中包含多个图，我们可将其视为一整个大图 $G_T$ （包含disconnected components，各组成部分对应单个图）

2.2 Determining Motif Significance

我们首先需要定义null-model 核心思想：在真实网络中比在随机网络中出现更频繁的subgraph有functional significance 图中论文（配图出处）：Zweig K A . Milo et al. (2002): Network Motifs: Simple Building Blocks of Complex Networks[M]. 2019.
定义随机图：Erdős–Rényi (ER) random graphs $G_{n,p}$ ： $n$ 个节点的无向图，每个边 $(u,v)$ 以频率 $p$ 独立同分布出现。可以是disconnected：
新模型：configuration model 目标：按照给定度数序列 $k_1,k_2,...\ k_N$ 生成随机图。作为网络的null model很有用，可以将真实图和与其具有相同度数序列的随机图作比。 configuration model流程如图所示：对节点上的边进行两两配对，得到最终的结果图（如果出现重边（multiple edges）⁶或自环的情况，由于其罕见，所以可以直接忽略。如图中A-B节点之间出现了double edge，但在最后的结果图中就忽略了，仅作为一条边来处理）
除了像上图那种节点辐条的做法，还可以使用switching方法：对给定图 $G$ ，重复switching步骤 $Q\cdot|E|$ 次： 1. 随机选取一对边A→B，C→D 2. 交换端点⁷ 使A→D，C→B（仅在无重边或自环被产生时进行交换操作）得到randomly rewired graph（与原图的节点度数相同） $Q$ 需足够大（如100）以使这个过程收敛⁸
motif significance overview 检验motif在真实网络中是否比在随机图中overrepresent的步骤：
1. step1：在真实图中数motif的个数
2. step2：产生多个与真实图有相同统计量（如节点数、边数、度数序列等）的随机图，在这些随机图中数motif的个数
3. step3：使用统计指标衡量每一motif是否显著用Z-score
Z-score for statistical significance $Z_i$ 捕获motif $i$ 的statistical significance： $Z_i=(N_i^{real}-\overline{N}_i^{rand})/\text{std}(N_i^{rand})$ 其中 $N_i^{real}$ 是真实图中motif $i$ 的个数， $\overline{N}_i^{rand}$ 是随机图实例中motif $i$ 的平均个数。 0就是说真实图和随机图中motif出现的一样多。绝对值大于2⁹ 时就算显著地多或者显著地少。
network significance profile (SP): $SP_i=Z_i\bigg/\sqrt{\sum_jZ_j^2}$ $SP$ 是归一化的Z-score向量，其维度取决于我们考虑的motif的数量。 $SP$ 强调subgraph的相对重要性：在比较不同大小的网络时很重要，因为一般来说，大图会出现更高的Z-score⁹。
significance profile 对每个subgraph，Z-score指标可用于分类subgraph significance：负数意味着under-representation，整数意味着over-representation。 network significance profile是具有所有subgraph大小上的值的feature vector。接下来就可以比较随机图和不同图上的profile了。不同网络举例：
1. 基因调控网络⁴
2. 神经网络（突触连接）
3. 万维网（网页超链接）
4. 社交网络
5. language networks (word adjacency)¹⁰
significance profile示例：相似领域的网络会具有相似的SP。可以通过motif frequency来判断网络功能。如社交网络中的subgraph6少、但是subgraph13多，因为一个人很难同时与两人保持紧密好友关系而这两个人不互相认识，每周还要出来喝2次咖啡。毕竟如果他们认识以后就可以一周出来只约1次咖啡了。
检测motif总结：
motif概念的变体：衍生：有向/无向，colored/uncolored（应该指的是节点类型，如下图中右上角045算motif出现、345不算motif出现的情况），动态/static motif 概念上的变体：不同的frequency概念、不同的significance指标、under-representation (anti-motifs)（如下图中右下角所示）、不同的null models
motif总结： subgraph和motif是图的组成部分，子图同构和技术问题是NP-hrad。理解数据集中motif的频繁或显著出现，可以使我们了解该领域的独有特征。使用随机图作为null model来通过Z-score衡量motif的显著性。

3. Neural Subgraph Matching / Representations

subgraph matching¹¹ 给出大的target graph（可以是disconnected），query graph（connected）问题：query graph是否是target graph中的子图？示例如下图（节点颜色表示映射关系）：
在本课程中我们不用combinatorial matching、逐边检查，而将该问题视作预测任务，使用机器学习方法来解决这一问题。直觉：利用嵌入域的几何形状来捕获子图同构的属性¹²
task setup 二元预测问题：返回query是否与target graph子图同构（注意在这里我们只关注该预测问题的最终决策，即是不是。而具体的节点之间如何一一对应的关系本课程中不讲）
overview of the approach 整体流程如图所示：将target graph拆成很多neighborhoods，嵌入neighborhoods和query，将每个neighborhood与query做匹配，判断其是否子图同构：
neural architecture for subgraphs
1. 我们将使用node-anchored定义，用anchor的嵌入来判断是否同构
2. 使用node-anchored neighborhoods：上图中应该是把右图的黄色点画成蓝色了
  用GNN基于anchor的邻居得到其嵌入，预测 $u$ 的邻居是否与 $v$ 的邻居同构（图中应该是用了二阶邻居的意思）：
为什么要使用anchor呢？回忆node-level frequency definition。这是因为我们可以用GNN来获得 $u$ 和 $v$ 对应的嵌入，从而可以得知 $u$ 的邻居是否与 $v$ 的邻居同构，这样就可以预测是否存在anchor的映射关系并识别出对应的特定节点。
将 $G_T$ 分解为neighborhoods：对 $G_T$ 中的每个节点（准anchor），获取其k跳邻居（可以通过BFS获取）。k是一个超参，k越大，模型代价越高。同样的过程也应用于 $G_Q$ ，同样得到neighborhoods。我们通过GNN得到anchor node基于neighborhood的embedding，也就是得到了这些neighborhoods的嵌入。
order embeddings space 将图 $A$ 映射到高维（如64维）嵌入域的点 $z_A$ ，使 $z_A$ 所有维度元素都非负。可以捕获partial ordering（关系可传递）（具体见图）：总之可以用嵌入各维元素全部小于等于的关系来表示subgraph⁹
subgraph order embedding space 如图：在order embedding space中全部维度小于target graph的anchor嵌入的就是其subgraph的anchor嵌入（在二维嵌入域中，就是在target graph的左下角）
为什么要使用order embedding space？因为subgraph isomorphism relationship可以很好地在order embedding space中编码，也就是说order embedding可以在向量域表示图域中subgraph的关系：
transitivity：图域：如果 $G_1$ 是 $G_2$ 的subgraph， $G_2$ 是 $G_3$ 的subgraph，则 $G_1$ 是 $G_3$ 的subgraph 嵌入域：（见图）
anti-symmetry：图域：如果 $G_1$ 是 $G_2$ 的subgraph， $G_2$ 是 $G_1$ 的subgraph，则 $G_1$ 与 $G_2$ 同构。嵌入域：（见图）
closure under intersection：图域：1个节点的图是所有图的subgraph 嵌入域：（见图） Corollary推论（必然的结果(或结论);）¹³ 这个valid embedding是啥东西我就没看懂
order constraint 我们用GNN来嵌入neighborhoods并保持其order embedding结构，因此我们需要学习一种可以学习反映subgraph关系的order embedding的损失函数。我们基于order constraint设计损失函数。order constraint规定了理想的可反映subgraph关系的order embedding属性： $\forall_{i=1}^Dz_q[i]\leq z_t[i]\ \text{iff}\ G_Q\subseteq G_T$ 图中应该是把t写成u了其中 $z_q$ 是query embedding， $z_u$ 是target embedding， $i$ 是embedding dimension（D应该是嵌入维度）， $\subseteq$ 是subgraph关系。 order constraint用max-margin loss¹⁴ 来训练。
损失函数 GNN嵌入通过最小化max-margin loss学得。定义 $E(G_q,G_t)=\sum_{i=1}^D\big(\max(0,z_q[i]-z_t[i])\big)^2$ 为图 $G_q$ 和 $G_t$ 之间的margin（penalty或violation）。当margin=0时 $z_q[i]\ ＜\ z_t[i]$ 恒成立，即 $G_q$ 是 $G_t$ 的subgraph；当margin＞0时 $G_q$ 不是 $G_t$ 的subgraph。
training neural subgraph matching 为了学习这种嵌入，需要约束训练集样本 $(G_q,G_t)$ 中一半 $G_q$ 是 $G_t$ 的subgraph，另一半不是。对正样本：最小化 $E(G_q,G_t)$ 对负样本：最小化 $\max\big(0,\alpha-E(G_q,G_t)\big)$ max-margin loss使正样本中 $z_q[i]\ ＜\ z_t[i]$ ，使全式为0；负样本中的 $E(G_q,G_t)$ 大于 $\alpha$ （使 $\alpha-E(G_q,G_t)$ 小于0，使全式为0）。但都不强化其差异，使嵌入向量之间不需要隔很远（slides中说这是degenerate strategy，我也没查到这又是个啥东西）
training example construction 我们需要从数据集 $G$ 中生成训练样本：query $G_Q$ 和target $G_T$
得到 $G_T$ ：随机选取anchor $v$ ，获取其全部 $K$ 阶邻居。
得到 $G_Q$ ：用BFS抽样，从 $G_T$ 中抽样induced subgraph： step1：初始化 $S=\{v\},V=\varnothing$ step2：使 $N(S)$ 为 $S$ 中节点的所有邻居。每一步抽样10%在 $N(S)\backslash V$ 中的节点放入 $S$ 中，并将其余节点放在 $V$ 中。 step3： $K$ 步后，获取 $G$ 的 induced by $S$ anchored at $v$ 的subgraph 图中应该是把v写成q了
对负样本（ $G_Q$ 不是 $G_T$ 的subgraph）：corrupt $G_Q$ ：增加/移动节点/边使它不再是一个subgraph
训练细节
1. 我们需要抽样出多少training examples？每次迭代，我们都需要抽样新的target-query对。这样做的好处是每次模型都会看到不同的subgraph例子，提升表现结果、避免过拟合（毕竟有指数级的可能subgraph可供抽样）
2. BFS抽样应该多深？这是一个需要平衡运行时间和表现结果的超参数，一般设置为3-5，也需要取决于数据集的大小。
在新图上预测一个图是否是另一个图的subgraph 已知：query graph $G_q$ anchored at node $q$ , target graph $G_t$ anchored at node $t$ 目标：输出query是否是target的node-anchored subgraph 过程：如果 $E(G_q,G_t)<\epsilon$ ，预测为真；反之预测为假（ $\epsilon$ 是超参数）为了检验 $G_Q$ 是否与 $G_T$ subgraph isomorphism，对所有 $q\in G_Q,t\in G_T$ 重复上述流程。此处的 $G_q$ 是 $q\in G_Q$ 附近的neighborhood。
neural subgraph matching总结 neural subgraph matching使用基于机器学习的方法学习NP-hard的subgraph isomorphism问题： 1. 将query和target图都嵌入order embedding space 2. 用这些嵌入向量计算 $E(G_q,G_t)$ 以判断query是否是target的subgraph
用order embedding space嵌入图使subgraph isomorphism可以高效表示并由图嵌入的相对位置进行检验。

4. Mining / Finding Frequent Motifs / Subgraphs

finding frequent subgraphs 找最频繁的大小为k的motif需要解决两个挑战：
1. 迭代所有大小为k的connected subgraph
2. 数每一类subgraph的出现次数
这个问题难在仅确定一个特定subgraph是否存在于图中就已是计算代价很高的问题（subgraph isomorphism是NP-complete），计算用时随subgraph指数级增长（因此传统方法可行的motif尺寸都相对较小，如3-7）。可以说这是两个指数级增长的问题梦幻联动（特定大小有多少motif（组合爆炸combinatorial explosion¹⁵）+找出每个motif的frequency（subgraph isomorphism and subgraph counting））
使用表示学习的方法来解决问题表示学习通过search space（每次增加一个节点，累积到size k的subgraph上，详情见下文。注意我们仅关心高频subgraph）解决组合爆炸问题，通过GNN预测解决subgraph isomorphism问题（就是本章第三节所讲述的知识）。
problem setup: frequent motif mining 给出target graph（数据集） $G_T$ ，subgraph大小参数 $k$ 所需结果数量 $r$ 目标：从所有大小为 $k$ 个节点的图中，识别出 $r$ 个在 $G_T$ 中出现频率最高的图。我们使用node-level subgraph定义。
SPMiner overview SPMiner：识别高频motifs的神经网络模型步骤：将输入图 $G_T$ decompose为重复的node-anchored neighborhoods，将subgraph嵌入到order embedding space（上述两步和neural subgraph matching是一样的）；然后进行search procedure，策略是不遍历所有可能的subgraph，而直接从2个节点的图开始增长出一个所需节点数的subgraph，在增长的同时尽量保证频率高。
SPMiner：核心思想 order embedding的核心优势：可以迅速找到特定subgraph $G_Q$ 的频率
motif频率估计已知：一系列 $G_T$ 的subgraph（node-anchored neighborhoods） $G_{N_i}$ （通过随机抽样得到）。核心思想：估计 $G_Q$ 的频率：通过数符合要求的 $G_{N_i}$ 的个数（ $z_Q\leq z_{N_i}$ ）这是由order embedding space的属性得到的结论：如下图所示，红色节点（motif）右上角的红色区域就是其super-graph region，红色节点是所有落在该区域的节点（ $G_T$ 的neighborhoods）的subgraph。这样的好处就是算得快。
SPMiner search procedure
1. 开始：从一个从target graph中随机选择的起点 $u$ 开始：设 $S=\{u\}$ 所有neighborhoods都在一个点的右上方区域，即都包含这个subgraph。
2. 迭代：每次选一个 $S$ 中节点的邻居，加到 $S$ 中，如此逐渐增长motif的尺寸。目标：在 $k$ 步后最大化红色阴影区域中的neighborhoods数目。
3. 停止：达到预期motif尺寸后，选取the subgraph of the target graph induced by $S$ 我们找到的motif就是预期尺寸备选subgraph嵌入中有最多target graph neighborhoods（蓝点）在红色区域的subgraph。
4. 每一步如何选取节点？定义subgraph $G$ 的total violation：不包含 $G$ 的neighborhoods数量。即不满足 $z_Q\preccurlyeq z_{N_i}$ 的neighborhoods $G_{N_i}$ 数量。最小化total violation就是最大化频率。我们采用贪心算法，每一步都选择时当前total violation最小的节点。
实验结果
1. 小motif ground truth：通过代价高昂的BF迭代算法（暴力破解）找到10个出现频率最高的motif。在大小为5-6的motif上，SPMiner可以识别出top 10中前9/8个，识别出的频率接近真实值：
2. 大motif SPMiner比baseline多识别10-100倍。
总结
1. subgraph和motif是可用于深入了解图结构的重要概念，对其计数可用作节点或图的特征。
2. 本章介绍一种预测subgraph isomorphism关系的神经网络方法。
3. order embeddings的特性可用于编码subgraph关系。
4. order embedding space上的neural embedding-guided search让我们有了一种比传统方法能识别出更高motif频率的机器学习模型。

我直接用百度查这个关键词没查到，只查到说离散数学中有matching problem（参考这篇文章：离散数学笔记6（Matching problems） - 知乎），但是感觉这个跟本章内容没什么关系？所以我觉得原文指的应该是传统的subgraph matching方法¹¹。对于这些传统的subgraph matching方法我没有做过了解，我就是直接学的本章课程。 ↩
Edge-induced subgraph | Article about edge-induced subgraph by The Free Dictionary 说edge-induced subgraph中的节点是全部在其边中出现过至少一次的原图节点。 Edge-Induced Subgraph -- from Wolfram MathWorld: An edge-induced subgraph is a subset of the edges of a graph G together with any vertices that are their endpoints⁷. 子图的概念_yyywww666的专栏-CSDN博客 networkx.classes.function.edge_subgraph — NetworkX 2.5 documentation dgl.edge_subgraph — DGL 0.6.1 documentation ↩
看图上原句……呃，不知道有没有理解并翻译正确。不过我觉得这点正误小问题不重要吧。 presence一词不知道在此处有无专业术语层面的引申义存在。 ↩
可参考基因调控网络 (Gene Regulatory Network) 01 - Bracer - 博客园因为不是我的专业领域，所以我没仔细看。 ↩ ↩²
我之前写的笔记：cs224w（图机器学习）2021冬季课程学习笔记2: Traditional Methods for ML on Graphs_诸神缄默不语的博客-CSDN博客中脚注4 ↩
当时我做这笔记的时候还能上英文维基，但是现在暂时不能了，所以我就把我当时用的参考资料列出来，具体内容只能看我之前笔记了，我没法再去原网页看一眼以作证实了： ① 一对节点之间连了多个边来源：zh.wikipedia.org/wiki/%E9%87… ② en.wikipedia.org/wiki/Dual_g… ↩
一条边的endpoint就是其对应的两个端点。无论有向还是无向都是两个端点。参考自OI Wiki。原话：图源：图论相关概念 - OI Wiki ↩ ↩²
所以这个玩意居然能收敛吗也是很神奇 ↩
没搞懂为什么 ↩ ↩² ↩³
这个好像挺复杂的，我也没仔细看。链接反正也发出来，以资参考：Networks based on words ↩
我在之前查到的对subgraph matching问题的定义：给出target graph和query graph，确定query graph是否是target graph的子图，如是确定其位置。参考¹⁶）其他可参考定义及来源： Subgraph matching problem is identifying a target subgraph in a graph. 来自 Subgraph Matching Using Graph Neural Network ↩ ↩²
我并没有这样的直觉，也没看懂这句话啥意思。 ↩
可以参考：区分定理(Theorem)、引理(Lemma)、推论(Corollary)等概念_cloudeagle_bupt的专栏-CSDN博客_corollary ↩
我是感觉max margin在网上我找到的定义跟这边这个不太一样啊……在这边的话就按照课程中讲的来吧。我在网上了解该概念所使用的资料也放出来，以资参考：损失函数：Hinge Loss（max margin） - 菜鸡一枚 - 博客园 ↩
其实我没看懂组合爆炸这个词是啥意思……就我能理解它随着size变大，subgraph变多，这件事。但是组合爆炸这个词本身是啥意思，我就没搞懂。我在stack overflow上看到了这样一个问题：Examples for combinatorial explosion in Java? - Stack Overflow 就是Java等编程语言里面的组合爆炸，这个我能看得懂，就……我也是程序员，我也是解耦人啊，我怎么能不懂interface的意义！但是这跟subgraph又有什么关系？ ↩
原文：Subgraph matching is the problem of determining the presence and location (s) of a given query graph in a large target graph. 来自 [2007.03092] Neural Subgraph Matching ↩

cs224w（图机器学习）2021冬季课程学习笔记15 Frequent Subgraph Mining with GNNs