CN114694832A

CN114694832A - 一种确定阿尔茨海默症的风险要素的电子设备

Info

Publication number: CN114694832A
Application number: CN202210412137.2A
Authority: CN
Inventors: 毕夏安; 王雅琴; 岳凌霞; 吴昊; 徐露允
Original assignee: Hunan Normal University
Current assignee: Hunan Normal University
Priority date: 2022-04-19
Filing date: 2022-04-19
Publication date: 2022-07-01

Abstract

本发明公开了一种确定阿尔茨海默症的风险要素的电子设备，包括存储器、处理器，所述处理器执行计算机程序时用于实现：将多个被试者的脑区‑基因网络通过结构信息聚集模型进行结构信息聚集并提取边的特征后输出每个被试者为阿尔茨海默症AD患者的概率的判断结果；根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区‑基因网络中各节点的重要度；根据节点的重要度，确定出所述脑区‑基因网络中作为AD风险要素的节点；其中，所述被试者包括多个正常人和AD患者。应用本发明可以更精确地确定AD的病变脑区和风险基因，以便助力靶向药物开发和针对性经颅磁治疗。

Description

一种确定阿尔茨海默症的风险要素的电子设备

技术领域

本发明涉及计算机技术领域，特别是指一种确定阿尔茨海默症的风险要素的电子设备。

背景技术

阿尔茨海默症(AD)以一系列显著的痴呆症状为特征，常见于老年人。AD的发病机制复杂，其发病因素横跨宏观到微观多个层次，例如，大脑功能与结构的病变与基因的异常表达密切相关，因此，找到大脑宏观病变与微观变异的关联模式有助于揭示阿尔茨海默症的多层次发病机制。一些研究人员构建网络，将基因或大脑区域抽象为节点，以研究致病因素之间的功能相关性。Huang等人定义了每个脑区的嵌入表示，并提出了一种新的度量评估两个区域之间的相似性，使测量结果更加合理。Lella等人在分析大脑网络时引入了沟通距离，提供了新的3种方法发现与AD发展相关的隐藏模式。

然而，现有技术基于目前方法分析得到的脑区网络和基因网络进行AD诊断的准确率比较低，有必要提供一种更为精准确定AD的风险基因和病变脑区方法，以便助力对应的靶向药物开发和针对性经颅磁治疗。

发明内容

有鉴于此，本发明的目的在于提出一种确定阿尔茨海默症的风险要素的电子设备，可以更为精准确定AD的风险要素，即AD的风险基因和病变脑区，以便助力对应的靶向药物开发和针对性经颅磁治疗。

基于上述目的，本发明提供一种确定阿尔茨海默症的风险要素的电子设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时用于实现：

将多个被试者的脑区-基因网络通过结构信息聚集模型进行结构信息聚集并提取边的特征后输出每个被试者为阿尔茨海默症AD患者的概率的判断结果；

根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各节点的重要度；

根据节点的重要度，确定出所述脑区-基因网络中作为AD风险要素的节点；

其中，所述被试者包括多个正常人和AD患者；以及所述结构信息聚集模型是以多个正常人和AD患者的脑区-基因网络作为样本预先训练得到的。

更优地，所述处理器执行所述计算机程序还用于实现：

根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各边的重要度；

根据边的重要度，确定出所述脑区-基因网络中作为AD风险要素的病变脑区和风险基因之间的关联关系。

其中，所述结构信息聚集模型包括：

生成器，用于对当前输入的被试者的脑区-基因网络进行结构信息的聚集后得到第一关键结构子网络，进而对第一关键结构子网络进行结构信息的扩散得到重构的脑区-基因网络；

判别器，用于对重构的脑区-基因网络进行结构信息的聚集后得到第二关键结构子网络，进而提取第二关键结构子网络的边的特征后输出该被试者为AD患者的概率的判断结果。

可选的，所述生成器中包括：

第一卷积模块，用于通过n阶双向十字卷积层对当前输入的被试者的脑区-基因网络进行结构信息的聚集后得到第一关键结构子网络；

反卷积模块，用于对第一关键结构子网络通过n阶扩散反卷积层进行结构信息的扩散得到重构的脑区-基因网络。

可选的，所述判别器中包括：

第二卷积模块，用于通过n阶双向十字卷积层对重构的脑区-基因网络进行结构信息的聚集后得到第二关键结构子网；

全连接层模块，用于提取第二关键结构子网络的边的特征后输出该被试者为AD患者的概率的判断结果。

较佳地，所述根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各节点的重要度，具体包括：

根据所述结构信息聚集模型输出的判断结果与真实情况的接近程度，以及所述结构信息聚集模型中第一/二关键结构子网络的权重矩阵，计算关键特征矩阵；

根据计算的关键特征矩阵，确定AD患者的脑区-基因网络中各节点的重要度；

根据节点的重要度，确定出所述脑区-基因网络中作为AD风险要素的节点。

较佳地，所述根据所述结构信息聚集模型输出的判断结果与真实情况的接近程度，以及所述结构信息聚集模型中第一/二关键结构子网络的权重矩阵，计算关键特征矩阵，具体包括：

根据如下公式17计算关键特征矩阵FW：

其中，KWⁱ表示所述结构信息聚集模型针对第i个被试者的脑区-基因网络生成的第一/二关键结构子网络的权重矩阵；Probⁱ表示所述结构信息聚集模型针对第i个被试者输出的判断结果和第i个被试者的真实情况的接近程度；Probⁱ×KWⁱ表示将第i个被试者的预测概率分数加权给权重矩阵KWⁱ的每个位置上；m表示被试的个数。

较佳地，所述根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各边的重要度，具体包括：

根据计算的关键特征矩阵，确定AD患者的脑区-基因网络中各边的重要度；

进一步，所述中央处理单元执行所述程序还实现：对所述所述结构信息聚集模型进行训练：

根据所述样本对所述生成器和判别器进行多次迭代的优化和训练；其中，在一次迭代过程中，包括：

固定所述生成器的参数后，根据所述样本对所述判别器进行优化；

固定所述判别器的参数后，根据所述样本对所述生成器进行训练；

在对所述生成器和判别器进行多次迭代的优化和训练后，得到由所述生成器和判别器构成的、训练后的结构信息聚集模型。

本发明还提供一种计算机可读存储介质，所述计算机可读存储介质内存储有计算机程序，所述计算机程序可被至少一个处理器所执行，以使所述至少一个处理器执行上述的电子设备的操作。

本发明的技术方案中，将多个被试者的脑区-基因网络通过结构信息聚集模型进行结构信息聚集并提取边的特征后输出每个被试者为阿尔茨海默症AD患者的概率的判断结果；根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各节点的重要度；根据节点的重要度，确定出所述脑区-基因网络中作为AD风险要素的节点；其中，所述被试者包括多个正常人和AD患者；以及所述结构信息聚集模型是以多个正常人和AD患者的脑区-基因网络作为样本预先训练得到的。由于结构信息聚集模型可以对输入的脑区-基因网络进行结构信息的聚集，使得正常人与AD患者的脑区-基因网络经过结构信息聚集后得到的关键结构子网络差异显著，从而能够基于关键结构子网络更为精准确定AD的风险要素，即AD的风险基因和病变脑区，以便助力对应的靶向药物开发和针对性经颅磁治疗。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提供的结构信息聚集过程的示意图；

图2为本发明实施例一提供的结构信息聚集模型的内部结构框图；

图3为本发明实施例一提供的生成器的内部结构示意图；

图4、5为本发明实施例一提供的n阶双向十字卷积的操作示意图；

图6为本发明实施例一提供的结构信息扩散和吸收之间的关系示意图；

图7、8为本发明实施例一提供的n阶双向扩散反卷积的操作示意图；

图9为本发明实施例一提供的判别器的内部结构示意图；

图10为本发明实施例一提供的全连接层的内部结构示意图；

图11为本发明实施例二提供的结构信息聚集模型的训练方法流程图；

图12为本发明实施例二提供的对所述生成器和判别器进行一次迭代的优化和训练的方法流程图；

图13为本发明实施例二提供的对判别器进行优化的方法流程图；

图14为本发明实施例二提供的对生成器进行训练的方法流程图；

图15为本发明实施例三提供的一种基于结构信息聚集模型的阿尔茨海默症的判断方法的流程图；

图16为本发明实施例四提供的一种确定阿尔茨海默症的风险要素的方法流程图；

图17为本发明实施例四提供的提取病变脑区和致病基因及其关联关系的过程示意图；

图18为本发明实施例四提供的电子设备硬件结构示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。

需要说明的是，除非另外定义，本发明实施例使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本公开中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也可能相应地改变。

本发明的发明人发现，在疾病演化为AD的这个过程中，大脑发生了质的变化，例如，部分脑区出现萎缩，功能连接出现大量异常。大脑的这些宏观病变与微观层面基因的异常表达密切相关，可以揭示从宏观到微观发病因素的关联机制。在基于脑影像(fMRI)数据和基因数据构建的脑区-基因网络中，关键结构子网络的结构信息表征着网络所独有的结构特征信息，主要体现在网络中边的权重，脑区或基因越异常，则连接它们的边就越能够代表该网络的结构特征信息。为了更好地探究正常人和AD患者的脑区-基因网络的关键结构子网络，发现正常人和AD患者的区别性特征，本发明设计了结构信息聚集模型，又称为结构聚集生成对抗网络(Structure GatheringGAN，SG-GAN)。模型的设计思想是：对于脑区-基因网络中的每条边，通过聚集其相邻边的结构信息，更新自身边结构信息，经过n次结构信息聚集后，边的结构信息不断更新可以得到关键结构子网络。由正常人和AD患者的脑区-基因网络分别得到的关键结构子网络具有更显著的差异，从而能够更为精准确定AD的风险要素，即AD的风险基因和病变脑区，以便助力对应的靶向药物开发和针对性经颅磁治疗。

下面结合附图详细说明本发明实施例的技术方案。

对脑区-基因网络进行结构信息聚集的过程，如图1所示。在依据被试者的脑影像数据和基因数据构建脑区-基因网络后，对所述脑区-基因网络中的边经过多次结构信息聚集能够不断获得邻域中边的高阶结构信息，使得网络的关键结构信息逐渐聚集到部分的边上；聚集的过程能够使代表边结构信息的权重增大或者减小，网络的拓扑结构也随之不断变化，原始的脑区-基因网络逐渐演变为关键结构子网络。也就是说，模型能够在脑区-基因网络不断演变的过程中将疾病的特征信息逐渐汇聚到关键结构子网络上，因此，根据关键结构子网络能够更为精准确定AD的风险要素，即AD的风险基因和病变脑区，以便助力对应的靶向药物开发和针对性经颅磁治疗。

下面针对脑区-基因网络建立结构信息聚集模型，模型的要素包括如下：

(1)脑区-基因网络：Gⁱ表示第i个被试者(NC或AD)的脑区-基因网络，由第i个被试者的脑影像数据和基因数据构建。

(2)边E_pq的权重：

边E_pq表示的是节点p和q之间的边，表示第i个被试者的脑区-基因网络初始状态时边E_pq的权重，通过计算节点p和q之间的皮尔逊相关系数得到。

表示第i个被试者第n次结构信息聚集后边E_pq的权重，权重越大表明节点p和q之间的联系程度越紧密。其中，N为脑区-基因网络中节点的个数。

(3)被试者的关键结构子网络：

KGⁱ表示第i个被试者的脑区-基因网络经过n次结构信息聚集得到第i个被试者的关键结构子网络，

表示第i个被试者的关键结构子网络中边E_pq的权重。通常某一类疾病有多个被试者，每个被试者都有对应的关键结构子网络，即一类疾病有多个关键结构子网络。

根据以上模型要素，被试者的脑区-基因网络经过结构信息聚集模型得到关键结构子网络可分为以下两个步骤：(1)脑区-基因网络中边的权重经过n次结构信息聚集得到更新的边权重；(2)根据更新的边权重构建被试的关键结构子网络。

第一步：脑区-基因网络中边的权重经过n次结构信息聚集得到更新的边权重。针对每一个被试的脑区-基因网络，在每次结构信息聚集过程中，脑区-基因网络中每条边都会聚集其邻域内边的结构信息，进而更新自身边的权重大小，经过多次结构信息聚集后，原始的脑区-基因网络逐渐演变为关键结构子网络。在第n次结构信息聚集过程中，第i个被试者的脑区-基因网络中边的权重W^(i，n)的更新公式如公式1所示：

其中，β是节点聚集邻域内边结构信息的比例，

表示边E_pq聚集节点p所连接边的结构信息，

表示边E_pq聚集节点q所连接边的结构信息；

表示第i个被试者的脑区-基因网络中边E_pq聚集其邻域中边的结构信息之和。

第二步：根据更新的边权重构建被试者的关键结构子网络。对于第i个被试的脑区-基因网络Gⁱ，初始边的权重W^(i，0)经过n次结构信息聚集后更新为W^(i，n)，根据W^(i，n)可以得到第i个被试者的关键结构子网络KGⁱ，如公式2所示：

KGⁱ＝Graph(W^(i，n)) (公式2)

其中，Graph(W^(i，n))表示根据边的权重W^(i，n)构建关键结构子网络。

经过n次结构信息聚集，模型能够捕获到正常人和AD患者脑区-基因网络之间关键结构子网络的差异，完成被试者的分类任务，如公式3所示：

Pⁱ＝Classify(KGⁱ) (公式3)

其中，KGⁱ表示第i个被试这的关键结构子网络；Classify(·)表示分类函数；Pⁱ表示第i个被试者的判断结果为AD患者或者正常人的概率，是一个大小为1×2的特征向量。

最后，根据判断结果提取正常人和病人之间具有显著性差异的病变脑区和致病基因以及它们的关联模式，其特征提取的公式如公式4所示：

(PG，PBR，AM)＝Feature_extract(Pⁱ，KGⁱ)i＝1，...，m (公式4)

其中，KGⁱ表示第i个被试的关键结构子网络；Pⁱ表示第i个被试者判断结果为AD患者或者正常人的概率；Feature_extract(·)表示特征提取函数；PG表示风险基因，PBR表示病变脑区，AM表示病变脑区和风险基因的关联模式。

综上所述，将公式1～4联立，用以下方程组来表示脑区-基因网络的结构信息聚集模型和应用：

公式1表示脑区-基因网络中的边聚集其邻域内边的结构信息之和，公式2表示根据第n次结构信息聚集后得到的边权重W^(i，n)构建第i个被试的关键结构子网络，公式3表示根据被试的关键结构子网络对正常人和AD患者分类，公式4表示利用分类结果和关键结构特征信息提取AD患者的病变脑区和致病基因以及它们的关联模式。

实施例一

基于上述的结构信息聚集模型的要素，本发明实施例一提供了一种基于生成对抗网络的结构信息聚集模型，其内部结构如图2所示，包括：生成器201和判别器202；

其中，生成器201用于对所述脑区-基因网络进行结构信息的聚集后得到第一关键结构子网络，进而对第一关键结构子网络进行结构信息的扩散得到重构的脑区-基因网络。

判别器202用于对重构的脑区-基因网络进行结构信息的聚集后得到第二关键结构子网络，进而提取第二关键结构子网络的边的特征后输出该被试者为AD患者的概率的判断结果。

在由生成器201和判别器202组成的结构信息聚集模型中，生成器201包含卷积部分和反卷积部分，卷积部分的双向十字卷积对真实脑区-基因网络进行卷积，从而提取关键结构子网络；针对关键结构子网络，反卷积部分中的双向扩散反卷积用于重构出生成脑区-基因网络。最后，判别器202对输入的生成脑区-基因网络和真实脑区-基因网络进行判断和分类，并根据其结果不断更新生成器或判别器的参数，直到真实脑区-基因网络和生成脑区-基因网络的差异最小并且判别器能够准确判断输入网络的类别，即整个生成对抗网络收敛。此时，训练后的结构信息聚集模型能够精准地对AD进行分类诊断。

在一个较优的实施方式中，如图3所示，生成器201包括卷积部分和反卷积部分；具体地，生成器201中包括：作为卷积部分的第一卷积模块211，以及作为反卷积部分的反卷积模块212。

作为卷积部分的第一卷积模块211用于对所述脑区-基因网络进行结构信息的聚集后得到第一关键结构子网络；

作为反卷积部分的反卷积模块212用于对第一关键结构子网络进行结构信息的扩散得到重构的脑区-基因网络。

如图3所示，生成器201由卷积部分和反卷积部分组成，卷积部分完成脑区-基因网络结构信息的聚集，进而得到关键结构子网络，反卷积部分完成关键结构子网络结构信息的扩散，从而重构出生成脑区-基因网络。下面设计卷积和反卷积的矩阵运算，描述输入的脑区-基因网络经过结构信息聚集得到关键结构子网络，以及根据关键结构子网络重构脑区-基因网络的过程。

生成器201的输入具体为脑区-基因网络的权重矩阵，该权重矩阵W⁽⁰⁾∈R^N×N的构建方式如下：通过皮尔逊相关系数计算得到，定义Pear_i，j为节点i和j的皮尔逊相关系数，N为脑区-基因网络的节点个数(N＝N₁+N₂)，N₁、N₂分别为脑区-基因网络中脑区节点个数、基因节点个数；W_ij表示节点i和j之间边的权重，则权重矩阵W⁽⁰⁾可如公式5所示：

也就是说，输入生成器201的脑区-基因网络的权重矩阵W⁽⁰⁾中，第i行和第j列的元素即为脑区-基因网络中节点i和j之间边的权重；而节点i和j之间边的权重是通过皮尔逊相关系数计算得到的Pear_i，j。

在一个较优的实施方式中，作为卷积部分的第一卷积模块211中包括n阶双向十字卷积层；每经过一次双向十字卷积层的卷积操作，脑区-基因网络的权重矩阵都会得到更新，进而使脑区-基因网络的拓扑结构发生变化。下面以k-1阶脑区-基因网络经过第k阶双向十字卷积层的卷积操作，得到k阶脑区-基因网络为例，描述双向十字卷积的具体过程。

其中，第k阶双向十字卷积层的输入为W^(k-1)，W^(k-1)是基于构建的脑区-基因网络的权重矩阵W⁽⁰⁾经过k-1阶双向十字卷积层的k-1次双向十字卷积后得到的权重矩阵，即k-1阶脑区-基因网络的权重矩阵；W^(k-1)经过第k阶双向十字卷积层的双向十字卷积操作后，每条边聚集了k阶邻域内边的结构信息，得到更新的权重矩阵W^(k)，即k阶脑区-基因网络的权重矩阵，为第k阶双向十字卷积层的输出；W^(k)的计算公式如公式6所示：

其中，b^(k)表示第k阶双向十字卷积层的偏置；A^(k)∈R^N×N表示第k阶双向十字卷积层的总核；r和c分别表示矩阵的行、列索引值；

表示权重矩阵W^(k-1)中第r行第i列的值和A^(k)中第r行第i列的值对应相乘；f(·)是激活函数ReLU(·)；

是

经过双向十字卷积层的卷积操作更新后的值。

亦即，第k阶双向十字卷积层的输入为k-1阶脑区-基因网络的权重矩阵W^(k-1)，输出为k阶脑区-基因网络的权重矩阵权重矩阵W^(k)。其中，n为大于1的自然数，k为小于等于n的自然数；0阶脑区-基因网络的权重矩阵W⁽⁰⁾，为输入生成器201的脑区-基因网络的权重矩阵W⁽⁰⁾。

图4展示了一个具有6个节点的脑区-基因网络解释n阶双向十字卷积操作。图4中的(a)图E₁₃为将要进行n阶双向十字卷积操作的边，E₁₂、E₁₄、E₁₆、E₃₂和E₃₄为该边一阶邻域内的边，此时这些边已经聚集了它们n-1阶邻域内边的结构信息。图4中的(b)图是根据即将卷积的边从卷积总核选取对应的双向十字卷积核。根据即将卷积的边E₁₃，选取n阶双向十字卷积总核A⁽ⁿ⁾∈R^6×6中的第1行第3列和第3行第1列作为双向十字卷积核。图4中的(c)图是对边E₁₃的n阶双向十字卷积操作。根据上一步选取的双向十字卷积核(第1行第3列和第3行第1列卷积核)，将它们分别作为权重一一对应赋予给权重矩阵W^(n-1)∈R^6×6的第1行第3列和第3行第1列，再进行加权求和得到边E₁₃相应的值，完成对边E₁₃的n阶双向十字卷积操作。以此类推，对每条边进行n阶双向十字卷积操作，得到更新的权重矩阵W⁽ⁿ⁾∈R^6×6。脑区-基因网络的权重矩阵W⁽⁰⁾经过n阶双向十字卷积层的卷积操作进行了n次更新，使脑区-基因网络的拓扑结构发生n次变化得到W⁽ⁿ⁾；如上述公式2所示，根据更新的W⁽ⁿ⁾即可得到关键结构子网络；也就是说，更新的权重矩阵W⁽ⁿ⁾体现了脑区-基因网络中每条边聚集的其n阶邻域内边的结构信息，体现了关键结构子网络的拓扑结构。

n阶双向十字卷积的意义：在n阶双向十字卷积过程中，n-1阶脑区-基因网络中的每条边能够聚集其一阶邻域内边的结构信息，而该边一阶邻域内的边已经聚集它们n-1阶邻域内边的结构信息，所以此时网络中每条边能够聚集其n阶邻域内边的结构信息。也就是说，每条边n阶邻域内代表疾病特征信息的结构信息通过节点传递给该边。此外，由于所构建的权重矩阵是无向的，而权重W_ij是由节点i和j之间的两条有向连接(i→j和j→i)共同决定的，双向十字卷积可以聚集存在无向边中的有向边结构信息，从而可以更准确地捕捉到关键结构信息。在AD患者的脑区-基因网络中，边将异常的脑区或基因连接起来形成疾病独有的结构信息，脑区或基因越异常，连接它们的边就越能够代表该网络特征的结构信息。本发明设计的n阶双向十字卷积能够很好地捕捉脑区-基因网络中边与其n阶邻域内边的联系，从而使每条边都能聚集其n阶邻域内边的结构信息。图5表示的是n-1阶脑区-基因网络的n阶双向十字卷积操作，此时网络中的每条边都聚集了n-1阶邻域内的结构信息。以图中边E₁₂为例，红色虚线框区域为边E₁₂的一阶邻域，经过n阶双向十字卷积操作，边E₁₂聚集了其n阶邻域内边的结构信息。

基于卷积部分的设计思路，本发明提供的作为反卷积部分的反卷积模块212中包括n阶双向扩散反卷积层，实现从关键结构子网络重构脑区-基因网络的过程。反卷积模块212输入的是W⁽ⁿ⁾，即关键结构子网络的权重矩阵；如图3所示，反卷积部分展示的是双向扩散反卷积的过程。每经过一次双向扩散反卷积操作，便重构一次脑区-基因网络，重构的脑区-基因网络的权重矩阵会更新，亦即重构的脑区-基因网络的拓扑结构发生变化。

下面以k-1阶重构脑区-基因网络经过k阶双向扩散反卷积得到k阶重构脑区-基因网络为例，描述双向扩散反卷积的具体过程。

在脑区-基因网络中，如果两个节点之间没有边连接，称其被一条潜在边所连接；对于网络中存在的边，称该边为真实边。关键结构子网络中每条真实边将其自身的结构信息扩散到其n阶邻域内的真实边和潜在边上，进而使得网络中边的权重发生变化。对于接收结构信息的真实边和潜在边，真实边吸收了其n阶邻域边扩散的结构信息，从而更新边的权重；潜在边通过连接的两个节点吸收了其n阶邻域内边扩散的结构信息，当吸收的信息达到一定阈值时，潜在边就会转换为真实边。

如图6所示，以关键结构子网络的一阶边结构信息扩散为例，展示网络中边结构信息的扩散和边吸收结构信息之间的关系。图6中(a)图展示的是边E₁₂一阶邻域内边的结构信息扩散的过程，图中标紫的边(E₂₆、E₂₃、E₁₆、E₁₅、E₂₄)为边E₁₂的一阶邻域内的真实边，当这些边的结构信息扩散到它们的一阶邻域上时，真实边E₁₂吸收了每条边扩散的结构信息，即黑色箭头所示。为了更好地展示扩散的过程，可以将真实边E₁₂一阶邻域内的所有边扩散到它的结构信息，看作为真实边E₁₂吸收其一阶邻域内边扩散的结构信息，如图6中(b)图所示。在结构信息扩散的过程中，待更新边的结构信息变化量等于该边一阶邻域内所有真实边扩散到该边的总和，由此可见该过程可以用卷积的矩阵运算来描述。

第k阶双向扩散反卷积层的输入为W^′(k-1)，W^′(k-1)是由上述经过n次双向十字卷积得到的权重矩阵W⁽ⁿ⁾经过k-1阶双向扩散反卷积层的k-1次双向扩散反卷积得到的，为k-1阶重构脑区-基因网络的权重矩阵。

权重矩阵W^′(k-1)经过第k阶双向扩散反卷积层的双向扩散反卷积操作后，每条边(真实边和潜在边)吸收了其k阶邻域内边扩散的结构信息，得到更新的权重矩阵W^′(k)，即k阶重构脑区-基因网络的权重矩阵；W^′(k)的计算公式如公式7所示：

其中，b^′(k)是第k阶双向扩散反卷积层的偏置；B^(k)∈R^N×N是第k阶双向扩散反卷积层的总核；r和c分别是矩阵的行、列索引值；

表示权重矩阵W^’(n-1)中第r行第i列的值和总核B^(k)中第r行第i列的值对应相乘；

表示权重矩阵W^′(k-1)中第i行第c列的值和总核B^(k)中第i行第c列的值对应相乘；f(·)是激活函数ReLU(·)；

是

轻过第k阶双向扩散反卷积层的双向扩散反卷积更新后的值。

亦即，第k阶双向扩散反卷积层的输入为k-1阶重构脑区-基因网络的权重矩阵W^′(k-1)，输出为k阶重构脑区-基因网络的权重矩阵W^(k)。

图7展示了一个具有6个节点的关键结构子网络解释n阶双向扩散反卷积操作。图7中(a)的E₁₂为将要进行n阶双向扩散反卷积操作的真实边，图中标黄的真实边E₂₃、E₂₄、E₁₄、E₁₅、E₁₆和潜在边E₁₃、E₂₅、E₂₆为该真实边的一阶邻域内的边，此时这些真实边和潜在边已经吸收了它们n-1阶邻域内边的结构信息。图7中(b)是根据即将反卷积的边从双向扩散反卷积总核选取对应的双向扩散反卷积核。根据即将反卷积的真实边E₁₂，选取n阶双向扩散反卷积总核B⁽ⁿ⁾中的第1行第2列的值K₁₂作为n阶双向扩散反卷积核。图7中(c)是对真实边E₁₂的n阶双向扩散反卷积操作。根据上一步选取的双向扩散反卷积核K₁₂，将其作为权重赋予给权重矩阵W^′(n-1)∈R^6×6的第1行第2列的值

再将其扩散到边E₁₂的一阶邻域内的边上，完成对边E₁₂的n阶双向扩散反卷积操作。以此类推，对每条边(真实边和潜在边)进行n阶双向扩散反卷积操作，得到更新的权重矩阵W^′(n)∈R^6×6。

n阶双向扩散反卷积的意义：在n阶双向扩散反卷积的过程中，n-1阶重构脑区-基因网络中的每条边已经吸收了n-1阶邻域内边扩散的结构信息，此时网络中的每条边将自身结构信息扩散到其一阶邻域内的边上，也就是网络中每条边能够吸收到其n阶邻域内边的结构信息。也就是说，n阶双向扩散反卷积使得网络中每条边n阶邻域内边的结构信息发生变化，即n阶邻域内边的权重增大或者减小。在关键结构子网络中，通过边将异常的脑区或基因连接起来形成疾病信息的关键拓扑结构，网络中的边将关键结构信息扩散到其n阶邻域内的边上，得到扩散了关键结构信息的n阶生成脑区-基因网络。本发明设计的n阶双向扩散反卷积能够很好捕捉关键结构子网络中边与其邻域内边的关系，从而使每条边都将自身边结构信息扩散到其n阶邻域内的边上。图8表示的是n-1阶重构脑区-基因网络的n阶双向扩散反卷积操作，此时网络中的每条边都吸收了n-1阶邻域内的结构信息。以图8中边E₁₂为例，红色虚线框区域为边E₁₂的一阶邻域，经过n阶双向扩散反卷积，边E₁₂将其自身的结构信息扩散到了其一阶邻域内的边上，也就是说，边E₁₂一阶邻域内的边吸收了其扩散的n阶结构信息。

在一个较优的实施方式中，如图9所示，判别器202可以包括卷积层和全连接层；具体地，判别器202中包括第二卷积模块213和全连接层模块214。

第二卷积模块213用于对重构的脑区-基因网络进行结构信息的聚集后得到第二关键结构子网；

全连接层模块214用于提取第二关键结构子网络的边的特征后输出该被试者为AD患者的概率的判断结果。

判别器202包括双向十字卷积层和全连接层，它的功能是判断输入的脑区-基因网络是真实的还是重构的以及是正常人的还是AD患者的。在训练结构信息聚集模型的过程时，生成器和判别器保持相互对抗的关系，判别器会对真实或重构的脑区-基因网络进行判断，根据结果不断迭代参数并持续训练直到模型收敛。此时，生成器能够提取到脑区-基因网络的关键结构子网络，判别器能够准确地对关键结构子网络进行分类。

判别器中设计双向十字卷积层用来捕获网络中边的低维表示，全连接层用来提取边的高级特征，并输出判断结果。如图9所示，判别器202中的第二卷积模块213包括n阶双向十字卷积层；判别器202中的第二卷积模块213通过n阶双向十字卷积层，对于输入的脑区-基因网络进行n次双向十字卷积操作，从而捕获脑区-基因网络中边的关键结构信息，得到第二关键结构子网络。判别器202中的第二卷积模块213输入的脑区-基因网络，既可以是真实的脑区-基因网络，也可以是生成器201重构的脑区-基因网络。

具体地，判别器202中的第二卷积模块213通过n阶双向十字卷积层，对于输入的脑区-基因网络的权重矩阵进行n次双向十字卷积操作，从而捕获脑区-基因网络中边的关键结构信息，得到第二关键结构子网络的权重矩阵。

第二卷积模块213中的n阶双向十字卷积层的输入为真实脑区-基因网络或生成器重构出的生成脑区-基因网络。第二卷积模块213中的n阶双向十字卷积层的目的是使输入的脑区-基因网络中的边能够聚集其邻域内边的关键结构信息。以第n阶双向十字卷积层为例展示卷积的过程，n阶双向十字卷积公式如公式8所示：

其中，c⁽ⁿ⁾表示第n阶双向十字卷积层的偏置；C⁽ⁿ⁾∈R^N×N表示第n阶双向十字卷积层的总核；r和c分别表示矩阵的行、列索引值；f(·)是激活函数ReLU(·)；

是

经过第n阶双向十字卷积层的双向十字卷积更新后的值。

第二卷积模块213中的n阶双向十字卷积层对输入的脑区-基因网络的权重矩阵进行n次双向十字卷积操作，与上述第一卷积模块211中的n阶双向十字卷积层相同，此处不再赘述。

全连接层模块214包括L层完全连接层，用于对第二卷积模块213输出的第二关键结构子网络提取边的特征后输出为AD患者的概率的判断结果；其中，L为大于2的自然数。

具体地，全连接层模块214对第二关键结构子网络的权重矩阵提取边的特征后输出为AD患者的概率的判断结果。

全连接层的作用是提取第二关键结构子网络中边的高级特征，并输出判断结果；如图10所示，全连接层模块214将第二卷积模块213输出的第二关键结构子网络的权重矩阵的二维数据拉直为一维数据后，对一维数据进行线性组合和非线性变换的计算，得到判断结果；其中，将二维数据拉直为一维数据的公式如公式9所示：

公式9中，Ran(·)代表按行拉直运算；拉直运算的输入为矩阵W⁽ⁿ⁾；拉直运算的输出为X⁰，即一维特征向量。

全连接层模块214进而对X⁰进行线性组合和非线性变换，最后经过激活函数Sigmoid(·)得到大小为1×2的向量，分别表示输入的脑区-基因网络为真或者为假的概率；经过Softmax(归一化)得到大小为1×2的向量，分别表示输入的脑区-基因网络为正常人或者为AD患者的概率。对一维数据进行线性组合和非线性变换的计算公式如公式10所示：

其中，W^l表示第l层全连接层的权重矩阵；X^l和X^l+1分别表示第l层和第l+1层全连接层的一维向量；bl表示第l层全连接层的偏置；Sigmoid(·)，ReLU(·)和Softmax(·)均为全连接层中的激活函数。

实施例二

在对基于生成对抗网络的结构信息聚集模型进行训练时，判别器的输入有两种：真实脑区-基因网络和生成器重构的脑区-基因网络；判别器的作用有两个，第一是判断输入的脑区-基因网络的真假性，第二是判断输入的脑区-基因网络是正常人的还是AD的。判别器和生成器互相训练，不断调优，直至生成对抗网络收敛。此时，关键结构子网络能够聚集被试脑区-基因网络的特征信息，从而根据关键结构子网络进行AD的判断，可以辅助于对阿尔茨海默症做精准的诊断。

本发明实施例二提供的一种训练上述基于生成对抗网络的结构信息聚集模型的具体方法，流程如图11所示，包括如下步骤：

步骤S1101：生成训练样本。

具体地，以正常人的脑影像(fMRI)数据和基因数据分别构建脑网络和基因网络，进而将构建的脑网络和基因网络根据现有方法进行融合，得到正常人的脑区-基因网络；具体地，以脑区和基因为节点，采用现有的皮尔逊相关分析方法计算节点之间边的权重，若边权重大于阈值则保留该边，反之去掉该边，以此构建脑区-基因网络，实现脑网络和基因网络的融合。

以AD患者的脑影像(fMRI)数据和基因数据分别构建脑网络和基因网络，进而将构建的脑网络和基因网络进行融合，得到AD患者的脑区-基因网络.

本步骤中，以多个正常人的脑影像(fMRI)数据和基因数据，分别构建多个LMCI脑区-基因网络作为样本；以及以多个AD患者的脑影像(fMRI)数据和基因数据，分别构建多个AD脑区-基因网络作为样本。

步骤S1102：根据所述样本对所述生成器和判别器进行多次迭代的优化和训练；

具体地，在模型训练过程中，判别器和生成器始终维持着对立与抗争关系，前者期盼重构的脑区-基因网络非常接近真实的脑区-基因网络，以此瞒过后者，而后者期盼能够精准判断生成的脑区-基因网络和真实的脑区-基因网络；由此，本步骤中对所述生成器和判别器进行多次迭代的优化和训练，使得二者不断进行对抗博弈，生成能力和判断能力分别得到增强，直到模型收敛，最终使得生成器能够更好地完成生成任务。

其中，对所述生成器和判别器进行一次迭代的优化和训练的过程，具体流程如图12所示包括如下子步骤：

子步骤S1201：固定所述生成器的参数后，根据所述样本对所述判别器进行优化；

子步骤S1202：固定所述判别器的参数后，根据所述样本对所述生成器进行训练。

步骤S1103：在对所述生成器和判别器进行多次迭代的优化和训练后，将由所述生成器和判别器组成的生成对抗网络作为训练得到的结构信息聚集模型。

具体地，上述子步骤S1201中提到的固定所述生成器的参数后，根据所述样本对所述判别器进行优化的具体方法流程，如图13所示，包括如下步骤：

步骤S1301：固定所述生成器的参数后，将样本中的正常人的脑区-基因网络或AD患者的脑区-基因网络输入到所述生成器中，由所述生成器根据输入的脑区-基因网络重构变化的脑区-基因网络；

步骤S1302：将所述生成器重构的脑区-基因网络或者样本中的正常人或AD患者的脑区-基因网络输入判别器；

步骤S1303：根据所述判别器输出的判断结果，对所述判别器的参数进行优化。

在根据所述判别器输出的判断结果与真实标签之间的差异，对所述判别器的参数进行优化时，是通过损失函数的计算结果进行判别器的参数优化。

例如，P_i是判别器根据输入其的第i个脑区-基因网络输出的判断结果；P_i为一个1×2的向量

表示判断结果为正常人的概率，

表示判断结果为AD患者的概率。

T_i是对应于P_i的真实标签；T_i为一个2×1的向量，值取(1，0)^T时即表示第i个脑区-基因网络为真实的正常人的脑区-基因网络，值取(0，1)^T时即表示第i个脑区-基因网络为真实的AD患者的脑区-基因网络，值取(0，0)^T时即表示第i个脑区-基因网络为重构的脑区-基因网络。

根据P_i与T_i之间的差异，对所述判别器的参数进行优化。

损失函数是评估模型的性能与预期的结果接近程度，损失函数值的大小与模型预测性能成反比。

本发明定义的总损失函数为Loss由两个值组成：重构脑区-基因网络的数据和真实脑区-基因网络的数据的差异损失值Loss1，生成器与判别器的博弈损失值Loss2。定义Loss1的公式如公式11所示：

其中，P_data(x)表示输入生成器的样本分布；

表示第x个样本的脑区-基因网络权重矩阵的第i行第j列特征值；

表示根据第x个样本真实脑区-基因网络的权重矩阵经过生成器得到的重构权重矩阵第i行第j列特征值。

构成生成器与判别器对抗训练的对抗损失函数Loss2，其计算公式如公式12所示：

其中，G(·)为生成器函数，D(·)为判别器函数，E(·)为求均值，P_data(x)表示输入生成器的脑区-基因网络的样本分布，P_z(z)表示生成器重构的脑区-基因网络的样本分布，G(z)表示生成器重构的第z个脑区-基因网络。

固定生成器优化判别器的参数，收敛函数表达形式如公式13：

判别器D的训练目标是正确区分输入的脑区-基因网络是真实的还是生成器重构的，以1、0分别代表真实的脑区-基因网络(真)、重构的脑区-基因网络(假)。D(x)和D(G(z))分别表示判别器将输入的真实脑区-基因网络和重构脑区-基因网络判断为真实脑区-基因网络的概率。因此优化的目标是D(x)＝1，D(G(z))＝0，也就是让V(D，G)越来越大。

具体地，上述子步骤S1202中提到的固定所述判别器的参数后，根据所述样本对所述生成器进行训练的具体方法流程，如图14所示，包括如下步骤：

步骤S1401：在固定所述判别器的参数后，将样本中的正常人或AD患者的脑区-基因网络输入到所述生成器中，由所述生成器根据输入的脑区-基因网络重构变化的脑区-基因网络；

步骤S1402：将所述生成器重构的脑区-基因网络输入到所述判别器；

步骤S1403：根据所述判别器输出的判断结果对所述生成器的参数进行调整。

本步骤中，根据判别器的判断结果与真实标签之间的差异，对所述生成器的参数进行调整，使生成器重构尽可能真实的脑区-基因网络来欺骗判别器。

固定判别器优化生成器的参数，收敛函数表达形式如公式14所示：

生成器G的训练目标是使重构的生成脑区-基因网络最大限度地接近真实脑区-基因网络，以此能够骗过判别器，因此D(G(z))的值越接近于1越好，此时优化的目标是D(G(z))＝1，也就是让V(D，G)越来越小。

通过反复迭代上述图13、图14所示的流程过程，判别器的判断能力和生成器的重构能力不断提高，直至生成对抗网络收敛，即总损失函数Loss最小，Loss的计算公式如公式15所示：

Loss＝Loss1+Loss2 (公式15)

实施例三

基于预先训练得到的结构信息聚集模型，本发明实施例三提供的一种阿尔茨海默症的判断方法，具体流程如图15所示，包括如下步骤：

步骤S1501：根据被试者的脑影像数据和基因数据，构建该被试者的脑区-基因网络。

步骤S1502：将该脑区-基因网络输入到通过上述方法预先训练得到的结构信息聚集模型中。

步骤S1503：通过结构信息聚集模型进行结构信息的聚集并提取边的特征后输出该被试者为阿尔茨海默症AD患者的概率的判断结果。

具体地，通过结构信息聚集模型中上述的生成器和判别器，对被试者的脑区-基因网络进行结构信息的聚集并提取边的特征后输出该被试者是否为AD患者的概率的判断结果。

例如，结构信息聚集模型输出的判断结果为P，是一个1×2的向量(P₁，sP₂)，P₁表示判断结果为正常人的概率，P₂表示判断结果为AD患者的概率。

实施例四

基于预先训练得到的结构信息聚集模型，本发明实施例三提供的一种确定阿尔茨海默症的风险要素的方法，具体流程如图16所示，包括如下步骤：

步骤S1601：将多个被试者的脑区-基因网络通过结构信息聚集模型进行结构信息聚集并提取边的特征后输出每个被试者为阿尔茨海默症AD患者的概率的判断结果。

本步骤中，将多个被试者的脑区-基因网络通过结构信息聚集模型进行结构信息聚集并提取边的特征；其中，所述多个被试者包括多个正常人和多个AD患者。

针对每个被试者，结构信息聚集模型输出该被试者为AD患者的概率的判断结果。

步骤S1602：根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各节点、边的重要度。

具体地，根据结构信息聚集模型的输出结果，计算输出的判断结果与真实情况的接近程度，如公式16所示：

Probⁱ＝P_i×T_i (公式16)

其中，P_i是结构信息聚集模型根据第i个被试者的脑区-基因网络输出的判断结果；T_i是第i个被试者的真实情况，即第i个被试者的脑区-基因网络的真实标签；

P_i为一个1×2的向置

表示判断结果为正常人的概率，

表示判断结果为AD患者的概率；T_i为一个2×1的向量，值取(1，0)^T时表示为正常人，值取(0，1)^T时表示为AD患者。因此，Prob为一个一维数值，可以反映判断结果和真实情况的接近程度，也可称为被试者的预测概率分数。若判断结果P_i为(0.9，0.1)，被试者为正常人即T_i为(1，0)^T。此时，Probⁱ＝P_i×T_i＝0.9，Probⁱ表示所述结构信息聚集模型针对第i个被试者输出的判断结果和第i个被试者的真实情况的接近程度，即预测概率分数比较大；若被试者为AD患者，T_i为(0，1)^T，此时，Probⁱ＝P_i×T_i＝0.1，表示判断结果和真实情况的接近程度，即预测概率分数比较小。

基于第一或第二关键结构子网络，以及计算的判断结果和真实情况的接近程度，求得关键特征矩阵FW，如公式17所示：

其中，KWⁱ表示所述结构信息聚集模型针对第i个被试者的脑区-基因网络生成的第一关键结构子网络或第二关键结构子网络的权重矩阵；Probⁱ表示第i个被试者的判断结果和真实情况的接近程度；Probⁱ×KWⁱ表示将第i个被试者的预测概率分数加权给权重矩阵KWⁱ的每个位置上；m表示被试的个数。其中，KWⁱ具体可以是生成器得到的第一关键子网络，也可以是判别器中得到的第二关键子网络。

将关键特征矩阵FW中第r行的每一列累加，得到节点v_r的重要度分数

从而能够得到每个节点(脑区或基因)的重要度分数，重要度分数值的大小与节点对分类的贡献度成正比，节点v_r重要度分数计算如公式18所示：

其中，FW_r，k表示关键特征矩阵FW中第r行k列的值；

表示节点v_r的重要度分数。

关键特征矩阵FW中的上三角的值即为对应边的重要度分数，例如，矩阵FW中第一行第三列的值即为节点v₁和节点v₃之间的边

的重要度分数

以此类推，可以得到每条边的重要度分数，边

重要度分数计算如公式19所示：

其中，FW_r，c表示表示关键特征矩阵FW中第r行c列的值；

表示边

的重要度分数。

步骤S1603：根据节点的重要度，确定出所述脑区-基因网络中作为AD风险要素的节点。

具体地，将计算得到的节点的重要度

进行降序排序得到节点特征集合Sub_Vertex，集合中越靠前的边和点表明其对分类的贡献越大，如公式20所示：

可以将若干个重要度最大的节点确定为AD风险要素的节点；其中，AD风险要素中的脑区节点确定为病变脑区，AD风险要素中的基因节点确定为风险基因。

作为一个更优的实施方式，还可采用增量搜索法进一步筛选出最优节点特征子集：首先对所有节点的重要度按照脑区和基因分别进行降序排序，得到一个大小为N₁＝45的基因特征集和一个大小为N₂＝116的脑区特征集。然后，从脑区特征集中选取前5个特征组成第1个脑区特征子集

以1为增量，依次选择前6，7，...，N₁个特征组成多个特征子集，最终得到(N₁-4)个脑区特征子集

同理可以得到(N₂-4)个基因特征子集

接着，将(N1-4)个脑区特征子集和(N₂-4)个基因特征子集两两任意组合得到(N₁-4)×(N₂-4)个组合，用结构信息聚集模型(SG-GAN)来测试这(N₁-4)×(N₂-4)个组合的分类准确率Acc，最高准确率对应的节点特征子集即为最后筛选的结果，寻找最优脑区特征子集和基因特征子集的函数如公式21所示：

该公式表示寻找最优脑区特征子集Sub^Brain和基因特征子集Sub^Gene，使得分类准确率Acc最大；s.t.表示受限条件，i表示脑区特征子集的编号，j表示基因特征子集的编号。最优脑区特征子集中的节点即为通过模型确定的AD风险要素；最优脑区特征子集中的节点对应的脑区即为AD风险要素中的病变脑区；最优脑区特征子集中的节点对应的基因即为模型捕获的致病基因。

步骤S1604：根据边的重要度，确定出所述脑区-基因网络中作为AD风险要素的病变脑区和风险基因之间的关联关系。

具体地，将计算得到的边的重要度

进行降序排序得到边特征集合Sub_Edge，集合中越靠前的边表明其对分类的贡献越大，如公式20所示。

进而采用增量搜索法进一步筛选出最优边特征子集：

首先，选取边特征集合前10个边特征组成第1个边特征子集

以1为增量，依次选择前6，...，N×(N-1)/2个特征组成多个边特征子集，最终得到N×(N-1)/2-10个边特征子集

用结构信息聚集模型(SG-GAN)来测试这N×(N-1)/2-10个子集的分类准确率Acc，最高准确率对应的边特征子集即为最后筛选的结果，寻找最优边特征子集的函数公式如下：

该公式表示寻找最优边特征子集Sub^Edge，使得分类准确率Acc最大；s.t.表示受限条件，i表示边特征子集的编号。最优边特征子集中的边即为模型捕获的病变脑区和致病基因之间的关联关系。

以具有4个节点的关键结构子网络为例解释3个AD患者特征提取的计算过程。如图17所示，展示的是病变脑区和致病基因以及关联模式的提取过程。假设AD患者的真实标签T_i＝(0，1)^T，i＝1，...，3，假设模型的预测结果P₁＝(0.1，0.9)，P₂＝(0.2，0.8)，P₃＝(0.5，0.5)，三个患者的关键结构子网络的权重矩阵为KW¹，KW²，KW³。首先，计算第i个AD患者的预测概率分数Probⁱ＝P_i×T_i，得到每个AD患者的预测概率分数：Prob¹，Prob²，Prob³，再计算3个AD患者关键结构子网络的权重矩阵KW中每条边的平均重要度，得到关键特征矩阵

其次，根据关键特征矩阵FW，计算重要度分数

和

并进行降序排序，得到节点特征集合Sub_Vertex和边特征集合Sub_Edge。最后使用增量搜索法寻找最优节点特征子集和最优边特征子集。寻找最优节点特征子集如下：首先，对所有节点的重要度按照脑区和基因分别进行降序排序，脑区排序结果为：v₃、v₁、v₄，基因排序结果为：v₂。选取第1个脑区组成第1个脑区特征子集，前2个脑区组成第2个脑区特征子集，前3个脑区组成第3个脑区特征子集；再选取第1个基因组成第1个基因特征子集。将3个脑区特征子集和1个基因特征子集两两任意组合得到3个组合。用SG-GAN来测试这3个组合的分类准确率，假设第2个组合的准确率最高，那么该组合中第2个脑区特征子集中的节点v₃和v₁对应的脑区即为模型捕获的病变脑区，第1个基因特征子集中的节点v₂对应的基因即为模型捕获的致病基因。寻找最优边特征子集如下：选取第1个边特征组成第1个边特征子集，前2个边特征组成第2个边特征子集，...。用模型来测试这所有边特征子集的分类准确率，假设第2个边特征子集的准确率最高，那么该子集中的边

和

即为模型捕获的病变脑区和致病基因的关联模式。

实施例五

图18示出了本实施例所提供的一种更为具体的电子设备硬件结构示意图，该设备可以包括：处理器1010、存储器1020、输入/输出接口1030、通信接口1040和总线1050。其中处理器1010、存储器1020、输入/输出接口1030和通信接口1040通过总线1050实现彼此之间在设备内部的通信连接。

处理器1010可以采用通用的CPU(Central Processing Unit，中央处理器)、微处理器、应用专用集成电路(Application Specific Integrated Circuit，ASIC)、或者一个或多个集成电路等方式实现，用于执行相关程序，以实现本说明书实施例所提供的确定阿尔茨海默症的风险要素的方法。

此外，所述处理器执行所述计算机程序时还可实现本说明书实施例所提供的结构信息聚集模型的训练方法。

存储器1020可以采用ROM(Read Only Memory，只读存储器)、RAM(Random AccessMemory，随机存取存储器)、静态存储设备，动态存储设备等形式实现。存储器1020可以存储操作系统和其他应用程序，在通过软件或者固件来实现本说明书实施例所提供的技术方案时，相关的程序代码保存在存储器1020中，并由处理器1010来调用执行。

输入/输出接口1030用于连接输入/输出模块，可以与非线性接收机相连，从非线性接收机接收信息，实现信息输入及输出。输入输出/模块可以作为组件配置在设备中(图中未示出)，也可以外接于设备以提供相应功能。其中输入设备可以包括键盘、鼠标、触摸屏、麦克风、各类传感器等，输出设备可以包括显示器、扬声器、振动器、指示灯等。

通信接口1040用于连接通信模块(图中未示出)，以实现本设备与其他设备的通信交互。其中通信模块可以通过有线方式(例如USB、网线等)实现通信，也可以通过无线方式(例如移动网络、WIFI、蓝牙等)实现通信。

总线1050包括一通路，在设备的各个组件(例如处理器1010、存储器1020、输入/输出接口1030和通信接口1040)之间传输信息。

需要说明的是，尽管上述设备仅示出了处理器1010、存储器1020、输入/输出接口1030、通信接口1040以及总线1050，但是在具体实施过程中，该设备还可以包括实现正常运行所必需的其他组件。此外，本领域的技术人员可以理解的是，上述设备中也可以仅包含实现本说明书实施例方案所必需的组件，而不必包含图中所示的全部组件。

本申请实施例还提供一种计算机可读存储介质，计算机可读存储介质其上存储有计算机程序，所述计算机程序可被至少一个处理器所执行，以使所述至少一个处理器执行实现实施例中的电子设备的操作。

本实施例中，计算机可读存储介质包括闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘等。在一些实施例中，计算机可读存储介质可以是计算机设备的内部存储单元，例如该计算机设备的硬盘或内存。在另一些实施例中，计算机可读存储介质也可以是计算机设备的外部存储设备，例如该计算机设备上配备的插接式硬盘，智能存储卡(Smart Media Card，简称为SMC)，安全数字(Secure Digital，简称为SD)卡，闪存卡(Flash Card)等。当然，计算机可读存储介质还可以既包括计算机设备的内部存储单元也包括其外部存储设备。本实施例中，计算机可读存储介质通常用于存储安装于计算机设备的操作系统和各类应用软件，例如实施例中语音视频连线方法的程序代码等。此外，计算机可读存储介质还可以用于暂时地存储已经输出或者将要输出的各类数据。

本实施例的计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。

所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本公开的范围(包括权利要求)被限于这些例子；在本发明的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，步骤可以以任意顺序实现，并存在如上所述的本发明的不同方面的许多其它变化，为了简明它们没有在细节中提供。

另外，为简化说明和讨论，并且为了不会使本发明难以理解，在所提供的附图中可以示出或可以不示出与集成电路(IC)芯片和其它部件的公知的电源/接地连接。此外，可以以框图的形式示出装置，以便避免使本发明难以理解，并且这也考虑了以下事实，即关于这些框图装置的实施方式的细节是高度取决于将要实施本发明的平台的(即，这些细节应当完全处于本领域技术人员的理解范围内)。在阐述了具体细节(例如，电路)以描述本发明的示例性实施例的情况下，对本领域技术人员来说显而易见的是，可以在没有这些具体细节的情况下或者这些具体细节有变化的情况下实施本发明。因此，这些描述应被认为是说明性的而不是限制性的。

尽管已经结合了本发明的具体实施例对本发明进行了描述，但是根据前面的描述，这些实施例的很多替换、修改和变型对本领域普通技术人员来说将是显而易见的。例如，其它存储器架构(例如，动态RAM(DRAM))可以使用所讨论的实施例。

本发明的实施例旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此，凡在本发明的精神和原则之内，所做的任何省略、修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims

1.一种确定阿尔茨海默症的风险要素的电子设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时用于实现：

2.根据权利要求1所述的电子设备，其特征在于，所述处理器执行所述计算机程序还用于实现：

3.根据权利要求1或2所述的电子设备，其特征在于，所述结构信息聚集模型包括：

4.根据权利要求3所述的电子设备，其特征在于，所述生成器中包括：

5.根据权利要求3所述的电子设备，其特征在于，所述判别器中包括：

6.根据权利要求3所述的电子设备，其特征在于，所述处理器执行所述计算机程序实现所述根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各节点的重要度时，具体为：

7.根据权利要求6所述的电子设备，其特征在于，所述处理器执行所述计算机程序实现所述根据所述结构信息聚集模型输出的判断结果与真实情况的接近程度，以及所述结构信息聚集模型中第一/二关键结构子网络的权重矩阵，计算关键特征矩阵时，具体为：

根据如下公式17计算关键特征矩阵FW：

8.根据权利要求6所述的电子设备，其特征在于，所述处理器执行所述计算机程序实现所述根据输出的判断结果与真实情况的接近程度，确定AD患者的脑区-基因网络中各边的重要度时，具体包括：

9.根据权利要求1-8任一所述的电子设备，其特征在于，所述中央处理单元执行所述程序还实现：对所述所述结构信息聚集模型进行训练：

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质内存储有计算机程序，所述计算机程序可被至少一个处理器所执行，以使所述至少一个处理器执行权利要求1-9中任意一项所述的电子设备的操作。