一种风险预测模型的构建方法、设备、终端及存储介质与流程

[0001]
本发明涉及风险预测模型技术领域，特别涉及一种风险预测模型的构建方法、设备、终端及存储介质。


背景技术：

[0002]
糖尿病合并冠心病，区别于普通冠心病患者，冠脉病变弥漫而且复杂，经常容易合并糖尿病自主神经病变，很多病人心脏缺血症状较为隐匿，不容易早期识别及干预，导致心血管死亡风险高，社会危害极大。随着本世纪冠心病治疗技术的发展，糖尿病合并冠心病患者的预后得到了显著改善，在一定情况下减少了由于心血管导致死亡的风险，已有的治疗方案包括：标准的药物治疗、在药物治疗基础上的冠脉介入治疗和冠脉搭桥等，但由于糖尿病合并冠心病涉及到的病变情况复杂，目前糖尿病合并冠心病患者的死亡率仍然很高。
[0003]
目前的治疗策略可以改善糖尿病合并冠心病患者的长期预后，但仍然有部分患者生存周期有限。在此情况下，如果能准确预测糖尿病合并冠心病患者心血管死亡风险，以及对死亡风险影响因素贡献大小作评估，有助于准确识别高危的糖尿病合并冠心病患者，做到及时有效干预。
[0004]
但现有的针对糖尿病合并冠心病患者的风险因素的研究集中于评估单个危险因素对冠心病变患者预后的影响，缺少多个危险因素综合评估的研究，且不能对相关危险因素贡献度大小综合评估，没有一个有效的心血管死亡预测机制。


技术实现要素：

[0005]
有鉴于此，本发明提出了一种风险预测模型的构建方法、设备、终端及存储介质，通过本方案建立了针对糖尿病合并冠心病人的心血管死亡风险预测模型，有助于建立个体化干预方案，降低冠冠心病患者的死亡风险，提高严重冠心病患者生存率。
[0006]
具体的，本发明提出了以下具体的实施例：
[0007]
本发明提出了一种风险预测模型的构建方法，包括：
[0008]
步骤s1、获取糖尿病合并冠心病患者的多个属性数据；
[0009]
步骤s2、基于信息增益评价函数对各所述属性数据进行评估，得到各所述属性数据的信息增益，并基于多个所述属性数据组成数据集；
[0010]
步骤s3、初始化迭代次数n，n＝1；
[0011]
步骤s4、第n次基于所述数据集在深度神经网络中进行训练生成模型；
[0012]
步骤s5、基于损失函数对所述模型进行性能评估，得到性能值；
[0013]
步骤s6、若所述性能值与初始模型的性能值的差值不大于性能阈值，从所述数据集剔除所述信息增益最小的一个或多个所述属性数据，并将所述迭代次数加1后作为n，返回步骤s4；所述初始模型为n＝1时基于所有所述属性数据形成的所述数据集在深度神经网络中进行训练生成的；
[0014]
步骤s7、若所述性能值与所述初始模型的性能值的差值大于性能阈值，将第n-1次
生成的模型作为用于预测糖尿病合并冠心病人的心血管死亡风险的最终模型。
[0015]
在一个具体的实施例中，所述属性数据关联有病例样本的数据；所述步骤s2，包括：
[0016]
基于信息增益评价函数与所述病例样本的数据确定第一信息熵，且基于信息增益评价函数、所述属性数据及所述病例样本的数据确定第二信息熵；所述第一信息熵为所述病例样本中心血管死亡事件的信息熵；所述第二信息熵为不同的所述属性数据取值下的心血管死亡事件的信息熵；
[0017]
基于所述第一信息熵与所述第二信息熵的差值确定各所述属性数据的信息增益。
[0018]
在一个具体的实施例中，所述信息增益评价函数为：
[0019]
g(d,a)＝h(d)-h(d|a)；
[0020]
其中，g(d,a)为所述信息增益；
[0021]
h(d)为第一信息熵；p
i
包括p0、p1；；n0和n1分别为所述属性数据对应的糖尿病合并冠心病患者中没有发生心血管死亡事件和发生心血管死亡事件的病例数；
[0022]
a为所述属性数据，a为所述属性数据的值；p(a＝a)为所述属性数据的值为a时的概率；h(d|a＝a)为当所述属性数据的值是a时的第二信息熵。
[0023]
在一个具体的实施例中，所述步骤s6，包括：
[0024]
若所述性能值与初始模型的性能值的差值不大于性能阈值，基于预设的剔除阈值从所述数据集剔除所述信息增益最小的一个或多个所述属性数据；
[0025]
其中，剔除的一个所述属性数据的所述信息增益小于所述剔除阈值，或剔除的多个所述属性数据的所述信息增益的和小于所述剔除阈值。
[0026]
在一个具体的实施例中，在将第n-1次生成的模型作为所述最终模型之后，该方法还包括：
[0027]
将第n-1次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第一属性数据；
[0028]
将第n-2次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第二属性数据；
[0029]
基于所述第一属性数据与所述第二属性数据两者的信息增益确定针对所述属性数据的信息增益阈值。
[0030]
在一个具体的实施例中，所述步骤s1，包括：
[0031]
获取糖尿病合并冠心病患者的病例样本；其中，所述病例样本包括没有发生心血管死亡的第一病例样本和发生心血管死亡的第二病例样本；所述第一病例样本的数量大于所述第二病例样本的数量；
[0032]
对所述第一病例样本进行随机下采样，以使所述第一病例样本的数量和所述第二病例样本的数量相同；
[0033]
在所述第一病例样本与第二病例样本中获取属性数据。
[0034]
在一个具体的实施例中，所述步骤s4，包括：
[0035]
基于所述属性数据确定输入值，并基于所述属性数据关联的病例样本的数据确定输出值；
[0036]
基于所述输入值与输出值确定所述深度神经网络中每个连接的权值以及每个神经元的偏差值；
[0037]
基于所述权值与所述偏差值建立针对糖尿病合并冠心病人的心血管死亡风险预测模型。
[0038]
本发明实施例还提出了一种风险预测模型的构建设备，包括：
[0039]
获取模块，用于获取糖尿病合并冠心病患者的多个属性数据；
[0040]
确定模块，用于基于信息增益评价函数对各所述属性数据进行评估，得到各所述属性数据的信息增益，并基于多个所述属性数据组成数据集；
[0041]
初始化模块，用于初始化迭代次数n，n＝1；
[0042]
建立模块，用于第n次基于所述数据集在深度神经网络中进行训练生成模型；
[0043]
评估模块，用于基于损失函数对所述模型进行性能评估，得到性能值；
[0044]
迭代模块，用于若所述性能值与初始模型的性能值的差值不大于性能阈值，从所述数据集剔除所述信息增益最小的一个或多个所述属性数据，并将所述迭代次数加1后作为n，启动所述建立模块；所述初始模型为n＝1时基于所有所述属性数据形成的所述数据集在深度神经网络中进行训练生成的；
[0045]
处理模块，用于若所述性能值与所述初始模型的性能值的差值大于性能阈值，将第(n-1)
[0046]
次生成的模型作为用于预测糖尿病合并冠心病人的心血管死亡风险的最终模型。
[0047]
在一个具体的实施例中，所述属性数据关联有病例样本的数据；所述确定模块，用于：
[0048]
基于信息增益评价函数与所述病例样本的数据确定第一信息熵，且基于信息增益评价函数、所述属性数据及所述病例样本的数据确定第二信息熵；所述第一信息熵为所述病例样本中心血管死亡事件的信息熵；所述第二信息熵为不同的所述属性数据取值下的心血管死亡事件的信息熵；
[0049]
基于所述第一信息熵与所述第二信息熵的差值确定各所述属性数据的信息增益。
[0050]
在一个具体的实施例中，所述信息增益评价函数为：
[0051]
g(d,a)＝h(d)-h(d|a)；
[0052]
其中，g(d,a)为所述信息增益；
[0053]
h(d)为第一信息熵；p
i
包括p0、p1；；n0和n1分别为所述属性数据对应的糖尿病合并冠心病患者中没有发生心血管死亡事件和发生心血管死亡事件的病例数；
[0054]
a为所述属性数据，a为所述属性数据的值；p(a＝a)为所述属性数据的值为a时的概率；h(d|a＝a)为当所述属性数据的值是a时的第二信息熵。
[0055]
在一个具体的实施例中，所述迭代模块，用于：
[0056]
若所述性能值与初始模型的性能值的差值不大于性能阈值，基于预设的剔除阈值从所述数据集剔除所述信息增益最小的一个或多个所述属性数据；
[0057]
其中，剔除的一个所述属性数据的所述信息增益小于所述剔除阈值，或剔除的多个所述属性数据的所述信息增益的和小于所述剔除阈值。
[0058]
在一个具体的实施例中，还包括：信息增益阈值模块，用于在将第(n-1)次生成的模型作为所述最终模型之后，将第(n-1)次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第一属性数据；
[0059]
将第(n-2)次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第二属性数据；
[0060]
基于所述第一属性数据与所述第二属性数据两者的信息增益确定针对所述属性数据的信息增益阈值。
[0061]
本发明实施例还提出了一种终端，包括用于执行上述方法的处理器。
[0062]
相较于现有技术，本方案具有以下技术效果：
[0063]
通过本方案建立了针对糖尿病合并冠心病人的心血管死亡风险预测模型，且该心血管死亡风险预测模型实现了对危险因素综合评估，以及对相关危险因素贡献度大小进行评估；因此本方案中的心血管死亡风险预测模型有助于建立个体化干预方案，降低冠冠心病患者的死亡风险，提高严重冠心病患者生存率。
附图说明
[0064]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0065]
图1为本发明实施例提出的一种风险预测模型的构建方法的流程示意图；
[0066]
图2为本发明实施例提出的一种风险预测模型的构建设备的结构示意图；
[0067]
图3为本发明实施例提出的一种终端的结构示意图。
具体实施方式
[0068]
在下文中，将更全面地描述本公开的各种实施例。本公开可具有各种实施例，并且可在其中做出调整和改变。然而，应理解：不存在将本公开的各种实施例限于在此公开的特定实施例的意图，而是应将本公开理解为涵盖落入本公开的各种实施例的精神和范围内的所有调整、等同物和/或可选方案。
[0069]
在本公开的各种实施例中使用的术语仅用于描述特定实施例的目的并且并非意在限制本公开的各种实施例。如在此所使用，单数形式意在也包括复数形式，除非上下文清楚地另有指示。除非另有限定，否则在这里使用的所有术语(包括技术术语和科学术语)具有与本公开的各种实施例所属领域普通技术人员通常理解的含义相同的含义。所述术语(诸如在一般使用的词典中限定的术语)将被解释为具有与在相关技术领域中的语境含义相同的含义并且将不被解释为具有理想化的含义或过于正式的含义，除非在本公开的各种
实施例中被清楚地限定。
[0070]
实施例1
[0071]
本发明实施例1公开了一种风险预测模型的构建方法，如图1以及图2所示，包括以下步骤：
[0072]
步骤s1、获取糖尿病合并冠心病患者的多个属性数据；
[0073]
在一个具体的实施例中，由于在实际获取到的数据中，没有发生心血管死亡的病例远多于发生心血管死亡的病例，由此步骤s1包括：获取糖尿病合并冠心病患者的病例样本；其中，所述病例样本包括没有发生心血管死亡的第一病例样本和发生心血管死亡的第二病例样本；所述第一病例样本的数量大于所述第二病例样本的数量；对所述第一病例样本进行随机下采样，以使所述第一病例样本的数量和所述第二病例样本的数量相同；在所述第一病例样本与第二病例样本中获取属性数据。
[0074]
步骤s2、基于信息增益评价函数对各所述属性数据进行评估，得到各所述属性数据的信息增益，并基于多个所述属性数据组成数据集；
[0075]
步骤s3、初始化迭代次数n，n＝1；
[0076]
步骤s4、第n次基于所述数据集在深度神经网络中进行训练生成模型；具体的，当n＝1时，数据集中包括所有属性数据，生成的是初始模型；后续不断进行s4-s6的迭代，直到满足条件执行s7确定最终模型，即为需要建立的模型；具体的n≥1且n为正整数，用于标识执行s4的次数；
[0077]
步骤s5、基于损失函数对所述模型进行性能评估，得到性能值；
[0078]
步骤s6、若所述性能值与初始模型的性能值的差值不大于性能阈值，从所述数据集剔除所述信息增益最小的一个或多个所述属性数据，并将所述迭代次数加1后作为n，返回步骤s4；所述初始模型为n＝1时基于所有所述属性数据形成的所述数据集在深度神经网络中进行训练生成的；
[0079]
步骤s7、若所述性能值与所述初始模型的性能值的差值大于性能阈值，将第n-1次生成的模型作为用于预测糖尿病合并冠心病人的心血管死亡风险的最终模型。
[0080]
具体的，通过对糖尿病合并冠心病患者进行多项指标的检查(在机器学习中，这些指标被称为属性)，得到患者的属性数据，属性数据对应的属性可以有多种，例如性别、年龄、身高(cm)、身高(cm)、体重(kg)、入院日期、入组日期即造影检查日期、治疗策略、高血压病史、陈旧心梗史、既往心梗史类型、既往心梗次数、既往血运重建史、溶栓史等等多种，具体的每个属性对应有属性的值，每个值均可通过数字来进行标识，即使是非数值类型的值也通过数字来进行标识，例如性别为男通过1来进行标识，性别为女则通过2来进行标识；在一个实施例中，具体的属性数据中属性的名称及对应的值可以如表1所示：
[0081]
表1
[0082]
[0083]
[0084][0085]
[0086]
具体的，糖尿病合并冠心病患者的属性数据有多个类型，例如可以有：(1)相关属性：对于诊断有帮助，可以提升机器学习算法的效果；(2)无关属性：对于诊断没有任何帮助，不会给算法的效果带来任何提升；(3)冗余属性：不会对诊断带来新的信息，或者这种属性的信息可以由其他的属性推断出。具体的，属性数据的类型是相关属性、无关属性还是冗余属性，是基于重要程度，也即是信息增益来进行判断的，具体的信息增益是对后续需要建立的模型带来性能提升的能力。
[0087]
此外，所述属性数据关联有病例样本的数据；所述步骤s2，包括：
[0088]
基于信息增益评价函数与所述病例样本的数据确定第一信息熵，且基于信息增益评价函数、所述属性数据及所述病例样本的数据确定第二信息熵；所述第一信息熵为所述病例样本中心血管死亡事件的信息熵；所述第二信息熵为不同的所述属性数据取值下的心血管死亡事件的信息熵；
[0089]
基于所述第一信息熵与所述第二信息熵的差值确定各所述属性数据的信息增益。
[0090]
进一步的，所述属性数据关联有病例样本的数据；所述病例样本的数据中包含所有病例样本中发生心血管死亡事件的概率、所有病例样本中没有发生心血管死亡事件的概率、各所述属性数据不同取值对应的发生心血管死亡事件的概率、各所述属性数据不同取值对应的没有发生心血管死亡事件的概率；具体的具体的，属性数据和与属性数据相关联的病例样本的数据组成了样本数据，在确定各个属性数据的信息增益时，需要用到样本数据。由此，所述步骤s2，包括：
[0091]
基于信息增益评价函数、所有病例样本中发生心血管死亡事件的概率、所有病例样本中没有发生心血管死亡事件的概率确定心血管死亡事件的信息熵；
[0092]
基于信息增益评价函数、各所述属性数据不同取值对应的发生心血管死亡事件的概率、各所述属性数据不同取值对应的没有发生心血管死亡事件的概率确定各所述属性数据不同取值下心血管死亡事件的信息熵；
[0093]
基于所述心血管死亡事件的信息熵与各所述属性数据不同取值下心血管死亡事件的信息熵的差值确定各所述属性数据的信息增益。
[0094]
其中，具体步骤s2的所述信息增益评价函数为：
[0095]
g(d,a)＝h(d)-h(d|a)；
[0096]
其中，g(d,a)为所述信息增益；
[0097]
h(d)为第一信息熵，也即心血管死亡事件的信息熵；p
i
包括p0、p1；n0和n1分别为所述属性数据中没有发生心血管死亡和发生心血管死亡的病例数；
[0098][0099]
a为所述属性，a为所述属性的值；p(a＝a)为所述属性的值为a时的概率；h(d|a＝a)为当所述属性的值是a时，基于没有发生心血管死亡事件的概率和发生心血管死亡事件的概率得到的信息熵，也即第二信息熵。
[0100]
具体的，以属性treatment strategy为例来进行对如何确定信息增益进行说明，
属性treatment strategy的重要性计算对应的心血管死亡事件的信息熵为：
[0101]
h(d)＝-(0.5*log(0.5)+0.5*log(0.5))＝1；
[0102]
具体的,基于步骤s1的具体流程，第一病例样本的数量和第二病例样本的数量相同，因为心血管死亡事件为0和1的样本相同，所以心血管死亡事件为0和1的概率为0.5。
[0103]
如表1所示，属性treatment strategy有三种取值(1，2，3)；
[0104]
1.当属性treatment strategy为1时，心血管死亡事件的信息熵为：
[0105]
h(d|treatment strategy＝1)＝-(0.6449*log(0.6449)+0.3551*log(0.3551))＝0.9385；其中treatment strategy为1时心血管死亡事件为0的概率为0.6449，treatment strategy为1时心血管死亡事件为1的概率为0.3551。
[0106]
2.当属性treatment strategy为2时，心血管死亡事件的信息熵为：
[0107]
h(d|treatment strategy＝2)＝-(0.7119*log(0.7119)+0.2881*log(0.2881))＝0.8663；其中treatment strategy为2时心血管死亡事件为0的概率为0.7119，treatment strategy为2时心血管死亡事件为1的概率为0.2881。
[0108]
3.当属性treatment strategy为3时，心血管死亡事件的信息熵为：
[0109]
h(d|treatment strategy＝3)＝-(0.2381*log(0.2381)+0.7619*log(0.7619))＝0.7919；其中treatment strategy为3时心血管死亡事件为0的概率为0.2381，treatment strategy为3时心血管死亡事件为1的概率为0.7619。
[0110]
综上所述，引入属性treatment strategy后，属性的三个值对应的心血管死亡事件的信息熵，也即第二信息熵为：
[0111]
h(d|treatment strategy)＝(0.4125*0.9385+0.1987*0.8663+.3889*0.7919)＝0.8671；其中，treatment strategy为1的概率为0.4125，treatment strategy为2的概率为0.1987，treatment strategy为3的概率为0.3889。
[0112]
因此属性treatment strategy的信息增益为：
[0113]
g(d,treatment strategy)＝h(d)-h(d|treatment strategy)＝1-0.8671
[0114]
也即g(d,treatment strategy)＝0.1329。
[0115]
参照计算上述属性treatment strategy的信息增益的方法，计算表1中所有临床变量的信息增益，由于信息增益为0代表无法帮助模型进行性能提升，因此出于减少数据量的目的，可以先将信息增益为0的变量删除，得到表2；
[0116]
表2
[0117]
[0118]
[0119][0120]
具体的，如表2所述的属性数据可以组成数据集，基于表2中的所有属性数据与深度神经网络建立的模型即为初始模型；在此情况下，基于步骤s6，不断进行迭代，每次迭代可以逐步删除信息增益最小的一个或多个属性数据，并将剩余的属性数据组成新的数据集来建立模型，并与初始模型进行比较，基于合适的性能阈值确定最终模型，具体的模型的性能阈值可以为0，或者根据性能要求选择其他数值作为性能阈值；若性能阈值为0代表删除了属性数据之后，训练得到的最终模型的性能相对初始模型不会降低。
[0121]
例如所有属性数据的数量为20个，其中最小的属性数据为属性数据1，在此情况下，当n＝1时，初始模型是基于包括20个属性数据的数据集在深度神经网络中训练得到的，这种情况下，由于训练生成的模型只有初始模型，性能值的差值为0，肯定不会大于性能阈值，需要剔除数据集中的属性数据1，对数据集进行更新，然后返回步骤s3，此时迭代次数加1，也即n＝n+1，得到n＝2，此时模型是基于包括19个属性数据的数据集在深度神经网络中训练得到的，得到的该模型与初始模型进行性能值得比较，以此迭代，直到满足s7的条件。
[0122]
例如当基于表3中对应的所有属性数据得到的模型为最终模型。
[0123]
表3
[0124]
[0125][0126]
进一步的，所述步骤s6，包括：
[0127]
若所述性能值与初始模型的性能值的差值不大于性能阈值，基于预设的剔除阈值从所述数据集剔除所述信息增益最小的一个或多个所述属性数据；
[0128]
其中，剔除的一个所述属性数据的所述信息增益小于所述剔除阈值，或剔除的多个所述属性数据的所述信息增益的和小于所述剔除阈值。
[0129]
具体的，剔除数据集中属性数据时，是基于剔除阈值来选择被剔除的属性数据的，例如当剔除阈值为0.0001，则每次剔除属性数据的信息增益的和需要小于但是靠近0.0001这个值，通过这种方式可以实现快速迭代，减少迭代的次数，提高得到最终模型的效率，具体的剔除阈值可根据经验进行设置与调整。
[0130]
此外，考虑到后续可能得到新的属性数据，在此情况下，需要根据新的属性数据来训练得到新的最终模型，且希望这个过程尽可能快，不需要经过本方案的多次迭代过程，为此，在将第n-1次生成的模型作为所述最终模型之后，该方法还包括：
[0131]
将第n-1次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第一属性数据；
[0132]
将第n-2次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第二属性数据；
[0133]
基于所述第一属性数据与所述第二属性数据两者的信息增益确定针对所述属性数据的信息增益阈值。
[0134]
例如，如表2与表3所示，所述最终模型，也即第n-1次生成的模型对应的所述数据集中信息增益最小的所述属性数据为“tni”、第n-2次生成的模型对应的所述数据集中信息增益最小的所述属性数据为“慢性肾功能不全病史”；
[0135]
在此情况下，可以设置信息增益阈值的值为0.009833261825842321-0.010435200421651913之间，例如可以设置信息增益阈值为0.01。
[0136]
根据本实施例中上述实际数据，例如可以设置信息增益阈值为0.01，即只保留信息增益大于0.01的属性数据作为训练模型的数据。如此，不需要经过多次迭代，可以基于该信息增益阈值直接从新的属性数据中筛选出训练模型用的数据，
[0137]
仍以上述表1，表2的数据为例，基于该信息增益阈值的设置，如表3所示，输入模型的属性从表2中的81个减少到表3中19个，经过对模型的性能进行评估，可知模型性能有提升，具体的，可以通过二值预测对模型的性能进行评估涉及到的参数如表4所示：
[0138]
表4二值预报结果的四种可能结果
[0139] 预测的正类预测的负类真实的正类正确的肯定(tp)错误的否定(fn)真实的负类错误的肯定(fp)正确的否定(tn)
[0140]
具体的，性能评估可以为二值预测问题，正类表示死亡，负类表示存活。对于一个二值预测系统，其预测结果存在四种可能性，如表4所示。本身是正类又被正确地预测为正类的样本被称为正确的肯定，本身是负类又被正确地预测为负类的样本被称为正确的否定，本身是正类又被错误地预测为负类的样本被称为错误的否定，本身是负类又被错误地预测为正类的样本被称为错误的肯定。
[0141]
基于表4所示的四类输出，定义如下两个指标刻画心血管死亡风险预测模型的性能：
[0142]
tprate＝n
tp
/(n
tp
+n
fn
)；其中，n
tp
是正确的肯定样本数，n
fn
是错误的否定样本数。tnrate＝n
tn
/(n
tn
+n
fp
)；n
tn
是正确的否定样本数，n
fp
是错误的肯定样本数。由此，tp rate和tn rate分别用于评价死亡预测的准确程度和存活预测的准确程度。通过10折交叉验证测试，本方案的心血管死亡风险预测模型对应的tp rate和tn rate均在70％以上，具体的性能评估如表5所示：
[0143]
表5
[0144]
性能指标全部属性属性选择tp rate0.5830.750tn rate0.6860.622auc(受试者工作特征曲线下的面积)0.7200.722
[0145]
以此，如表5所示，输入模型的属性从表2中的81个减少到表3中19个，但是模型性
能有较大提升。
[0146]
此外，在一个具体的实施例中，
[0147]
所述步骤s4，包括：基于所述属性数据确定输入值，并基于所述属性数据关联的病例样本的数据确定输出值；基于所述输入值与输出值确定所述深度神经网络中每个连接的权值以及每个神经元的偏差值；基于所述权值与所述偏差值建立针对糖尿病合并冠心病人的心血管死亡风险预测模型。
[0148]
本方案使用深度神经网络建立糖尿病合并冠心病人群得心血管死亡风险预测模型。本发明使用的深度神经网络可以由4层组成，输入层19个节点，两个隐藏层分别为100个节点和10个节点，输出层为2个节点。输入端的网络称为输入层，输出端的网络称为输出层，输入层和输出层中间的所有网络都称为隐藏层。隐藏层数目比较多的神经网络，常常被称为多层感知机；为了应对隐藏层数较多带来的梯度消失问题引入了新的激活函数(如relu、maxout等函数)，而神经网络建模，可以通过已知模型的输出y和输入x，求解每个连接的权值w和每个神经元上的偏差值b来进行。具体的，可以通过反向传播和梯度下降相结合来求解w和b的值。
[0149]
此外，在数据样本量有限的条件下，层数较多的深度神经网络容易发生过拟合现象，为此，可以采用如下方案解决深度神经网络的过拟合问题：(1)及早停止训练。一旦发现在验证集合上性能(性能评价指标可以使用auc值，auc也即受试者工作特征曲线下的面积)下降，立即停止训练；(2)在损失函数(作用是计算预测值与真实值得差异，并利用该差异修正网络参数)上添加正则项。所谓正则项就是添加对模型复杂度的惩罚项，连接权值越大，惩罚越大。这样模型的优化目标就不仅仅是要预测偏差尽量小，同时还要预测所使用的模型尽量简单；(3)使用dropout。就是每次随机选择一些神经元不参与训练，只有在预测的时候这些神经元才生效。这个神经元的输出结果要乘以一个概率值，降低贡献；(4)神经元权重加正则项。限制每个神经元的联结的权重的范数在一个阈值内。
[0150]
通过本方案建立了针对糖尿病合并冠心病人的心血管死亡风险预测模型，且该心血管死亡风险预测模型实现了对危险因素综合评估，以及对相关危险因素贡献度大小进行评估；因此本方案中的心血管死亡风险预测模型有助于建立个体化干预方案，降低冠冠心病患者的死亡风险，提高严重冠心病患者生存率。
[0151]
实施例2
[0152]
为了对本发明进行进一步的说明，本发明实施例2还公开了一种风险预测模型的构建设备，如图2所示，包括
[0153]
获取模块201，用于获取糖尿病合并冠心病患者的多个属性数据；
[0154]
确定模块202，用于基于信息增益评价函数对各所述属性数据进行评估，得到各所述属性数据的信息增益，并基于多个所述属性数据组成数据集；
[0155]
初始化模块203，用于初始化迭代次数n，n＝1；
[0156]
建立模块204，用于第n次基于所述数据集在深度神经网络中进行训练生成模型；
[0157]
评估模块205，用于基于损失函数对所述模型进行性能评估，得到性能值；
[0158]
迭代模块206，用于若所述性能值与初始模型的性能值的差值不大于性能阈值，从所述数据集剔除所述信息增益最小的一个或多个所述属性数据，并将所述迭代次数加1后作为n，启动所述建立模块；所述初始模型为n＝1时基于所有所述属性数据形成的所述数据
集在深度神经网络中进行训练生成的；
[0159]
处理模块207，用于若所述性能值与所述初始模型的性能值的差值大于性能阈值，将第(n-1)次生成的模型作为用于预测糖尿病合并冠心病人的心血管死亡风险的最终模型。
[0160]
在一个具体的实施例中，所述属性数据关联有病例样本的数据；所述确定模块202，用于：
[0161]
基于信息增益评价函数与所述病例样本的数据确定第一信息熵，且基于信息增益评价函数、所述属性数据及所述病例样本的数据确定第二信息熵；所述第一信息熵为所述病例样本中心血管死亡事件的信息熵；所述第二信息熵为不同的所述属性数据取值下的心血管死亡事件的信息熵；
[0162]
基于所述第一信息熵与所述第二信息熵的差值确定各所述属性数据的信息增益。
[0163]
在一个具体的实施例中，所述信息增益评价函数为：
[0164]
g(d,a)＝h(d)-h(d|a)；
[0165]
其中，g(d,a)为所述信息增益；
[0166]
h(d)为心血管死亡事件的信息熵；p
i
包括p0、p1；；n0和n1分别为所述属性数据对应的糖尿病合并冠心病患者中没有发生心血管死亡事件和发生心血管死亡事件的病例数；
[0167]
a为所述属性数据，a为所述属性数据的值；p(a＝a)为所述属性数据的值为a时的概率；h(d|a＝a)为当所述属性数据的值是a时的第二信息熵。
[0168]
在一个具体的实施例中，
[0169]
所述迭代模块206，用于：
[0170]
若所述性能值与初始模型的性能值的差值不大于性能阈值，基于预设的剔除阈值从所述数据集剔除所述信息增益最小的一个或多个所述属性数据；
[0171]
其中，剔除的一个所述属性数据的所述信息增益小于所述剔除阈值，或剔除的多个所述属性数据的所述信息增益的和小于所述剔除阈值。
[0172]
还包括：信息增益阈值模块，用于在将第(n-1)次生成的模型作为所述最终模型之后，将第(n-1)次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第一属性数据；
[0173]
将第(n-2)次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第二属性数据；
[0174]
基于所述第一属性数据与所述第二属性数据两者的信息增益确定针对所述属性数据的信息增益阈值。
[0175]
在一个具体的实施例中，所述获取模块201，用于：
[0176]
获取糖尿病合并冠心病患者的病例样本；其中，所述病例样本包括没有发生心血管死亡的第一病例样本和发生心血管死亡的第二病例样本；所述第一病例样本的数量大于
所述第二病例样本的数量；
[0177]
对所述第一病例样本进行随机下采样，以使所述第一病例样本的数量和所述第二病例样本的数量相同；
[0178]
在所述第一病例样本与第二病例样本中获取属性数据。
[0179]
在一个具体的实施例中，所述建立模块204，用于：
[0180]
基于所述属性数据确定输入值，并基于所述属性数据关联的病例样本的数据确定输出值；
[0181]
基于所述输入值与输出值确定所述深度神经网络中每个连接的权值以及每个神经元的偏差值；
[0182]
基于所述权值与所述偏差值建立针对糖尿病合并冠心病人的心血管死亡风险预测模型。
[0183]
实施例3
[0184]
本发明实施例3还公开了一种终端，包括用于执行实施例1中所述方法的处理器。具体的，本发明实施例3还公开有其他相关特征，具体的相关特征请参见实施例1中的记载。
[0185]
本发明实施例还提供了【实体装置，例如终端设备】，如图3所示，为了便于说明，仅示出了与本发明实施例相关的部分，具体技术细节未揭示的，请参照本发明实施例方法部分。该终端可以为包括、平板电脑、pda(personal digital assistant，个人数字助理)、电脑等任意终端设备，以终端为电脑为例：
[0186]
图3示出的是与本发明实施例提供的终端相关的电脑的部分结构的框图。参考图3，电脑包括：存储器1520、输入单元1530、显示单元1540、传感器1550、音频电路1560、无线保真(wireless fidelity，wifi)模块1570、处理器1580、以及电源1590等部件。本领域技术人员可以理解，图3中示出的结构并不构成对电脑的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
[0187]
下面结合图3对电脑的各个构成部件进行具体的介绍：
[0188]
存储器1520可用于存储软件程序以及模块，处理器1580通过运行存储在存储器1520的软件程序以及模块，从而执行电脑的各种功能应用以及数据处理。存储器1520可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据电脑的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器1520可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
[0189]
输入单元1530可用于接收输入的数字或字符信息，以及产生与电脑的用户设置以及功能控制有关的键信号输入。具体地，输入单元1530可包括触控面板1531以及其他输入设备1532。触控面板1531，也称为触摸屏，可收集用户在其上或附近的触摸操作(比如用户使用手指、触笔等任何适合的物体或附件在触控面板1531上或在触控面板1531附近的操作)，并根据预先设定的程式驱动相应的连接装置。可选的，触控面板1531可包括触摸检测装置和触摸控制器两个部分。其中，触摸检测装置检测用户的触摸方位，并检测触摸操作带来的信号，将信号传送给触摸控制器；触摸控制器从触摸检测装置上接收触摸信息，并将它转换成触点坐标，再送给处理器1580，并能接收处理器1580发来的命令并加以执行。此外，
可以采用电阻式、电容式、红外线以及表面声波等多种类型实现触控面板1531。除了触控面板1531，输入单元1530还可以包括其他输入设备1532。具体地，其他输入设备1532可以包括但不限于物理键盘、功能键(比如音量控制按键、开关按键等)、轨迹球、鼠标、操作杆等中的一种或多种。
[0190]
显示单元1540可用于显示由用户输入的信息或提供给用户的信息以及电脑的各种菜单。显示单元1540可包括显示面板1541，可选的，可以采用液晶显示器(liquid crystal display，lcd)、有机发光二极管(organic light-emitting diode,oled)等形式来配置显示面板1541。进一步的，触控面板1531可覆盖显示面板1541，当触控面板1531检测到在其上或附近的触摸操作后，传送给处理器1580以确定触摸事件的类型，随后处理器1580根据触摸事件的类型在显示面板1541上提供相应的视觉输出。虽然在图3中，触控面板1531与显示面板1541是作为两个独立的部件来实现电脑的输入和输入功能，但是在某些实施例中，可以将触控面板1531与显示面板1541集成而实现电脑的输入和输出功能。
[0191]
音频电路1560、扬声器1561，传声器1562可提供用户与电脑机之间的音频接口。音频电路1560可将接收到的音频数据转换后的电信号，传输到扬声器1561，由扬声器1561转换为声音信号输出；
[0192]
wifi属于短距离无线传输技术，电脑通过wifi模块1570可以帮助用户收发电子邮件、浏览网页和访问流式媒体等，它为用户提供了无线的宽带互联网访问。虽然图3示出了wifi模块1570，但是可以理解的是，其并不属于电脑的必须构成，完全可以根据需要在不改变发明的本质的范围内而省略。
[0193]
处理器1580是电脑的控制中心，利用各种接口和线路连接整个电脑的各个部分，通过运行或执行存储在存储器1520内的软件程序和/或模块，以及调用存储在存储器1520内的数据，执行电脑的各种功能和处理数据，从而对电脑进行整体监控。可选的，处理器1580可包括一个或多个处理单元；优选的，处理器1580可集成应用处理器和调制解调处理器，其中，应用处理器主要处理操作系统、用户界面和应用程序等，调制解调处理器主要处理无线通信。可以理解的是，上述调制解调处理器也可以不集成到处理器1580中。
[0194]
电脑还包括给各个部件供电的电源1590(比如电池)，优选的，电源可以通过电源管理系统与处理器1580逻辑相连，从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。
[0195]
在本发明实施例中，该终端所包括的处理器1580还具有执行以下步骤的功能：
[0196]
步骤s1、获取糖尿病合并冠心病患者的多个属性数据；
[0197]
步骤s2、基于信息增益评价函数对各所述属性数据进行评估，得到各所述属性数据的信息增益，并基于多个所述属性数据组成数据集；
[0198]
步骤s3、初始化迭代次数n，n＝1；
[0199]
步骤s4、第n次基于所述数据集在深度神经网络中进行训练生成模型；
[0200]
步骤s5、基于损失函数对所述模型进行性能评估，得到性能值；
[0201]
步骤s6、若所述性能值与初始模型的性能值的差值不大于性能阈值，从所述数据集剔除所述信息增益最小的一个或多个所述属性数据，并将所述迭代次数加1后作为n，返回步骤s4；所述初始模型为n＝1时基于所有所述属性数据形成的所述数据集在深度神经网络中进行训练生成的；
[0202]
步骤s7、若所述性能值与所述初始模型的性能值的差值大于性能阈值，将第n-1次生成的模型作为用于预测糖尿病合并冠心病人的心血管死亡风险的最终模型。
[0203]
进一步的，所述属性数据关联有病例样本的数据；所述步骤s2，包括：
[0204]
基于信息增益评价函数与所述病例样本的数据确定第一信息熵，且基于信息增益评价函数、所述属性数据及所述病例样本的数据确定第二信息熵；所述第一信息熵为所述病例样本中心血管死亡事件的信息熵；所述第二信息熵为不同的所述属性数据取值下的心血管死亡事件的信息熵；
[0205]
基于所述第一信息熵与所述第二信息熵的差值确定各所述属性数据的信息增益。
[0206]
在一个具体的实施例中，所述信息增益评价函数为：
[0207]
g(d,a)＝h(d)-h(d|a)；
[0208]
其中，g(d,a)为所述信息增益；
[0209]
h(d)为第一信息熵；p
i
包括p0、p1；；n0和n1分别为所述属性数据对应的糖尿病合并冠心病患者中没有发生心血管死亡事件和发生心血管死亡事件的病例数；
[0210]
a为所述属性数据，a为所述属性数据的值；p(a＝a)为所述属性数据的值为a时的概率；h(d|a＝a)为当所述属性数据的值是a时的第二信息熵。
[0211]
所述步骤s6，包括：
[0212]
若所述性能值与初始模型的性能值的差值不大于性能阈值，基于预设的剔除阈值从所述数据集剔除所述信息增益最小的一个或多个所述属性数据；
[0213]
其中，剔除的一个所述属性数据的所述信息增益小于所述剔除阈值，或剔除的多个所述属性数据的所述信息增益的和小于所述剔除阈值。
[0214]
在一个具体的实施例中，在将第n-1次生成的模型作为所述最终模型之后，该方法还包括：
[0215]
将第n-1次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第一属性数据；
[0216]
将第n-2次生成的模型对应的所述数据集中信息增益最小的所述属性数据作为第二属性数据；
[0217]
基于所述第一属性数据与所述第二属性数据两者的信息增益确定针对所述属性数据的信息增益阈值。
[0218]
在一个具体的实施例中，所述步骤s1，包括：
[0219]
获取糖尿病合并冠心病患者的病例样本；其中，所述病例样本包括没有发生心血管死亡的第一病例样本和发生心血管死亡的第二病例样本；所述第一病例样本的数量大于所述第二病例样本的数量；
[0220]
对所述第一病例样本进行随机下采样，以使所述第一病例样本的数量和所述第二病例样本的数量相同；
[0221]
在所述第一病例样本与第二病例样本中获取属性数据。
[0222]
在一个具体的实施例中，所述步骤s4，包括：
[0223]
基于所述属性数据确定输入值，并基于所述属性数据关联的病例样本的数据确定输出值；
[0224]
基于所述输入值与输出值确定所述深度神经网络中每个连接的权值以及每个神经元的偏差值；
[0225]
基于所述权值与所述偏差值建立针对糖尿病合并冠心病人的心血管死亡风险预测模型。
[0226]
实施例4
[0227]
本发明实施例4还公开了一种存储介质，所述存储介质中存储有程序，所述程序运行时执行实施例1-3中任意一项所述方法。具体的，本发明实施例4还公开有其他相关的技术特征，具体其他相关的技术特征请参见实施例1-3中的记载，在此不再进行赘述。
[0228]
在具体的实施过程中，存储介质包括：通用串行总线闪存盘(universal serial bus flash drive，usb)、移动硬盘、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、磁碟或者光盘等各种可以存储程序代码的存储介质。
[0229]
在本发明实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。例如，以上所描述的设备实施例仅仅是示意性的，例如，所述单元或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性或其它的形式。
[0230]
在本发明实施例中的各功能单元可以集成在一个处理单元中，或者各个单元也可以均是独立的物理模块。
[0231]
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备，例如可以是个人计算机，服务器，或者网络设备等，或处理器(processor)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：通用串行总线闪存盘(universal serial bus flash drive)、移动硬盘、rom、ram、磁碟或者光盘等各种可以存储程序代码的介质。
[0232]
本领域技术人员可以理解附图只是一个优选实施场景的示意图，附图中的模块或流程并不一定是实施本发明所必须的。
[0233]
本领域技术人员可以理解实施场景中的装置中的模块可以按照实施场景描述进行分布于实施场景的装置中，也可以进行相应变化位于不同于本实施场景的一个或多个装置中。上述实施场景的模块可以合并为一个模块，也可以进一步拆分成多个子模块。
[0234]
上述本发明序号仅仅为了描述，不代表实施场景的优劣。
[0235]
以上公开的仅为本发明的几个具体实施场景，但是，本发明并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。