基于多维分层漂移扩散模型的认知决策评估方法及系统与流程

1.本发明涉及一种基于多维分层漂移扩散模型的认知决策评估方法，同时也涉及相应的认知决策评估系统，属于认知心理学技术领域。


背景技术：

2.个体在衰老过程中会经历一定程度的认知衰退。轻度认知障碍（mild cognitive impairment，简写为mci）是阿尔兹海默症的一种前驱状态，介于正常衰老与痴呆的中间状态，可作为阿尔兹海默症的“预报器”。随着计算机建模技术的发展，漂移扩散模型（drift diffusion model，简写为ddm）逐渐兴起，并迅速被应用到心理学，尤其是认知决策任务中。ddm通过分离个体行为反应数据所包含的潜在决策成分，模拟和提炼人类大脑在决策任务中的神经动态过程。ddm主要应用于单个知觉范式、决策范式，建模健康人的神经决策过程。
3.为比较不同条件或情境中个体认知决策敏感性的特点，分层漂移扩散模型（hierachical drift diffusion model，简写为hddm）应运而生。hddm采用基于贝叶斯框架的分层设计，能够综合考虑不同水平上决策成分的特征。然而，现有的hddm还仅限于从个体和群组水平上建模单一认知决策任务和范式中受试者的决策过程，并未扩展到聚类多种认知决策任务的情境中，无法整合评估个体不同认知域（区分为感知觉、高级认知和社会认知三类）的敏感性阈值和总体元认知加工效率。
4.在专利号为zl 201510357365.4的中国发明专利中，公开了一种认知障碍数据处理方法，包括以下步骤：建立认知下降程度评估工具，根据认知障碍的判断标准以及认知障碍的认知特点，编制认知测试任务，采集数据，形成各个认知域的认知测量数据以及相应认知障碍程度的数据库；利用因子分析的方法对数据库中的全部认知测验任务分别进行因子分析；对所形成的6个公因子的因子得分建立认知障碍的评估与预测模型。通过该系统或方法处理的认知障碍数据能从某种程度上预测患者得认知障碍的风险。


技术实现要素：

5.本发明所要解决的首要技术问题在于提供一种基于多维分层漂移扩散模型的认知决策评估方法，用于实现对mci患者受损能力的精准评估，促进对受损认知域内决策能力的早发现、早干预和早治疗。
6.本发明所要解决的另一技术问题在于提供一种基于多维分层漂移扩散模型的认知决策评估系统。
7.为了实现上述目的，本发明采用以下的技术方案：根据本发明实施例的第一方面，提供一种基于多维分层漂移扩散模型的认知决策评估方法，包括以下步骤：对于健康受试者，在群组水平上建立多维分层漂移扩散模型，通过拟合优度最佳的模型，获得所述健康受试者在每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布，作为健康人常模；
对于每一位mci患者，在个体水平上建立多维分层漂移扩散模型，通过拟合优度最佳的模型，获得所述mci患者在每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布，作为待衡量指标；比较所述mci患者的待衡量指标在所述健康人常模中的相对位置，评估所述mci患者在各认知域和总体元认知处理效率上的受损情况。
8.其中较优地，获得所述健康受试者每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布，具体包括：利用马尔科夫链蒙特卡罗方法，采用metropolis采样算法，根据每个认知域对应的目标概率密度函数，输入多维分层漂移扩散模型中的所有参数的联合后验分布作为抽样样本的初始状态，生成马尔科夫链状态转移矩阵；基于状态转移矩阵从已知分布中采样得到样本，以得到转移序列；所述转移序列在相同的转移矩阵作用下，不断生成新的转移序列，直至转移序列不再发生变化趋于稳定状态；当所述转移序列在第n步收敛稳定时，对符合平稳分布的对应样本集进行蒙特卡罗模拟求和，得到模型参数的后验分布；其中，n为转移序列趋于稳定状态时的作用次数；其中，所述模型参数的后验分布为群组水平的参数分布，所述群组水平的参数分布包括健康受试者在群组水平上的每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布。
9.其中较优地，获得所述mci患者每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布具体包括：利用马尔科夫链蒙特卡罗方法，采用metropolis采样算法，根据每个认知域对应的目标概率密度函数，输入多维分层漂移扩散模型中的所有参数的联合后验分布作为抽样样本的初始状态，生成马尔科夫链状态转移矩阵；基于状态转移矩阵从已知分布中采样得到样本，以得到转移序列；所述转移序列在相同的转移矩阵作用下，不断生成新的转移序列，直至转移序列不再发生变化趋于稳定状态；当所述转移序列在第n步收敛稳定时，对符合平稳分布的对应样本集进行蒙特卡罗模拟求和，得到模型参数的后验分布；其中，n为转移序列趋于稳定状态时的作用次数；其中，所述模型参数的后验分布为个体水平的参数分布，所述个体水平的参数分布包括每一位mci患者在个体水平上的每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布。
10.其中较优地，比较所述mci患者的待衡量指标在所述健康人常模中的相对位置，评估所述mci患者在各认知域和总体元认知处理效率上的受损情况具体包括：对参数后验密度分布进行显著性检验，通过贝叶斯假设检验衡量得到的各个参数的健康人常模分布与每一位mci患者的后验密度分布中提取的样本百分比之间重叠的面积；根据两个分布之间的重叠面积评估mci患者在每一个认知域敏感性阈值及其总体元认知加工效率上的受损情况。
11.其中较优地，所述健康人常模分布与所述mci患者的后验密度分布之间的重叠面
积越小，均值的差距越大，则所述mci患者的认知受损越严重。
12.其中较优地，通过比较所述mci患者的每个认知域的敏感性阈值及其概率分布，在所述健康受试者的每个认知域的敏感性阈值及其概率分布中的相对位置，评估所述mci患者在每个认知域上的认知受损情况。
13.其中较优地，通过比较所述mci患者的总体元认知处理效率的敏感性阈值及其概率分布，在所述健康受试者的总体元认知处理效率的敏感性阈值及其概率分布中的相对位置，评估所述mci患者在总体元认知处理效率上的认知受损情况。
14.其中较优地，所述认知域至少包括：感知觉决策认知域、高级认知决策认知域和社会认知决策认知域。
15.其中较优地，所述多维分层漂移扩散模型基于受试者在不同认知域下的行为反应数据建立。
16.根据本发明实施例的第二方面，提供一种基于多维分层漂移扩散模型的认知决策评估系统，包括处理器和存储器，所述处理器读取所述存储器中的计算机程序，用于执行上述的认知决策评估方法。
17.本发明所提供的基于多维分层漂移扩散模型的认知决策评估方法，采用基于多维分层漂移扩散模型计算得到的各认知域敏感性阈值和总体元认知处理效率的后验分布，将mci患者和正常人常模对比，能够更为准确地评估mci患者各认知域敏感性阈值和总体元认知加工效率在健康人常模中的相对位置，获得mci患者的整体受损情况和各个子认知域的受损状况，实现对mci患者受损能力的精准评估，促进对受损认知域内决策能力的早发现、早干预和早治疗。
附图说明
18.图1为多维分层漂移扩散模型的建模过程示意图；图2为模拟受试者在二元认知决策过程中信念累积的轨迹示意图；图3为多维分层漂移扩散模型的分层结构示意图；图4为本发明实施例提供的基于多维分层漂移扩散模型的认知决策评估方法的流程图；图5为本发明实施例提供的基于多维分层漂移扩散模型的认知决策评估系统的结构示意图。
具体实施方式
19.下面结合附图和具体实施例对本发明的技术内容进行详细具体的说明。
20.图1为本发明实施例中，多维分层漂移扩散模型的建模过程示意图。该建模过程具体包括以下步骤：s1：选择和确定需要实施的认知决策任务和范式，获取受试者在各个认知决策任务上的行为反应数据。
21.具体的，在本发明实施例中，认知决策任务至少包括：判断动点运动方向或空间位置的感知觉决策任务；基于选项相对价值做出选择的记忆、推理、执行控制的高级认知决策任务；学习他人属性、评估社会信息的社会认知决策任务。由此，根据任务内容及属性，可以
将不同的决策任务标记为感知觉决策、高级认知决策和社会认知决策这三个认知域下的任务。
22.同时，获取受试者在各个认知决策任务上的行为反应数据具体包括：获取受试者在不同认知决策任务中的反应选择及做出决策所需要的反应时间。因此，在本发明实施例中，第i个认知域下第j个任务范式中的第k个决策试次的决策时间（response time，简写为rt）和选项（decision option，简写为do）被编码进观测数据xi，j，k中，即（rt，do）～xi，j，k。i＝1，2或3。
23.s2：对受试者在不同认知决策任务上的行为反应数据进行拟合，建立多维分层漂移扩散模型。
24.具体的，包括步骤s21～s24：s21：针对受试者在单个试次中的行为反应观测数据构建偏移扩散模型。
25.该偏移扩散模型通过如下公式进行构建：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）其中，e表示一个选项优于另一个选项的累积证据平衡，提供了两个选项的偏差，即决策前受试者对某个选项的偏好，当受试者对两种决策的已有经验相同时e＝0；et表示在第t个试次（即第t个具体的任务回合）时，两个选项的累积偏差；t表示思考开始的时间，包括对信息编码与做出按键反应的时间；δ表示在每个决策时间中支持一个选项的增量证据；d表示受试者的漂移率系数，表示在每个时间单位内个体积累的关于刺激的相对信息量，表征了信息积累的速度和个体感知敏感性的测量；μi是选项i的主观价值，决策阈值d（μ1‑
μ2），代表两个决策选项之间的差异，表征决策信息收集的边界/阈值，当累积信息量达到该阈值时受试者即做出决策；σ2是扩散噪声系数，表示证据积累和/或比较系统中的噪声；n（ ） 表示正态分布，指决策轮次间的增量证据在具有特定方差（σ2）的决策阈值（μ1
‑
μ2）周围呈正态分布δ～（n）（（d（μ1
‑
μ2），σ2）。
26.s22：分析计算受试者的选择概率和预期响应的反应时间；具体的，通过以下方程式进行分析计算：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）如图2所示，该模型用以模拟受试者在二元认知决策（即从两个选项中做出决策）过程中信念累积的示例轨迹。选择概率（p）和反应时间（rt）是漂移率系数d、扩散噪声系数σ2和响应边界θ的高度的函数。其中θ是触发决策的证据积累边界的高度（此上限为对选项 1 的选择），p（ch＝1）是达到上限而非下限的概率，rt是累积过程结束的预期时间（到达任一边界）。
27.s23：计算偏移扩散模型的参数估计值，其中，偏移扩散模型的参数包括：受试者的决策阈值、相对起始点、漂移率和非决策时间；具体的，输入受试者的决策行为反应数据包括决策时间（rt）和选项（do）作为反应
时间（pt）和选择概率（p），根据公式1和公式2，计算ddm模型的参数估计值。
28.s24：将受试者在不同认知域内多个决策范式中的反应数据整合到同一模型框架中，形成由低至高的决策试次
‑
决策任务
‑
认知域
‑
总体元认知加工效率的多维分层漂移扩散模型。
29.其中，在多维分层漂移扩散模型中，低层次水平的参数估计会受到更高层次水平参数分布的约束。以决策试次
‑
决策任务这一层次举例，对于服从某种似然函数f|θ分布的每一个受试者的观测数据xi，j，k而言，假设该受试者在特定认知决策任务（感知觉、高级认知和社会认知三类）中决策轮次的参数θk在具有特定方差（σ2）的平均值（μ）周围呈正态分布（n）（λ ＝（μ，σ）），而这些任务水平参数又是从给定超先验（hyper
‑
prior）g0的数据中估计得到的，因此可以得到这样的生成性描述：μ，σ～g0（）；θk～n（μ，σ2）；xi，j，k～f（θk）。
30.对试次水平参数θ的固定先验是由任务水平参数λ参数化的随机变量，据此得到θ的后验概率公式：。
31.其中，p（θ，λ|x）为似然函数，表示在给定观测数据x的情况下试次水平参数为θ、任务水平参数为λ的可能性，p（x|θ）为给定试次水平参数θ的情况下观测数据等于x的可能性，p（θ|λ）表示任务水平参数为λ的情况下试次水平参数等于θ的可能性，p（λ）和p（x）表示任务水平参数为λ、观测数据等于x的概率。总体而言，当分层模型的层次越多，估计参数的层次越高，其后验概率估算公式越复杂。
32.如图３所示，xi，j，k
∼
l（ai， zi， vi， ti， sv， st， sz）。其中，xi， j，k表示该受试者在第i个认知域（field，简写为f）的第j个范式任务下（paradigm，简写为p）的第k个决策试次（trial，简写为t）中，包括决策时间（rt）和选项（do）在内的行为观测数据；l表示由每个认知域对应的ddm模型的似然函数；ai、zi、vi和ti分别表示第i个认知域的决策阈值、相对起始点、漂移率和非决策时间分布。其他参数μ表示均值；σ表示标准差；ak、vk、zk和tk代表在第k个决策试次中决策阈值、漂移率、相对起始点和非决策时间。
33.s3：对所述多维分层漂移扩散模型进行优化，以得到拟合优度最佳的模型。
34.通过评估拟合优度，找到收敛性达到预定值的多维分层漂移扩散模型，当拟合效果不佳时，不断调整模型参数设定从而得到合适的模型。
35.具体的，采用gelman
‑
rubin r统计量来评估当前模型的收敛性，其值越接近1，说明模型收敛性越好；计算当前模型的偏差信息准则，在考虑模型复杂度的情况下评估模型的拟合优度，其值越小说明模型拟合越好；或者，进行后验预测检验，从估计出的模型中抽取新的数据集，对比这些模拟数据集与实际观测数据的一致性，以评估当前模型对观测数据的拟合效果；当模型拟合优度不佳时，不断调整模型参数设定从而得到合适的多维分层漂移扩散模型以及相应的参数估计值和后验概率分布。
36.由此，通过上述方式可建立多维分层漂移扩散模型，该多维分层漂移扩散模型即为对mci患者的认知决策进行评估的多维分层漂移扩散模型，用于辅助了解mci患者的认知受损情况。
37.请参照图4所示，为本发明实施例提供的一种基于多维分层漂移扩散模型的认知决策评估方法，具体包括步骤s10～s30：s10：利用第一实施例得到的多维分层漂移扩散模型，建立健康人常模。
38.具体的，通过第一实施例中所述的建模方法，对于所有的健康受试者，在群组水平上建立多维分层漂移扩散模型，获得所述健康受试者在每个认知域的敏感性阈值及其概率分布，以及总体元认知处理效率的敏感性阈值及其概率分布，将此两个分布模型作为健康人常模。
39.其中，可以理解的是，所述健康受试者在每个认知域的敏感性阈值及其概率分布用于作为评估mci患者各认知域受损情况的标准；所述健康受试者在总体元认知处理效率的敏感性阈值及其概率分布用于作为评估mci患者的总体元认知处理效率受损情况的标准。
40.s20：获取每一位mci患者的待衡量指标。
41.具体的，通过第一实施例中所述的建模方法，对于每一位mci患者，在个体水平上建立多维分层漂移扩散模型，通过拟合优度最佳的模型，获得所述mci患者在每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布，作为待衡量指标。
42.其中，可以理解的是，该待衡量指标既包括mci患者在每个认知域的敏感性阈值及其概率分布，还包括mci患者在总体元认知处理效率的敏感性阈值及其概率分布，从而用于评估mci患者在各认知域及其总体认知处理效率上的受损情况。
43.s30：比较mci患者的待衡量指标在健康人常模中的相对位置，评估mci患者在各认知域和总体元认知处理效率上的受损情况。
44.其中，比较mci患者的待衡量指标在健康人常模中的相对位置具体包括步骤s301～s302：s301：计算各个参数的健康人常模的后验密度分布，以及mci患者的敏感性阈值指标的后验密度分布；s302：根据步骤301中得到的健康人常模的后验密度分布与mci患者的敏感性阈值指标的后验密度分布，提取样本百分比之间重叠的面积；对参数后验密度分布进行显著性检验，通过贝叶斯假设检验衡量得到的各个参数的健康人常模分布与每一位mci患者的后验密度分布中提取的样本百分比之间的重叠面积。
45.s303：根据重叠面积评估mci患者在各认知域的敏感性阈值及其总体元认知加工效率上的受损情况。
46.其中，健康人常模的后验密度分布与mci患者的敏感性阈值指标的后验密度分布，这两个分布之间的重叠面积越小，均值的差距越大，表示mci患者的受损越严重。
47.可以理解的是，针对各认知域而言，需要比较的是同一个认知域下，健康人常模的后验密度分布与mci患者的敏感性阈值指标的后验密度分布之间的重叠面积；针对总体元认知加工效率而言，需要比较的是两个整体分布之间的重叠面积。所述整体分布是指：根据各认知域下得到的后验密度分布，所得到的总体分布。
48.由此，通过上述方法能够实现对mci患者受损能力的精准评估，促进对受损认知域
内决策能力的早发现、早干预和早治疗。
49.在上述实施例中，对于任一受试者而言，获得每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布的方式相同，只是对于健康受试者而言，所获得的每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布作为健康人常模；对于mci患者而言，所获得的每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布作为待衡量指标。
50.下面对如何获得每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布进行说明，具体包括步骤s401～s404：s401：利用马尔科夫链蒙特卡罗方法，采用metropolis采样算法，根据每个认知域对应的目标概率密度函数，输入多维分层漂移扩散模型中的所有参数的联合后验分布作为抽样样本的初始状态，生成马尔科夫链状态转移矩阵；s402：基于状态转移矩阵从已知分布中采样得到样本，以得到转移序列；s403：所述转移序列在相同的转移矩阵作用下（作用次数为1， 2，
ꢀ…
， n），不断生成新的转移序列，直至转移序列不再发生变化趋于稳定状态；s404：当所述转移序列在第n步收敛稳定时，对符合平稳分布的对应样本集进行蒙特卡罗模拟求和，得到模型参数的后验分布；其中，n为转移序列趋于稳定状态时的作用次数。
51.其中，所述群组水平的参数分布包括健康受试者在群组水平上的每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布，所述个体水平的参数分布包括每一位mci患者在个体水平上的每个认知域的敏感性阈值及其概率分布和总体元认知处理效率的敏感性阈值及其概率分布。
52.在上述认知决策评估方法的基础上，本发明进一步提供一种基于多维分层漂移扩散模型的认知决策评估系统。如图5所示，该认知决策评估系统包括一个或多个处理器21和存储器22。其中，存储器22与处理器21耦接，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器21执行，使得所述一个或多个处理器21实现如上述实施例中的认知决策评估方法。
53.其中，处理器21用于控制该认知决策评估系统的整体操作，以完成上述认知决策评估方法的全部或部分步骤。该处理器21可以是中央处理器（cpu）、图形处理器（gpu）、现场可编程逻辑门阵列（fpga）、专用集成电路（asic）、数字信号处理（dsp）芯片等。存储器22用于存储各种类型的数据以支持在该认知决策评估系统的操作，这些数据例如可以包括用于在该认知决策评估系统上操作的任何应用程序或方法的指令，以及应用程序相关的数据。该存储器22可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器（sram）、电可擦除可编程只读存储器（eeprom）、可擦除可编程只读存储器（eprom）、可编程只读存储器（prom）、只读存储器（rom）、磁存储器、快闪存储器等。
54.在一个示例性实施例中，认知决策评估系统具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现，用于执行上述的认知决策评估方法，并达到如上述方法一致的技术效果。一种典型的实施例为计算机。具体地说，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的
任何设备的组合。
55.在另一个示例性实施例中，本发明还提供一种包括程序指令的计算机可读存储介质，该程序指令被处理器执行时实现上述任意一个实施例中的认知决策评估方法的步骤。例如，该计算机可读存储介质可以为上述包括程序指令的存储器，上述程序指令可由认知决策评估系统的处理器执行以完成上述的认知决策评估方法，并达到如上述方法一致的技术效果。
56.综上所述，本发明提供的基于多维分层漂移扩散模型的认知决策评估方法及系统，通过数学建模的方式，模拟和提炼大脑进行认知决策过程中刺激表征、信息捕获、信念积累和做出选择的动态过程，采用基于贝叶斯概率框架的分层设计，整合个体在不同认知决策任务中的反应数据，并获取个体在不同认知域（感知觉、高级认知和社会认知三类）的敏感性阈值的指标和总体元认知处理效率。
57.与现有技术相比较，本发明具有以下优点：1） 多维分层漂移扩散模型通过数学建模的手段能够精确模拟大脑的决策过程，相比于传统的采用观察记录和零假设检验分析行为反应数据的方式，多维分层漂移扩散模型能够更直接地表征和提炼大脑进行认知决策过程中刺激表征、信息捕获、信念积累和做出选择的动态神经决策过程；2） 多维分层漂移扩散模型同时整合和建模了个体在认知决策中的抉择以及做出该选择的反应时两种反应指标，综合考虑了认知任务中受试者准确性和速度的权衡过程，能够更为充分的反映受试者决策加工过程中各种认知成分相互影响和竞争的本质特征。
58.3） 多维分层漂移扩散模型在同一个数学模型内整合受试者在不同的认知域（感知觉、高级认知和社会认知三类）下的多种认知决策任务中的行为反应数据，在个体和群组两个水平上建立决策试次
‑
决策任务
‑
认知域
‑
总体元认知加工效率的分层模型，获得大脑整体加工效率，以及后验的各个认知域的敏感性阈值水平，一方面能够显著减少模型拟合的计算过程和时间，有效提升运行效率，另一方面能够在同一个能力空间内实现对受试者不同认知域敏感性阈值和总体元认知加工效率的精准定位，便于在个体水平对比同一受试者不同能力的受损情况和在群体水平上对比不同受试者相同认知决策能力的差异。
59.4） 采用基于多维分层漂移扩散模型计算得到的各认知域敏感性阈值和总体元认知处理效率的后验分布，将mci患者和正常人常模对比，能够更为准确地评估mci患者各认知域（感知觉、高级认知和社会认知三类）敏感性阈值和总体元认知加工效率在健康人常模中的相对位置，获得mci患者的整体受损情况和各个子认知域的受损状况，实现对mci患者受损能力的精准评估，促进对受损认知域内决策能力的早发现、早干预和早治疗。
60.上面对本发明所提供的基于多维分层漂移扩散模型的认知决策评估方法及系统进行了详细的说明。对本领域的一般技术人员而言，在不背离本发明实质内容的前提下对它所做的任何显而易见的改动，都将构成对本发明专利权的侵犯，将承担相应的法律责任。