《群体脑网络的建模方法及其拓扑属性分析方法.pdf》由会员分享,可在线阅读,更多相关《群体脑网络的建模方法及其拓扑属性分析方法.pdf(17页完整版)》请在专利查询网上搜索。
1、10申请公布号CN104050394A43申请公布日20140917CN104050394A21申请号201410317904722申请日20140704G06F19/12201101G06F17/30200601G06N3/0220060171申请人北京师范大学地址100875北京市海淀区新街口外大街19号72发明人朱朝喆段炼戴瑞娜74专利代理机构北京汲智翼成知识产权代理事务所普通合伙11381代理人陈曦董烨飞54发明名称群体脑网络的建模方法及其拓扑属性分析方法57摘要本发明提供了一种群体脑网络建模方法,包括下列步骤1根据群体脑成像数据定义群体脑网络的节点;2根据步骤1中定义的节点计算不同节。
2、点之间的连接矩阵作为边;3根据步骤1中定义的节点和步骤2中定义的边构建群体脑网络模型。本发明同时提供了基于上述群体脑网络建模方法实现的群体脑网络拓扑属性分析方法。通过使用上述群体脑网络建模方法及其拓扑属性分析方法,可以判断出被研究群体的交互模式特点,并进一步分析出群体中是否存在关键成员;从而可以利用群体神经活动模式结果预测群体行为。51INTCL权利要求书2页说明书13页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书13页附图1页10申请公布号CN104050394ACN104050394A1/2页21一种群体脑网络建模方法,其特征在于包括下列步骤1根据群体脑成像。
3、数据定义群体脑网络的节点;2根据步骤1中定义的节点计算不同节点之间的连接矩阵作为边;3根据步骤1中定义的节点和步骤2中定义的边构建群体脑网络模型。2如权利要求1所述的群体脑网络建模方法,其特征在于在所述步骤1中,根据脑功能整合的不同空间尺度以及所对应的社会交互过程不同,从下列三种定义方式中任选一种定义节点11局部功能脑区作为节点;一个节点为单个测量单位;12局部脑功能系统作为节点;一个节点为一个或多个感兴趣区;13整体脑功能系统作为节点;个体脑的全部测量单位作为整体当作网络中的节点。3如权利要求2所述的群体脑网络建模方法,其特征在于在所述步骤2中,采用下述三种定义方式中的一种定义群体脑网络节点。
4、之间的交互性,作为群体脑网络的边21基于群体脑网络节点的时域交互性构建连接矩阵;22基于群体脑网络节点的空间活动模式交互性构建连接矩阵;23基于群体脑网络节点的时间空间联合维度交互性构建连接矩阵。4如权利要求3所述的群体脑网络建模方法,其特征在于在所述定义方式21中,计算多种度量节点时间序列交互性的指标,选择一种交互性指标或者多种交互性指标的组合表征节点间信息的交互性,定义为群体脑网络中的边;所述度量节点时间序列交互性的指标是以下指标中的任意一种统计特性的差异、回归系数、皮尔森相关系数、相干系数、小波相干、欧几里得距离、明科夫斯基距离、互信息。5如权利要求3所述的群体脑网络建模方法,其特征在于。
5、在所述定义方式22中,计算多种度量节点空间模式交互性的指标,选择一种节点空间模式交互性指标或多种节点空间模式交互性指标的组合表征节点间信息的交互性,定义为群体脑网络中的边;所述度量群体脑网络节点的空间模式交互性的指标是以下指标中的任意一种回归系数、皮尔森相关系数、相干系数、小波相干、欧几里得距离、明科夫斯基距离、互信息。6如权利要求3所述的群体脑网络建模方法,其特征在于在所述定义方式23中,采取基于数据驱动的主成分分析和独立成分分析方法,找出与群体交互有关的群体脑网络空间模式及其动态变化过程;或者,通过整合群体脑成像数据中的时间和空间信息获得时间空间联合维度的交互性。7如权利要求3所述的群体脑。
6、网络建模方法,其特征在于当所述步骤1中,采用定义方式11定义群体脑网络的节点时,在所述步骤2中,采用定义方式21定义群体脑网络的节点;当所述步骤1中,采用定义方式12或13定义群体脑网络的节点时,在所述步骤2中,采用定义方式21、22、23中的任一种定义群体脑网络的节点。8一种群体脑网络拓扑属性分析方法,基于权利要求1所述的群体脑网络建模方法实现,其特征在于包括所述步骤1、所述步骤2和所述步骤3,权利要求书CN104050394A2/2页3此外,还包括步骤4对所述步骤3构建的群体脑网络模型分析其整体拓扑属性和/或局部拓扑属性。9如权利要求8所述的群体脑网络拓扑属性分析方法,其特征在于在所述步骤。
7、4中,使用下列一种指标或多种指标,分析群体脑网络的整体拓扑属性计算网络密度获取群体成员之间的交互紧密程度;计算全局效率和最短路径获取群体交互信息的整体传输效率;计算聚类系数获取群体交互信息的局部传输效率;计算模块化指标判断是否存在子群体;计算中心势判断群体交互过程是否呈现集中趋势。10如权利要求8或9所述的群体脑网络拓扑属性分析方法,其特征在于在所述步骤4中,使用下列一种指标或多种指标,分析群体脑网络的局部拓扑属性计算节点的度衡量一个成员与其他成员交互的多少;计算节点的最短路径长度和效率衡量一个成员与其他成员交互的效率;计算节点的中心度衡量某一群体成员是否是群体交互中的关键人物;计算节点的聚类。
8、系数,判断所有与该节点代表的成员存在直接交互的成员间彼此交互的紧密程度。权利要求书CN104050394A1/13页4群体脑网络的建模方法及其拓扑属性分析方法技术领域0001本发明涉及一种群体脑网络建模方法,尤其涉及一种基于图论的群体脑网络建模方法,同时涉及一种基于上述群体脑网络建模方法实现的群体脑网络拓扑属性分析方法。背景技术0002人类大多以群体作为基本生活方式。小至家庭、单位,大至民族、国家,人们总是以各种方式组织在一起,产生了形形色色的群体。群体社会行为通常远较个体社会行为复杂和丰富,并常常体现出一些独特的群体心理现象,例如群体凝聚力、群体智慧、群体极化等等。0003近年来,多人交互同。
9、步记录HYPERSCANNING技术的出现使同时观测交互过程中的群体脑活动变为现实,为获取群体心理行为的神经活动数据提供了全新的脑成像技术。在获得群体脑成像数据之后,如果能够实现群体脑成像数据的建模,那么就可以对群体神经活动进行解码,从而利用群体神经活动模式预测群体行为。例如,通过对团队神经活动模式的建模,解码出不同团队凝聚力程度的神经活动模式,从而预测某一团队是否具有高凝聚力。但是,目前群体脑成像数据的建模尚为空白。0004现有技术中的个体脑区间建模及分析方法,具体可参见公开号为CN1626031A的中国发明申请中提供的“基于图论的脑区间功能连接的检测方法”。这种检测方法基于磁共振成像技术,。
10、在图论的框架下,通过对人脑不同脑区之间的关系进行数学建模,来检测脑区间的功能连接强度,进而评价某一脑区在其所属的功能网络中所起的作用。这种脑区间功能连接的检测方法可以用于临床医学中脑疾病的诊断,病人用药前后脑区连接差异的评估,脑区功能网络的研究,对脑区发育变化的认识以及对人脑工作机理的理解等。0005群体脑网络与个体脑网络存在根本不同。从信息处理的角度来看,群体中的每个个体其交互行为必须依赖群体中其他成员的信息。在交互过程中,每个个体不断地接收群体提供的信息输入并做出处理,再将反应结果信息输出到群体环境中为他人所获取。群体社会交互即是所有参与交互的个体同时地不断地输入、处理、输出交互信息的过程。
11、。换言之,群体社会交互过程可以看作是每个个体的信息不断地在其他群体成员脑中表达的连续过程。因此,群体交互行为与普通的个体认知行为的一个根本区别在于,群体交互行为的结果取决于各个群体成员脑中信息表达之间的关系,而不是仅由各个群体成员脑中信息表达本身所决定。这种群体成员脑间信息表达的关系称为“脑间连接”。0006上述群体成员间的脑间连接可以使用网络模型来描述。将群体中每名成员的大脑看作网络中的节点,将群体成员间的脑间连接作为网络的边。其中,某些节点间存在直接相连的边,则这些成员的大脑彼此间可以通过脑间连接进行直接的信息映射。某些节点间不存在直接相连的边,而需要借助经过其他节点的通路连通,则这些成员。
12、大脑彼此之间的信息需要通过中间媒介进行传输。整个群体的脑间连接网络即群体脑网络的拓扑属性可能是影响群体行为的关键所在。而群体脑网络的动态变化则会反映在群体交互的动态变化上。某一个体的神经活动模式可以通过群体脑网络进行广泛的传播,从而导致行为或心说明书CN104050394A2/13页5理的传播。当群体脑网络具有足够高的信息传递效率时,群体成员间就可以产生彼此间神经活动信息的快速映射,从而使群体体现出某些行为表现如群体成员行为保持高度一致等。0007因此,群体脑网络和个体脑网络具有以下差别第一,这两类脑网络表征的心理行为过程不同个体脑网络表征的是个体心理过程,如个体的感知觉、注意状态、情绪状态等。
13、,群体脑网络表征的是群体成员之间进行交互的过程。第二,网络建模方法及复杂性不同个体脑网络中的节点是个体脑中的不同脑区,边是脑区之间的结构或功能连接。群体脑网络中的节点是每个个体的脑,边是个体脑之间的脑间连接,因此,群体脑网络建模具有个体脑网络建模所不具备的复杂性。第三,在群体脑网络中,每个个体脑可以作为一个整体当做节点,产生较为简单的平面网络;每个个体脑也可以分为多个脑区提供多个节点,相当于将网络的节点从全部个体扩展到全部个体的全部脑区,此时产生的群体脑网络将体现出立体的层次结构,每一层包含所有个体相同脑区的节点;层内的连接表示不同个体相同脑区之间的神经活动同步性,层间的连接表示不同脑区之间的。
14、神经活动同步性。0008此外,目前已有研究者利用脑电图EEG,ELECTROENCEPHALOGRAM同时观测两个人在社会决策中的神经活动,构建两脑的网络,并分析该网络拓扑属性和两人决策行为之间的关系具体参见DEVICOFALLANIF,NICOSIAV,SINATRAR,ASTOLL,CINCOTTIF,ETAL2010DEFECTINGORNOTDEFECTINGHOWTO“READ”HUMANBEHAVIORDURINGCOOPERATIVEGAMESBYEEGMEASUREMENTSPLOSONE512E14187DOI101371/JOURNALPONE0014187。但是,需要注意。
15、的是,对群体脑成像数据的建模与两人脑成像数据的建模也具有极大的不同。首先,群体脑网络和两人脑网络的意义,即其所表示的两人交互过程与群体交互过程存在明显不同。正如PERC等人2013年对结构化群体的演化动力学研究综述中所指出的“结构化群体的群体交互要远远地多于群体中两两配对交互的总和。无论对于规则网络还是复杂网络,其交互策略的复杂性、不同的公共利益函数和演化过程都为使用统计物理方法进行研究提供了丰富的内容具体参见PERCM,GOMEZGARDENESJ,SZOLNOKIA,FLORIALM,MORENOY2013EVOLUTIONARYDYNAMICSOFGROUPINTERACTIONSONS。
16、TRUCTUREDPOPULATIONSAREVIEWJRSOCINTERFACE1020120997。”例如,在群体交互过程中,虽然每个群体成员往往无法同时与每一个其他群体成员都进行有效的沟通,但整个群体最终却能够达到一定的同步性。这说明群体成员间的信息传递并不全都是基于“点对点”PEERTOPEER的传输,也需要借助由其他群体成员介导的“信息通路”。因此,群体成员间脑间连接的构建,可能同时存在“点对点”的直接连接和“借助第三方介导”的间接连接。其次,群体脑网络具有两人脑网络所没有的复杂性。如,群体脑网络中脑间连接即网络的边的可能数量随着群体人数的增加成倍增长;群体脑网络可能具有树形、星形等。
17、多样的结构,表征群体交互中不同的信息传递模式;群体脑网络也可以具有模块化的特点,表征群体在交互过程中产生的子群体;群体脑网络的局部属性能够用来刻画群体成员多样的角色如领导者、不合群者等。0009综上可知,由于群体脑网络与个体脑网络和两人脑网络存在根本的区别,因此,无法直接将现有建模方法借用到群体脑网络的建模中,为此,需要研究出新的适用于群体脑网络的建模方法。说明书CN104050394A3/13页6发明内容0010本发明所要解决的首要技术问题在于提供一种群体脑网络建模方法。0011本发明所要解决的另一技术问题在于提供一种基于上述群体脑网络建模方法实现的群体脑网络拓扑属性分析方法。0012为了实。
18、现上述的发明目的,本发明采用下述的技术方案0013一种群体脑网络建模方法,包括下列步骤00141根据群体脑成像数据定义群体脑网络的节点;00152根据步骤1中定义的节点计算不同节点之间的连接矩阵作为边;00163根据步骤1中定义的节点和步骤2中定义的边构建群体脑网络模型。0017其中较优地,在所述步骤1中,根据脑功能整合的不同空间尺度以及所对应的社会交互过程不同,从下列三种定义方式中任选一种定义节点001811局部功能脑区作为节点;一个节点为单个测量单位;001912局部脑功能系统作为节点;一个节点为一个或多个感兴趣区;002013整体脑功能系统作为节点;个体脑的全部测量单位作为整体当作网络中。
19、的节点。0021其中较优地,在所述步骤2中,采用下述三种定义方式中的一种定义群体脑网络节点之间的交互性,作为群体脑网络的边002221基于群体脑网络节点的时域交互性构建连接矩阵;002322基于群体脑网络节点的空间活动模式交互性构建连接矩阵;002423基于群体脑网络节点的时间空间联合维度交互性构建连接矩阵。0025其中较优地,在所述定义方式21中,计算多种度量节点时间序列交互性的指标,选择一种交互性指标或者多种交互性指标的组合表征节点间信息的交互性,定义为群体脑网络中的边;0026所述度量节点时间序列交互性的指标是以下指标中的任意一种统计特性的差异、回归系数、皮尔森相关系数、相干系数、小波相。
20、干、欧几里得距离、明科夫斯基距离、互信息。0027其中较优地,在所述定义方式22中,计算多种度量节点空间模式交互性的指标,选择一种节点空间模式交互性指标或多种节点空间模式交互性指标的组合表征节点间信息的交互性,定义为群体脑网络中的边;0028所述度量群体脑网络节点的空间模式交互性的指标是以下指标中的任意一种回归系数、皮尔森相关系数、相干系数、小波相干、欧几里得距离、明科夫斯基距离、互信息。0029其中较优地,在所述定义方式23中,采取基于数据驱动的主成分分析和独立成分分析方法,找出与群体交互有关的群体脑网络空间模式及其动态变化过程;或者,通过整合群体脑成像数据中的时间和空间信息获得时间空间联合。
21、维度的交互性。0030其中较优地,当所述步骤1中,采用定义方式11定义群体脑网络的节点时,在所述步骤2中,采用定义方式21定义群体脑网络的节点;0031当所述步骤1中,采用定义方式12或13定义群体脑网络的节点时,在所述步骤2中,采用定义方式21、22、23中的任一种定义群体脑网络的节点。0032一种群体脑网络拓扑属性分析方法,基于上述群体脑网络建模方法实现,包括上说明书CN104050394A4/13页7述步骤1、步骤2和步骤3,此外,还包括步骤4,对所述步骤3构建的群体脑网络模型分析其整体拓扑属性和/或局部拓扑属性。0033其中较优地,在所述步骤4中,使用下列一种指标或多种指标,分析群体脑。
22、网络的整体拓扑属性0034计算网络密度获取群体成员之间的交互紧密程度;0035计算全局效率和最短路径获取群体交互信息的整体传输效率;0036计算聚类系数获取群体交互信息的局部传输效率;0037计算模块化指标判断是否存在子群体;0038计算中心势判断群体交互过程是否呈现集中趋势。0039其中较优地,在所述步骤4中,使用下列一种指标或多种指标,分析群体脑网络的局部拓扑属性0040计算节点的度衡量一个成员与其他成员交互的多少;0041计算节点的最短路径长度和效率衡量一个成员与其他成员交互的效率;0042计算节点的中心度衡量某一群体成员是否是群体交互中的关键人物;0043计算节点的聚类系数,判断所有与。
23、该节点代表的成员存在直接交互的成员间彼此交互的紧密程度。0044本发明提供的群体脑网络建模方法,通过选择群体脑成像数据定义节点,并通过计算不同节点之间的连接矩阵定义边;从而构建群体脑网络模型。本发明提供的群体脑网络拓扑属性分析方法,基于上述群体脑网络建模方法实现。通过使用上述群体脑网络建模方法及其拓扑属性分析方法对群体脑成像数据进行处理,可以判断出群体的交互模式特点,并进一步分析出群体中是否存在关键成员;从而可以利用群体神经活动模式预测群体行为,并对群体行为进行针对性的指导。附图说明0045图1是本发明所提供的群体脑网络建模方法及其拓扑属性分析方法的逻辑框图。具体实施方式0046下面结合附图和。
24、具体实施例对本发明的技术内容进行详细具体的说明。0047本发明提供的群体脑网络拓扑属性分析方法,包括基于群体脑成像数据的群体脑网络建模过程,还包括在构建群体脑网络模型之后,对群体脑网络的拓扑属性进行分析的过程。在群体脑网络构建过程中,首先需要选择与社会认知相关的脑区作为节点,根据脑功能分化理论和脑功能整合理论,可以从不同的空间尺度选择节点。其次需要计算节点之间的关系作为网络的边,根据群体脑活动具有的时间变化规律和空间活动模式,可以从时间维度、空间维度以及更为复杂的时间空间联合维度定义群体脑网络节点之间的交互性,作为群体脑网络的边。在构建群体脑网络之后,通过从群体脑网络的整体层面和局部层面对群体。
25、脑网络进行分析,可以分析出群体的神经信息传输模式和效率,并可以获得群体中重要个体的神经信息传输模式。0048下面结合图1所示的逻辑框图对本群体脑网络建模方法及其拓扑属性分析方法进行介绍。本发明提供的群体脑网络建模方法包括下列步骤1根据群体脑成像数据定说明书CN104050394A5/13页8义群体脑网络的节点;2根据步骤1中定义的节点计算不同节点之间的连接矩阵作为边,从而构建群体脑网络;3根据步骤1中定义的节点和步骤2中定义的边构建群体脑网络模型。0049下面对群体脑网络建模过程中的详细步骤进行介绍。0050根据心理理论和镜像神经元理论的现有研究,人类大脑中负责社会认知的专门区域,被称为“社会。
26、脑”,包括前额叶、颞顶联合区、颞极、顶下小叶、颞上沟后部等区域,它们与感知、理解和推断他人意图有密切关系,对人类能够完成正常的社交活动具有重要的意义。因此,在群体脑网络构建过程中,我们可以选择单个个体的“社会脑”作为网络节点。在步骤1中,定义节点的方式有三种。根据脑功能整合的不同空间尺度以及所对应的社会交互过程不同,从下列三种定义方式中任选一种定义节点0051A局部功能脑区作为节点。一个局部功能脑区负责社会认知功能中的一个子功能即社会交互中的一个具体环节,如理解他人意图,该定义适用于对群体成员间交互过程的一个特定环节进行建模。此时,一个节点为单个测量单位,包括FMRI中,一个体素;脑电中,一个。
27、通道;NIRS中,一个测量导;0052B局部脑功能系统作为节点,局部脑功能系统反映的是相邻若干局部功能脑区的功能整合,共同表征了特定种类的交互过程如面对面交谈。该定义适用于对特定的一种交互过程进行建模。此时,一个节点为一个或多个感兴趣区ROI,REGIONOFINTEREST一个ROI包括FMRI中,多个体素组成的一块脑区、多个脑区;脑电中,多个通道;NIRS中,多个测量导;0053C整体脑功能系统作为节点。整体脑功能系统表征了人类全部的社会认知功能,这种情况适用于将交互过程作为一个整体进行建模。此时,即个体脑的全部测量单位作为整体当作网络中的节点,包括FMRI中,全部体素;脑电中,全部通道;。
28、NIRS中,全部测量导。0054在步骤1中确定群体脑网络节点之后,需要在步骤2中对每两个节点之间关系进行定义,从而构建群体成员神经活动的连接矩阵。根据群体脑活动具有的时间变化规律和空间活动模式,可以从时间维度、空间维度以及更为复杂的时间空间联合维度定义群体脑网络节点之间的交互性,作为群体脑网络的边。下面分别对21基于群体脑网络节点的时域交互性构建连接矩阵;22基于群体脑网络节点的空间活动模式交互性构建连接矩阵;23基于群体脑网络节点的时间空间联合维度交互性构建连接矩阵进行具体说明。0055根据群体脑网络中节点的不同定义方式,群体中每名成员可以为群体脑网络提供一个或多个节点。设节点A、B为群体脑。
29、网络中的两个节点,以下都将以A、B两个节点为例说明如何构建群体脑网络中的边。00561基于群体脑网络节点的时域交互性构建连接矩阵0057由于社会交互是一个连续的动态过程,所以可以从时间维度构建群体脑网络的边。该定义适用于节点定义方式A、B和C。0058群体脑成像会获取一段时间内的群体交互过程中所有群体成员的脑活动,在所观测的时间尺度内,可以对群体脑网络节点的时间序列之间的关系进行度量,构建基于时域交互性的连接矩阵。0059具体来说对群体进行一段时间的群体脑成像,群体脑网络中的每个节点都存说明书CN104050394A6/13页9在一条时间序列。节点A与节点B的时间序列可表示为N维向量X和Y,分。
30、别为XX1,XN,YY1,YN,N表示数据采样点数。通过计算多种度量节点时间序列交互性的指标,选择一种交互性指标或者多种交互性指标的组合表征节点间信息的交互性,定义为群体脑网络中的边。0060度量节点时间序列交互性的指标可以是以下任意一种统计特性的差异、回归系数、皮尔森相关系数、相干系数、小波相干、欧几里得距离、明科夫斯基距离、互信息。指标计算方法见下文00612基于群体脑网络节点的空间活动模式交互性构建连接矩阵0062对于预测某一时刻T1的交互状态如群体将要做出一个怎样的决定,需要利用T1之前的群体脑活动信息。越临近T1的信息,对于T1的交互状态预测力越强,但同时也削减了时间维度所包含的信息。
31、。因此,可以利用T1之前某一个或某几个时刻的空间信息对群体脑网络进行建模,从而实现预测某一时刻的交互状态。但值得注意的是,只有节点定义方式B、C能够提供某一时刻的空间信息。0063即将一个或多个感兴趣区或全部测量单元作为群体脑网络中的节点,则每个节点都包括多个测量单元。每次采样,节点的多个测量单元都各自有一个采样值,共同表示了节点在这个采样时刻的空间活动模式SPATIALPATTERN。设群体脑网络中的两个节点A、B,在某个采样时刻T,节点A、B的空间活动模式分别表示为XX1,XM、YY1,YM,M表示节点所包含的测量单元数目。通过计算多种度量节点空间模式交互性的指标,并选择一种或多种节点空间。
32、模式交互性指标的组合表征节点间信息的交互性,定义为群体脑网络中的边。0064度量群体脑网络节点的空间模式交互性的指标可以是以下指标中的任意一种回归系数、皮尔森相关系数、相干系数、小波相干、欧几里得距离、明科夫斯基距离、互信息。0065上述指标的具体计算方法如下0066A统计特性的差异0067与表示两个节点时间序列的均值,X与Y表示两个节点时间序列的标准差。节点时间序列统计特性差异有均值的差比值标准差的差值XY、比值X/Y。节点A与节点B的边可以是有向的,节点A到节点B的边表示为均值的差或标准差的差XY;节点B到节点A的边表示为均值的差或标准差的差YX。节点A与节点B的边也可以是无向的,可以定义。
33、为均值差的绝对值标准差的差的绝对值|XY|、均值的比值标准差的比值X/Y。0068B回归系数0069节点间时间序列间的、空间活动模式间的线性关系表示为0070YABX。其中A为回归常数,B为回归系数。0071C皮尔森相关系数说明书CN104050394A7/13页100072节点间时间序列间的、空间活动模式间的皮尔森相关系数表示为其中与为两个节点时间序列的均值。0073D相干系数0074节点间时间序列间的、空间活动模式间的相干系数表示为其中CX,Y为X与Y的互功率谱,PX,X与PY,Y分别为X与Y的自功率谱。0075E小波相干0076节点间时间序列间的、空间活动模式间的小波相干系数表示为007。
34、7其中SWXXT,F、SWXYT,F为在频率F和时刻T时,信号的光滑小波谱和交叉小波谱。0078F欧几里得距离EUCLIDEANDISTANCE0079节点间时间序列间的、空间活动模式间的欧几里得距离表示为0080G明科夫斯基距离MINKOWSKIDISTANCE0081节点间时间序列间的、空间活动模式间的明氏距离表示为这里的P值是一个变量。0082当P1时,数学上为L1范数,表示的是曼哈顿距离MANHATTANDISTANCE,是将多个维度上的距离进行求和后的结果0083当P2时,数学上为L2范数,表示的即是欧几里得距离。0084当P趋于无穷时,数学上为无穷范数,表示的是切比雪夫距离0085。
35、H互信息0086节点间时间序列间的、空间活动模式间的互信息表示为0087IX;YHXHYHX,Y,其中HX与HY为离散序列的信息熵,HX,Y为联合熵。0088群体脑网络中连接节点A与节点B的边可以定义为一种节点A,B时间序列间的交互性指标;或者,是在采样时刻T,群体脑网络中连接节点A与节点B的边可以定义为一种说明书CN104050394A108/13页11节点A,B空间活动模式的交互性指标。上述指标可表示为LINKTIMEA,BP,P可以是统计特性的差异仅限时域交互性、回归系数、皮尔森相关系数、相干系数、小波相干、欧几里得距离、明科夫斯基距离、互信息中的任意一种;此外,也可以定义为节点A,B时。
36、间序列多个相似性指标的组合LINKTIMEA,BFTIMEP1,PN,P1,PN代表节点A,B时间序列多个交互性指标,FTIME代表组合方式可以是平均,也可以自定义不同指标的权值系数;或是节点A,B空间模式多个相似性指标的组合LINKSPATIALA,BFSPATIALP1,PN,P1,PN代表节点A,B空间模式多个交互性指标,FSPATIAL代表组合方式可以是平均,也可以自定义不同指标的权值系数。00893基于群体脑网络节点的时间空间联合维度交互性构建连接矩阵。0090对于后两种节点定义方式B、C,由于其包含了群体脑活动在某一时刻的空间信息,因此,当使用这两种节点定义方式时,可以对交互的全过。
37、程及其每一时刻的交互状态同时进行建模。即从时间空间联合维度进行边的定义。0091在每个采样时刻,每一个节点所包含的多个测量单元均有一个采样值,共同构成该时刻的空间活动模式;在一段采样时间内,这些空间活动模式以一定的时间序列进行动态变化。因此,可以采取基于数据驱动的主成分分析和独立成分分析方法,找出与群体交互有关的群体脑网络空间模式及其动态变化过程;此外,还可以通过整合群体脑成像数据中的时间和空间信息来获得时间空间联合维度的交互性。0092A基于数据驱动的主成分分析和独立成分分析0093对一个群体脑网络,有P个节点,采样点数为N。群体脑网络数据可记为00940095其中J1,2,P为第J个节点数。
38、据。0096对于这样一个群体脑网络,主成分分析可以有效地找到其中最重要的群体脑活动模式。主成分分析是将P个观测变量综合成为P个新的变量,即FAX,其中,00970098F1为第一主成分,F2为第二主成分,依此类推,有P个主成分;A称为主成分系数矩阵。每个主成分对应群体脑网络中一种时间空间活动模式。0099ICA将群体脑网络中每个节点信号按照统计独立原则,分解为若干独立成分。P个说明书CN104050394A119/13页12节点记录的信号看成是来自P个信号源的信号SS1,S2,SP按照加权方式A线性组合而成0100XAS,ICA分解是上式逆过程是源信号的估计,解混矩阵W是混合矩阵A的估计,01。
39、01代表一系列空间模式,对应群体脑网络中的分离出的独立的不同空间模式。代表这些空间模式的时间活动,可以分析小群体空间、时间特性确定群体中不同成员的角色等。0102B基于时间窗口的时间空间联合维度交互性算法0103如果将一个或多个感兴趣区或全部测量单元作为群体脑网络中的节点,在时间窗口T内进行群体脑成像,可以在时间尺度定义群体脑网络的边,进而构建一个群体脑网络;对于时间窗口T内的每个采样时刻,也可以在空间尺度构建群体脑网络,假设时间窗口T包含N个采样点,则时间窗口T内包含N个空间维度的群体脑网络。0104在时间空间联合尺度,群体脑网络中的边可定义为0105LINKTIMESPATIALA,BFT。
40、IMESPATIALLINKTIMEA,B,LINKSPATIALNA,B0106其中LINKTIMEA,B为在时间窗口T,时间尺度上连接节点A与节点B的边;LINKSPATIALNA,B为在时间窗口T,空间尺度上连接节点A与节点B的边。在时间窗口T内,空间尺度连接节点A与节点B的边定义为表示时间窗口T内N个采样时刻节点A、B边的均值。FTIMESPATIAL为组合函数。0107上面对群体脑网络建模过程中节点和边的定义方式进行了详细介绍,在具体使用时,可以根据实际测量的社会交互过程不同,任意选择节点的定义方式,并根据节点定义方式计算两个节点间的连接矩阵作为边。具体来说,当步骤1中采用定义方式A。
41、定义群体脑网络的节点时,在所述步骤2中,只能采用定义方式21定义群体脑网络的节点;当所述步骤1中采用定义方式B或C定义群体脑网络的节点时,在所述步骤2中,可以采用定义方式21、22、23中的任一种定义群体脑网络的节点。0108此外,在步骤3中,在使用上述方法得到连接矩阵之后,即可构建无权的或加权的群体脑网络。下文介绍群体脑网络指标的数学定义时,用AAIJ和WWIJ分别代表无权网络和加权网络的连接矩阵。对一些网络指标,无权网络和加权网络的定义不同,用标号B和W分别代表无权网络和加权网络。如无特别说明,下文所有公式都基于GN,K,网络G包括N个节点和K条边。0109下面对基于上述群体脑网络建模方法。
42、实现的群体脑网络拓扑属性分析方法进行介绍。群体脑网络拓扑属性分析方法,包括上述步骤1、步骤2和步骤3的脑体脑网络建模过程,还包括步骤4对步骤3构建的群体脑网络模型分析其整体拓扑属性和/或局部拓扑属性。下面对步骤4中的分析过程进行介绍。0110如前所述,群体社会交互过程可以看作是每个个体的信息不断地在其他群体成员脑中表达的连续过程。这一过程的特点可以用群体脑网络的拓扑属性来表达。整个群体的交互模式特点由群体脑网络的整体拓扑属性来表示;群体交互中是否存在关键成员等信息说明书CN104050394A1210/13页13则由群体脑网络的局部拓扑属性来表示。01111群体脑网络的整体拓扑属性分析0112。
43、在群体交互过程中,整个群体的交互特点可以用该群体脑网络的整体拓扑属性来刻画,例如所有群体成员之间的交互紧密程度可以群体脑网络的网络密度来表示;群体交互信息的整体传输效率可以由群体脑网络的全局效率和最短路径长度来表示,即群体中每个成员最少可以通过几个人与其他任意成员进行信息交换;群体交互信息的局部传输效率可以由群体的聚类系数来表示,即群体中任意一个成员和他临近成员进行信息交互的效率;群体中是否存在子群体SUBGROUP可以由群体脑网络的模块化指标来表示,即子群体成员内部信息交换频繁,而与子群体之外的成员联系很少;群体交互过程是否呈现集中趋势可以由群体脑网络的中心势来表示。0113群体脑网络的整体。
44、拓扑属性指标计算如下0114A网络密度0115网络密度为或ES为网络可能的最大边数,0116B网络的聚类系数与标准聚类系数0117网络所有节点聚类系数的平均值为网络聚类系数0118标准聚类系数为,其中为N个节点构成的完全随机网络的聚类系数。0119C全局效率0120全局网络效率定义为DIJ为节点I与J之间最短路径长度。0121D最短路径长度与标准最短路径长度0122网络最短路径长度为任意两个节点间最短路径长度的平均值,DIJ为点I与J之间最短路径长度。0123标准最短路径长度为,其中为N个节点构成的完全随机网络的特征路径长度。0124E模块化0125对一个给定的网络分割度P,模块化指数Q定义为。
45、说明书CN104050394A1311/13页140126或NM为模块数量,LW为网络总连接数总权重,LSWS为模块S中节点间连接的总边数总权重,DSWS为模块S中节点的节点度和值。0127F网络的中心势0128绝对度数中心势定义为0129相对度数中心势定义为0130中间中心势定义为0131接近中心势定义为01322群体脑网络的局部拓扑属性分析0133在群体交互过程中,往往存在一些特殊的群体成员,比如领导者、不合群者、子群体等,可以用群体脑网络的局部拓扑属性定义出这些特殊成员,并分析他们与其他群体成员的交互模式。群体脑网络中节点的度可以衡量一个成员与其他成员交互的多少,节点的最短路径长度和效率。
46、可以衡量一个成员与其他成员交互的效率,如果一个成员的神经活动与他人的交互很少或者效率很低,可能该成员不愿与其他成员交互,是群体中的不合群者;节点的中心度可以衡量某一群体成员是否是群体交互中的关键人物如领导者;节点的聚类系数用来判断所有与该节点代表的成员存在直接交互的成员间彼此交互的紧密程度。0134群体脑网络的局部拓扑属性指标计算如下0135A节点度0136节点I的节点度定义为或AIJWIJ是AW二值加权网络邻接矩阵中第I,J个元素。0137B节点的最短路径长度0138节点I与J的最短路径长度记为DIJ,对二值网络,两个节点的最短路径长度为连接两个节点的最小边数;对加权网络,两个节点的最短路径。
47、长度为两个节点间最短的距离,加权网络中相邻节点距离为边的倒数例如节点I和J距离为0139C节点的效率说明书CN104050394A1412/13页150140节点I的效率定义为DIJ为节点I和J的最短路径长度。0141D节点的中心度0142节点的度数中心度节点I的绝对度数中心度与节点I直接相连的节点数目,CADI。节点I的相对度数中心度节点的绝对度数中心度与节点最大可能的节点度之比,表示为0143节点的中间中心度节点I的绝对中间中心度其中JK是从节点I到节点K的所有最短路径的数量,JKI是这些最短路径通过节点I的数量。节点的相对中间中心度其中CMAX为节点可能的最大中间中心度只有星型网络下可能。
48、,CMAXN23N2/2。0144节点的接近中心度节点I的绝对接近中心度是该点与图中所有其他节点的最短路径长度之和,记为0145节点I的相对接近中心度只有在星型网中才可能取最小值N1。0146E节点的聚类系数0147节点I聚类系数定义为或E为节点I邻居节点间连接的边数。0148上面对基于群体脑网络建模方法的群体脑网络拓扑属性分析方法进行了介绍。基于上述分析方法,通常可以一并分析群体脑网络的整体拓扑属性和局部拓扑属性。此外,在对群体脑网络的拓扑属性进行分析时,当对群体整体交互效率及模式不太关注的情况下,可以根据构建的群体脑网络模型直接计算其局部拓扑属性,对群体中的个体的社会性格进行分析,而该结论。
49、的获得并不依赖于整体拓扑属性的计算。基于上述分析结果,可以判断出群体的交互模式特点,并进一步分析出群体中是否存在关键成员;从而可以利用群体神经活动模式预测群体行为,并对群体行为进行针对性的指导。例如,判断群体凝聚力,并对群体成员进行针对性的指导,从而提高群体凝聚力。0149下面,结合一个实例对该群体脑网络建模方法及其拓扑属性分析方法的具体应用进行举例说明。这是一个对9人群体敲鼓交互过程中群体脑网络进行建模以及对其群体拓扑属性进行分析的实验。说明书CN104050394A1513/13页160150受试群体的实验任务是进行不同交互程度的群体敲鼓任务,分别为条件1随意敲鼓低交互、条件2自发协同敲鼓高交互和条件3跟随节拍器敲鼓低交互。实验流程包括3个敲鼓任务条件以及任务之间的休息,其中每次敲鼓250秒,中间休息约90秒。实验使用ETG4000设备观测4名受试者,使用LABNIRS设备观测5名受试者。每名受试者覆盖了1个ROI区域共4个观测通道,该ROI覆盖于背侧前额叶皮层DPFC,该脑区是典型的社会认知相关区域,被认为与推断他人意图有关。01511群体脑网络中节点的定义0152选择一个感兴趣区ROI,REGIONOFINTEREST背侧前额叶皮层DPFC作为群体脑网络的节点。群体中每名成员为群体脑。