一种基于弥散张量成像的抑郁症用药决策模型的构建方法技术

本发明专利技术公开了一种基于弥散张量成像的抑郁症用药决策模型的构建方法。该方法包括：构建全脑结构连接网络SN；构建分类器特征集合；构建双层stacking集成模型；药物选择标记物识别。该方法以客观的影像学指标来辅助抑郁症患者的临床用药决策，此过程没有人工参与判断，避免了主观因素带来的误差，模型找到不同用药被试的特异性影像学标记，为临床抗抑郁药物决策提供帮助。策提供帮助。策提供帮助。

【技术实现步骤摘要】
一种基于弥散张量成像的抑郁症用药决策模型的构建方法


[0001]本专利技术涉及用药模型构建方法，具体涉及一种基于弥散张量成像的抑郁症用药决策模型的构建方法。

技术介绍

[0002]抑郁症是一种异质性很强的精神类疾病，个体抑郁症患者在治疗反应上存在着显著的差别，治疗效果的异质性往往使得临床用药决策更为复杂。目前，精神科医师对抑郁症患者进行抗抑郁药物治疗大多依据患者的临床表现和主诉症状等，具有很强的主观性，对临床医生的问询经验具有极高的要求。因此，即使应用了最佳实践指南，抑郁症治疗的反应率和缓解率仍然很低。大多数患者都需经过几个疗程或者“反复试错”才能达到缓解，然而反复的临床试药不仅增加了医疗系统的负担，同时也增加了患者的治疗成本和治疗时间。
[0003]此外，目前的一些基线模型大多为纯数据驱动，在一定程度上忽略了模型的生理学可解释性。且多数仅仅预测用药4周或者8周的治疗反应，未考虑过抑郁症患者早期治疗疗效。研究抑郁症治疗早期的大脑恢复机制，将有利于调整治疗手段、减少医疗资源浪费的作用。

技术实现思路

[0004]专利技术目的：本专利技术的模型以客观的影像学指标来辅助抑郁症患者临床用药决策，解决目前临床抑郁症用药的不足。
[0005]技术方案：本专利技术提供的一种基于弥散张量成像的抑郁症用药决策模型的构建方法，包括如下步骤：
[0006](1)构建全脑结构连接网络SN(Structural Network)：利用自动解剖标记模板和弥散张量图像，计算每个建模样本的全脑结构连接矩阵；<br/>[0007](2)构建分类器特征集合：按照建模样本的用药情况以及患者的汉密尔顿量表17项在抗抑郁药物治疗2周后的减分率，将整个建模样本集分为SSRI改善组和SNRI改善组；提取先验脑区的对外结构连接矩阵，构造基分类器特征集合，并进行基分类器的训练；
[0008](3)构建双层stacking集成模型：依据每个先验脑区对外结构连接矩阵，构建第一层基分类器，并用留一法LOOCV(Leave one out cross
‑
validation)预测得到一个最终的输出；采用stacking集成方法，将第一层基分类器的输出作为第二层输入特征集合，从而构建得到双层集成分类器模型。
[0009](4)药物选择标记物识别：在基模型中，依据敏感度曲线和特异度曲线分别得到最优的特征个数，依据敏感度最优特征和特异度最优特征分别重新构建基模型，并进行集成。在每次集成模型训练的过程中依次剔除一个基模型，计算相对于原始敏感性/特异性最优集成模型的敏感性和特异性变化值，变化最大的即为药物选择标记物。
[0010]进一步地，步骤(1)所述大脑结构网络构建过程如下：使用自动解剖标记模板(AAL结构模板)将建模样本中每个样本被试的大脑分成90个脑区，每个脑区作为脑结构网络的
节点。对建模样本中每个被试的弥散张量数据进行纤维束示踪，两脑区对之间的纤维束连接作为脑结构网络的边，计算弥散各向异性指数FA(Fractional Anisotropy)为边的权重。
[0011]进一步地，步骤(2)所述构建基分类器特征集合，包括用药标签集合和影像特征集合两种，方法如下：
[0012](1)用药标签：SSRI改善组和SNRI改善组为患者经过SSRI或SNRI进行2周治疗后HAMD量表17项总分的减分率达到20％，其中，经过单一SSRI用药的患者为SSRI改善组，经过单一SNRI用药的患者记为SNRI改善组。其中SSRI类药物包括：舍曲林，帕罗西汀，氟伏沙明，氟西汀和艾司西酞普兰，SNRI类药物包括文拉法辛和度洛西汀。
[0013](2)影像特征集合：先验脑区是通过统计以往研究中对抗抑郁药物SSRI和SNRI具有显著治疗反应的脑区。具体包括背外侧额上回、眶部额上回、额中回、三角部额下回、内侧额上回、脑岛、前扣带和旁扣带脑回、后扣带回、海马、海马旁回、杏仁核、舌回、枕上回、枕中回、枕下回、梭状回、顶上回、尾状核、豆状壳核、豆状苍白球、丘脑、颞中回。
[0014]对于脑区m，被试a在脑区m的对外连接向量为外连接向量为被试者b在脑区m的对外连接向量为被试者b在脑区m的对外连接向量为将所有N个被试在脑区m的对外连接矩阵进行堆叠，得到对应脑区m的基分类器特征合集M：
[0015][0016]进一步地，步骤(3)所述构建双层stacking集成模型方法如下：基模型训练采用线性支持向量机SVM(support vector machine)进行训练，线性核函数公式为：
[0017]K(x,x
i
)＝x
·
x
i
[0018]特征筛选方法采用支持向量机递归特征消除算法SVM
‑
RFE(support vector machine recursive feature elimination)，该算法是一个基于SVM的最大间隔原理的序列后向选择算法。该算法通过SVM训练样本，然后对每个特征得分进行排序，去掉最小特征得分的特征，然后用剩余的特征再次训练模型，进行下一次迭代，最后选出需要的特征数。SVM
‑
RFE算法原理如下：
[0019]输入：训练样本
[0020]输出：特征排序集R
[0021]1)初始化原始特征集合S＝{1,2,
…
,D}，特征排序集R＝[][0022]2)循环以下过程直至S＝[][0023]获取带候选特征集合的训练样本；
[0024]用训练SVM分类器，得到ω；
[0025]用k＝1,2,
…
,|S|计算排序准则得分；
[0026]找出排序得分最小的特征p＝argmin
k
c
k
；
[0027]更新特征集R＝[p,R]；
[0028]在S中去除此特征：S＝S/p
[0029]基模型的性能采用留一法进行评估，该验证方法在每一轮随机保留一个样本作为测试数据，其他剩余样本作为训练数据进行模型训练，并训练得到预测标签，该预测标签作为第二层的输入特征合集。
[0030]该模型为双层stacking模型，第一层的预测标签作为第二层的输入，依旧采用SVM
‑
RFE进行特征筛选，同样采用线性支持向量机进行模型的训练，最后得到一条不同特征数量对应的准确率曲线。
[0031]模型评价指标包括准确度、特异度、灵敏度三个指标。准确度计算公式：模型评价指标包括准确度、特异度、灵敏度三个指标。准确度计算公式：特异度计算公式：灵敏度计算公式：灵敏度计算公式：通过这三个指标，我们可以对模型性能进行有效评估。
[0032]进一步地，步骤(4)所述药物选择标记物识别过程如下：
[0033]在每一个基分类器中，采用SVM
‑
RFE进行特征筛选，得到敏感度和特异度曲线，依据这两条曲线分别得到基分类器的最优特征个数，并找到对应的最优特征集合，依据本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于弥散张量成像的抑郁症用药决策模型的构建方法，其特征在于，包括如下步骤：(1)构建全脑结构连接网络SN：利用自动解剖标记模板和弥散张量图像，计算每个建模样本的全脑结构连接矩阵；(2)构建分类器特征集合：按照建模样本的用药情况以及患者的汉密尔顿量表17项在抗抑郁药物治疗2周后的减分率，将整个建模样本集分为SSRI改善组和SNRI改善组；提取先验脑区的对外结构连接矩阵，构造基分类器特征集合，并进行基分类器的训练；(3)构建双层stacking集成模型：依据每个先验脑区对外结构连接矩阵，构建第一层基分类器，并用留一法LOOCV预测得到一个最终的输出；采用stacking集成方法，将第一层基分类器的输出作为第二层输入特征集合，从而构建得到双层集成分类器模型；(4)药物选择标记物识别：在基模型中，依据敏感度曲线和特异度曲线分别得到最优的特征个数，依据敏感度最优特征和特异度最优特征分别重新构建基模型，并进行集成，在每次集成模型训练的过程中依次剔除一个基模型，计算相对于原始敏感性/特异性最优集成模型的敏感性和特异性变化值，变化最大的即为药物选择标记物。2.根据权利要求1所述的基于弥散张量成像的抑郁症用药决策模型的构建方法，其特征在于：所述步骤(1)所述全脑结构连接网络SN构建过程如下：使用自动解剖标记模板将建模样本中每个样本被试的大脑分成90个脑区，每个脑区作为脑结构网络的节点，对建模样本中每个被试的弥散张量数据进行纤维束示踪，两脑区对之间的纤维束连接作为脑结构网络的边，计算弥散各向异性指数FA为边的权重。3.根据权利要求1所述的基于弥散张量成像的抑郁症用药决策模型的构建方法，其特征在于：所述步骤(2)所述构建基分类器特征集合，包括用药标签集合和影像特征集合两种，方法如下：用药标签集合：SSRI改善组和SNRI改善组为患者经过SSRI或SNRI进行2周治疗后HAMD量表17项总分的减分率达到20％，其中，经过单一SSRI用药的患者为SSRI改善组，经过单一SNRI用药的患者记为SNRI改善组；影像特征集合：先验脑区是通过统计以往研究中对抗抑郁药物SSRI和SNRI具有显著治疗反应的脑区，具体包括背外侧额上回、眶部额上回、额中回、三角部额下回、内侧额上回、脑岛、前扣带和旁扣带脑回、后扣带回、海马、海马旁回、杏仁核、舌回、枕上回、枕中回、枕下回、梭状回、顶上回、尾状核、豆状壳核、豆状苍白球、丘脑、颞中回，对于脑区m，被试a在脑区m的对外连接向量为对于脑区m，被试a在脑区m的对外连接向量为被试者b在脑区m的对外连接向量为被试者b在脑区m的对外连接向量为将所有N个被试在脑区m的对外连接矩阵进行堆叠，得到对应脑区m的基分类器特征合集M：4.根据权利要求1所述的基...

【专利技术属性】
技术研发人员：卢青，姚志剑，薛莉，
申请(专利权)人：东南大学，
类型：发明
国别省市：