CN115394448B - 冠心病运动反应性预测模型的建模方法、模型、设备 - Google Patents
冠心病运动反应性预测模型的建模方法、模型、设备 Download PDFInfo
- Publication number
- CN115394448B CN115394448B CN202211330829.9A CN202211330829A CN115394448B CN 115394448 B CN115394448 B CN 115394448B CN 202211330829 A CN202211330829 A CN 202211330829A CN 115394448 B CN115394448 B CN 115394448B
- Authority
- CN
- China
- Prior art keywords
- heart disease
- coronary heart
- prediction model
- recurrence
- motion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 208000029078 coronary artery disease Diseases 0.000 title claims abstract description 136
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000009257 reactivity Effects 0.000 title claims description 15
- 238000012549 training Methods 0.000 claims abstract description 61
- 230000000694 effects Effects 0.000 claims abstract description 44
- 230000004044 response Effects 0.000 claims abstract description 32
- 238000007781 pre-processing Methods 0.000 claims abstract description 18
- 230000004043 responsiveness Effects 0.000 claims abstract description 18
- 238000012216 screening Methods 0.000 claims abstract description 16
- 230000008859 change Effects 0.000 claims abstract description 14
- 230000006870 function Effects 0.000 claims description 53
- 238000012360 testing method Methods 0.000 claims description 12
- 230000000306 recurrent effect Effects 0.000 claims description 10
- 238000013528 artificial neural network Methods 0.000 claims description 9
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 8
- 238000002790 cross-validation Methods 0.000 claims description 8
- 238000013527 convolutional neural network Methods 0.000 claims description 7
- 238000007473 univariate analysis Methods 0.000 claims description 4
- 238000010200 validation analysis Methods 0.000 claims description 4
- 230000036541 health Effects 0.000 abstract description 2
- 238000013058 risk prediction model Methods 0.000 abstract 1
- 238000011176 pooling Methods 0.000 description 14
- 230000002457 bidirectional effect Effects 0.000 description 4
- 239000003814 drug Substances 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 230000004931 aggregating effect Effects 0.000 description 2
- 238000002554 cardiac rehabilitation Methods 0.000 description 2
- 230000002612 cardiopulmonary effect Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000007477 logistic regression Methods 0.000 description 2
- 230000002503 metabolic effect Effects 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000000284 resting effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 206010007558 Cardiac failure chronic Diseases 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007211 cardiovascular event Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006806 disease prevention Effects 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004220 muscle function Effects 0.000 description 1
- 230000010346 psychosocial stress Effects 0.000 description 1
- 210000002027 skeletal muscle Anatomy 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Biophysics (AREA)
- Databases & Information Systems (AREA)
- Artificial Intelligence (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明涉及运动健康领域,公开一种冠心病运动反应性预测模型的建模方法,包括S101:采集冠心病患者的医疗数据和运动信息数据,筛选出影响冠心病复发的显著特征参数;S102:对显著特征参数进行预处理得到特征数据;S103:建立第一预测模型,对冠心病复发风险预测模型进行训练;S104:基于第一预测模型,结合运动信息数据观察冠心病复发风险变化情况,基于该变化情况评价运动康复方案效果反应;S105:建立第二预测模型,对第二预测模型进行训练。本发明经过两次建模,有针对性地预测患者冠心病复发风险改善效果。
Description
技术领域
本发明涉及生物医学工程技术与运动健康领域,具体涉及一种冠心病运动反应性预测模型及建模方法、电子设备。
背景技术
心脏康复是治疗稳定期心血管疾病以及预防再发心血管事件的重要手段。运动是预防和治疗冠心病最具成本效益的生活方式干预。运动训练和以运动为基础的心脏康复对无症状运动能力、心血管和骨骼肌功能、一般健康生活方式以及减轻抑郁症状和社会心理压力的有益作用现在已得到公认。科学运动能够降低冠心病发病率。然而,哪种运动训练频率、强度、时间、类型的特征是最能降低冠心病复发率的,这在很大程度上仍不清楚。
欧洲心血管疾病预防和康复协会就冠心病患者运动期间的频率、强度、时间、运动类型以及安全性方面制定了建议,建议指出:对于冠状动脉疾病或慢性心脏衰竭患者,每天30-60分钟(每周3-5天)耐力运动训练可能减低复发风险。
而目前医生开具的运动康复方案因遗传及环境因素不同,运动干预存在个体反应差异。根据同一运动处方制定的不同运动方案,有的表现为运动有反应,有的表现为无反应甚至运动抵抗,运动康复方案的效果无法估计。
因此,亟需开发一种冠心病运动反应性预测模型及建模方法、设备。
发明内容
为了使医生对开具的冠心病运动康复方案来预测该运动康复方案对冠心病复发的改善效果,本发明的提供一种冠心病运动反应性预测模型及建模方法、设备。
根据本发明的第一方面,本发明公开一种冠心病运动反应性预测模型的建模方法,包括如下步骤:
S101:采集冠心病患者的医疗数据和运动信息数据,所述医疗数据包括冠心病患者的各项医疗指标和预定时间段内的冠心病的复发情况,通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数,所述冠心病的复发情况为复发或不复发;
S102:对显著特征参数进行预处理得到特征数据;
S103:建立第一预测模型:以特征数据作为第一预测模型的输入,以冠心病的复发情况作为第一预测模型的输出,对第一预测模型进行训练;
S104:基于第一预测模型,在给定运动周期开始到给定运动周期结束期间内,结合运动信息数据观察冠心病复发风险变化情况,基于该变化情况评价运动康复方案效果反应,所述运动康复方案效果反应为有反应或无反应;
S105:建立第二预测模型:以显著特征参数、运动信息数据和第一预测模型的输出作为第二预测模型的输入,以运动康复方案效果反应作为第二预测模型的输出,对第二预测模型进行训练。
作为本发明一示例实施方式,步骤S101中,所述运动信息数据包括运动量、运动强度、运动时间、运动类型、运动类型个数、运动次数和运动时机选择。
作为本发明一示例实施方式,步骤S101中,所述通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数的方法包括:基于单变量分析筛选法分析冠心病患者的各项医疗指标,对于服从正态分布的数据采用独立T检验,对于不服从正态分布的数据采用Wilconxon秩和检验。
作为本发明一示例实施方式,步骤S102中,所述对显著特征参数进行预处理得到特征数据的方法包括:将显著特征参数输入卷积神经网络,得到特征数据。
作为本发明一示例实施方式,步骤S103中,所述对第一预测模型进行训练的方法包括基于递归神经网络算法对第一预测模型进行训练。
作为本发明一示例实施方式,所述递归神经网络采用分类交叉损失函数和带有热重启的Adam优化器进行学习。
作为本发明一示例实施方式,所述分类交叉损失函数采用公式1:
其中,N表示训练样本的数量,M表示分类结果的个数,表示分类结果的真实值,表示预测为该分类结果的概率。Li指的分类信息熵,反映的是分类结果和预测为该分类结果概率的乘积之和,i为自然数,表示是当前是第几个样本,i的取值范围从1~N。
作为本发明一示例实施方式,所述带有热重启的Adam优化器进行学习的方法包括:在每个训练周期后使用验证集验证模型,直到连续训练周期的分类交叉损失函数的数值不再减少。
作为本发明一示例实施方式,步骤S104中,所述基于该变化情况评价运动方案效果的方法包括:
给定运动周期结束时的复发风险低于给定运动周期开始时的复发风险,则认为运动康复方案有反应,否则认为无反应。
作为本发明一示例实施方式,步骤S105中,所述对第二预测模型进行训练的方法包括:
采用XGBoost算法生成一个集合,集合内包括多棵回归树,每棵回归树利用十折交叉验证法进行训练,将对应每棵回归树的叶子节点得分相加,对第i个样本,计算出预测估计值,采用如下公式2:
其中,fk表示每一棵回归树的预测结果,K为回归树的数量,xi表示第i个样本的特征,yi表示第i个样本的第二预测估计值,公式2表示给定一个输入值xi,输出值为K棵回归树的第二预测估计值yi;k为1到K的自然数,i为1到m的自然数,m表示训练用的冠心病患者的样本数量;φ表示求和符号,是对每棵树生成的预测结果的求和,在本实施例中,树的个数为K个。
建立第二迭代函数,第t次第二迭代函数采用如下公式3:
其中,l是一个可以微分的损失函数,yi,yi ~(t-1)表示第t-1次迭代中的第i个样本的实际值yi ~(t-1)与第二预测估计值yi之间的差异,ft(xi)是经过第t次迭代后第i个样本的预测结果,Ω(ft)是复杂度的惩罚函数,t为大于1的自然数。m表示训练样本的个数。m为大于或等于1000的自然数。
作为本发明一示例实施方式,所述建模方法还包括S106:设置第二预测模型的输出,将第二预测模型输出经四舍五入后输出0或1,0代表运动康复方案无反应,1代表运动康复方案有反应。
根据本发明的第二个方面,提供一种冠心病运动反应性预测模型,采用所述的冠心病运动反应性预测模型的建模方法进行构建。
根据本发明一示例实施方式,所述反应性预测模型包括:医疗数据获取模块、显著特征参数获取模块、预处理模块、第一预测模型、运动信息数据获取模块、第二预测模型和反应效果获取模块;
所述医疗数据获取模块用于获取冠心病患者的医疗数据,所述医疗数据包括冠心病患者的各项医疗指标和预定时间段内的冠心病的复发情况,所述冠心病的复发情况为复发或不复发;
所述显著特征参数获取模块用于通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数;
所述预处理模块用于对显著特征参数进行预处理得到特征数据;
所述第一预测模型用于以特征数据作为输入,输出冠心病的复发情况;
所述运动信息数据获取模块用于获取冠心病患者的运动信息数据;
所述反应效果获取模块用于通过第一预测模型和运动信息数据获取模块的数据评价康复方案反应效果。
所述第二预测模型用于以显著特征参数、运动信息数据和冠心病的复发情况作为输入,输出运动康复方案效果反应,所述运动康复方案效果反应为有反应或无反应。
根据本发明的第三个方面,提供一种电子设备,包括:处理器、存储器和总线,
所述处理器和所述存储器通过所述总线完成相互间的通信,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令,能够执行所述的冠心病运动反应性预测模型的建模方法。
本发明提供了一种冠心病运动反应性预测模型及建模方法,第一预测模型采用神经网络算法构建,可以获得冠心病复发风险的信息,第二模型采用XGBoost算法构建,可以获得运动方案反应的信息。通过本方案的模型及建模方法,能够预测患者经过一个给定运动周期的干预后,细化的运动方案带来的效果是否能够使患者冠心病复发风险降低,从而辅助医生形成更有效的运动方案,提高运动干预效果。
附图说明
通过参照附图详细描述其示例实施例,本发明的上述和其它目标、特征及优点将变得更加显而易见。
图1是冠心病运动反应性预测模型的建模方法流程图。
图2是冠心病运动反应性预测模型的结构图。
图3是电子设备的结构图。
图4是卷积神经网络和递归神经网络的连接关系图。
具体实施方式
现在将参考附图更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本发明的描述将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。附图仅为本发明的示意性图解,并非一定是按比例绘制。
此外,所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中,提供许多具体细节从而给出对本发明的实施方式的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而省略所述特定细节中的一个或更多,或者可以采用其它的方法、组元、步骤等。在其它情况下,不详细示出或描述公知结构、方法、实现或者操作以避免喧宾夺主而使得本发明的各方面变得模糊。
本发明的第一种实施方式公开了一种冠心病运动反应性预测模型的建模方法,如图1所示,包括以下步骤:
步骤101:采集冠心病患者的医疗数据和运动信息数据,医疗数据包括冠心病患者的医疗指标和预定时间段内的冠心病的复发情况,通过分析冠心病患者的各项医疗指标和预定时间段内的复发情况筛选出影响冠心病复发的显著特征参数。
因为是在建模阶段,采集的数据均为训练用的数据。训练用的冠心病患者的样本量为n例,n≥1000。
医疗指标包括:一般指标、综合评价指标、运动风险和代谢指标。一般指标包括:年龄、性别、受教育程度和疾病史,一般指标在给定运动周期开始前采集。综合评估指标包括运动心肺能力、静息心电图指标和超声心动指标,综合评估指标在给定运动周期开始前和给定运动周期结束后采集。运动风险包括低危、中危和高危,运动风险通过问卷在给定运动周期开始前进行评估。采集冠心病患者的指标数据时,对所有数据进行归一化处理,保证所有数值的范围在0-1之间。冠心病的复发情况为复发或不复发,预定时间段为一年。
运动信息数据包括:运动量、运动强度、运动时间、运动类型、运动类型个数、运动次数和运动时机选择。
通过分析冠心病患者的各项医疗指标和预定时间段内的复发情况筛选出影响冠心病复发的显著特征参数,具体方法为:基于单变量分析筛选法分析冠心病患者的各项医疗指标,对于服从正态分布的数据采用独立T检验,对于不服从正态分布的数据采用Wilconxon秩和检验。显著特征参数为冠心病患者的各项医疗指标中的一项或多项。
步骤S102:对显著特征参数进行预处理得到特征数据。
对显著特征参数进行预处理得到特征数据的方法包括:将显著特征参数输入卷积神经网络,得到特征数据。
如图4所示,卷积神经网络包括四个一维卷积层、一个最大池化层和一个全局平均池化层。四个一维卷积层用于提取局部特征集,多层的卷积层提取到的特征集就越全局化,其中第四层一维卷积层提取的特征为全局化特征。第一层的一维卷积核大小设置为21,步幅设置为5;第二层的一维卷积核大小设置为21,步幅设置为1;其余两个一维卷积核大小设置为5,步幅设置为1。第一层和第二层的卷积滤波器数为64个,第三层和第四层的卷积滤波器数为128个。池化层用于对局部特征集进行抽样或聚合,可以降低数据的计算复杂度。最大池化层用于对局部特征中的最大值作为最终特征值。最大池化层的参数设置为池大小等于2,步幅等于2。全局平均池化层用于将局部特征中的平均值作为最终特征值。全局平均池化层的参数设置为池大小等于1,步幅等于2。
步骤S103:建立第一预测模型:以特征数据作为第一预测模型的输入,以冠心病的复发情况作为第一预测模型的输出,对第一模型进行训练。
对第一预测模型进行训练的方法包括基于递归神经网络算法采用十折交叉验证法对第一预测模型进行训练。如图4所示,递归神经网络包括一个dropout层、一个双向门控循环单元(GRU)层和一个SoftMax层。Dropout层通过高斯函数,即以概率为0.3随机去掉部分神经元,防止过拟合。
Dropout层采用下述公式4:
其中,其中c为标准方差,反映高斯函数的宽度,b为高斯函数尖峰中心的坐标,0.3表示高斯函数的最高值。c和d都为常数。
该公式4为高斯函数。丢失率为0.3。
双向门控循环单元(GRU)层是为了解决梯度衰竭的问题。双向门控循环单元层的前向步骤中的丢失率为0.3,循环步骤中的丢失率为0.5。
SoftMax层用于形成最终的冠心病复发概率,其采用分类交叉损失函数和带有热重启的Adam优化器进行学习。学习率的初始范围为0.001-0.0001,通过使用优化器估计学习率的最佳范围。
分类交叉损失函数采用公式1:
其中,N表示训练样本的数量,M表示分类结果的个数,表示分类结果的真实值,表示预测为该分类结果的概率。Li指的分类信息熵,反映的是分类结果和预测为该分类结果概率的乘积之和,i为自然数,表示是当前是第几个样本,i的取值范围从1~N。
带有热重启的Adam优化器进行学习的方法包括:在每个训练周期后使用验证集验证模型,直到连续训练周期的分类交叉损失函数的数值不再减少。
每计算一次分类交叉损失函数即为一个训练周期。
可以设置分类交叉损失函数的计算频率,如每次在纳入10个或50个样本后再计算分类交叉损失函数。
通过对第一预测模型进行训练,可以知晓显著特征参数与冠心病复发风险的关系,经过训练后的第一预测模型,当输入显著特征参数时,输出的数值为0-1的数值,该数值反应冠心病复发风险,即为复发风险概率。
步骤S104:基于第一预测模型,在给定运动周期开始到给定运动周期结束期间内,结合运动信息数据观察冠心病复发风险变化情况,基于该变化情况评价运动康复方案效果反应,运动康复方案效果反应为有反应或无反应。
所述观察冠心病复发风险变化情况的方法包括:
分别在给定运动周期开始时和给定运动周期结束时用第一预测模型评价冠心病的复发风险,第一预测模型的输出结果即为复发风险的概率,如40%、50%。
基于该变化情况评价运动方案效果的方法包括:
给定运动周期结束时的复发风险低于给定运动周期开始时的复发风险,则认为运动康复方案有反应,否则认为无反应。
如果给定运动周期结束时的复发风险低于给定运动周期开始时的复发风险,表示在这个运动干预周期内导致冠心病在预定时间段内复发风险变低,即为运动有反应。
给定运动周期指的是在预定时间段内用户执行运动的干预周期,例如一个月、两个月。类似吃药,吃一个疗程的药可以管一年。
S105:建立第二预测模型:以显著特征参数、运动信息数据和第一预测模型的输出作为第二预测模型的输入,以运动康复方案效果反应作为第二预测模型的输出,对第二预测模型进行训练。
第二预测模型的输入为m例训练用的与冠心病复发风险有关的显著特征参数、运动信息数据和第一预测模型的输出,样本量为m,m为大于或等于1000的自然数。采集上述显著特征参数时,对所有数据进行归一化处理,保证所有数值的范围在0-1之间。为了提高计算效率,m例数据不需要重新采集,可以在第一次采集n例数据中获取。
基于Xgboost算法采用十折交叉验证法对第二预测模型进行训练,具体为:
采用Xgboost算法生成一个集合,集合内包括多棵回归树。
每棵回归树利用十折交叉验证法进行训练,将对应每棵回归树的叶子节点得分相加,对第i个样本,计算出第二预测估计值,采用如下公式2:
其中,fk表示每一棵回归树的预测结果,K为回归树的数量,xi表示第i个样本的特征,yi表示第i个样本的第二预测估计值,公式2表示给定一个输入值xi,输出值为K棵回归树的第二预测估计值yi;k为1到K的自然数;i为1到m的自然数,m表示训练用的冠心病患者的样本数量。φ表示求和符号,是对每棵树生成的预测结果的求和,在本实施例中,树的个数为K个。
建立第二迭代函数,第t次第二迭代函数采用如下公式3:
其中,l是一个可以微分的损失函数,yi,yi ~(t-1)表示第t-1次迭代中的第i个样本的实际值yi ~(t-1)与第二预测估计值yi之间的差异,ft(xi)是经过第t次迭代后第i个样本的预测结果,i为1到m的自然数,Ω(ft)是复杂度的惩罚函数,t为大于1的自然数。m表示训练样本的个数。m为大于或等于1000的自然数。
以上公式2和公式3采用的训练参数为:将回归树的最大深度的缺省值设置为6,取值范围为1到正无穷大。回归树的深度即为回归树的层数,通常回归树的层数越深,模型越复杂。将学习率设置为0.3,学习率表示每次调整参数需要走的距离,即为模型参数的最优解,适当的学习率可以准确找到损失函数的最小值。将学习目标进行定义,学习目标即为模型的输出,运动康复方案是否有反应为二分类问题,有反应为1,无反应为0。将迭代次数默认设置为50棵,迭代次数表示模型会用到多少棵树。
将训练参数(最大深度、学习率、学习目标)输入到函数2中进行设置:
其中,max_depth表示树的深度,eta表示学习率,slient表示静默模式,objective表示学习目标参数,binary:logistic表示模型中输出二分类的逻辑回归的概率。
接着调用Xgboost中的train函数进行第二次迭代训练,输入以下代码:
其中,bst表示最后train训练后的结果,50表示迭代次数。
步骤S106:设置第二预测模型的输出,输出的数据是预测运动康复方案反应效果的百分数,如果输出的数值≥0.5则认为运动康复方案有效,如果输出的数值<0.5则认为运动康复方案无效。
本运动康复方案反应性预测模型将输出的数值进行四舍五入,转换成0或1,使得医疗工作人员更明显地看出运动康复方案的效果。
本发明的第二种实施方式公开了一种冠心病运动反应性预测模型,采用第一种实施方式的冠心病运动反应性预测模型的建模方法进行构建。
如图2所示,该反应性预测模型包括:医疗数据获取模块1、显著特征参数获取模块2、预处理模块3、第一预测模型4、运动信息数据获取模块5、第二预测模型6和反应效果获取模块7。
医疗数据获取模块1用于获取冠心病患者的医疗数据,医疗数据包括冠心病患者的各项医疗指标和预定时间段内的冠心病的复发情况,冠心病的复发情况为复发或不复发。医疗指标包括:一般指标、综合评价指标、运动风险和代谢指标。一般指标包括:年龄、性别、受教育程度和疾病史,一般指标在给定运动周期开始前采集。综合评估指标包括运动心肺能力、静息心电图指标和超声心动指标,综合评估指标在给定运动周期开始前和给定运动周期结束后采集。运动风险包括低危、中危和高危,运动风险通过问卷在给定运动周期开始前进行评估。采集冠心病患者的指标数据时,对所有数据进行归一化处理,保证所有数值的范围在0-1之间。冠心病的复发情况为复发或不复发,预定时间段为一年。
显著特征参数获取模块2用于通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数。所述通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数包括:基于单变量分析筛选法分析冠心病患者的各项医疗指标,对于服从正态分布的数据采用独立T检验,对于不服从正态分布的数据采用Wilconxon秩和检验。
预处理模块3用于对显著特征参数进行预处理得到特征数据。对显著特征参数进行预处理得到特征数据包括:将显著特征参数输入卷积神经网络,得到特征数据。如图4所示,卷积神经网络包括四个一维卷积层、一个最大池化层和一个全局平均池化层。四个一维卷积层用于提取局部特征集,多层的卷积层提取到的特征集就越全局化,其中第四层一维卷积层提取的特征为全局化特征。第一层的一维卷积核大小设置为21,步幅设置为5;第二层的一维卷积核大小设置为21,步幅设置为1;其余两个一维卷积核大小设置为5,步幅设置为1。第一层和第二层的卷积滤波器数为64个,第三层和第四层的卷积滤波器数为128个。池化层用于对局部特征集进行抽样或聚合,可以降低数据的计算复杂度。最大池化层用于对局部特征中的最大值作为最终特征值。最大池化层的参数设置为池大小等于2,步幅等于2。全局平均池化层用于将局部特征中的平均值作为最终特征值。全局平均池化层的参数设置为池大小等于1,步幅等于2。
第一预测模型4用于建立特征数据与冠心病复发风险的关系,以特征数据作为输入,输出冠心病的复发情况。建立第一预测模型4:以特征数据作为第一预测模型的输入,以冠心病的复发情况作为第一预测模型4的输出,对第一模型4进行训练。对第一预测模型4进行训练包括基于递归神经网络算法采用十折交叉验证法对第一预测模型4进行训练。如图4所示,递归神经网络包括一个dropout层、一个双向门控循环单元(GRU)层和一个SoftMax层。Dropout层通过高斯函数,即以概率为0.3随机去掉部分神经元,防止过拟合。
Dropout层采用下述公式4:
其中,其中c为标准方差,反映高斯函数的宽度,b为高斯函数尖峰中心的坐标,0.3表示高斯函数的最高值。c和b都为常数。
该公式4为高斯函数。丢失率为0.3。
双向门控循环单元(GRU)层是为了解决梯度衰竭的问题。双向门控循环单元层的前向步骤中的丢失率为0.3,循环步骤中的丢失率为0.5。
SoftMax层用于形成最终的冠心病复发概率,其采用分类交叉损失函数和带有热重启的Adam优化器进行学习。学习率的初始范围为0.001-0.0001,通过使用优化器估计学习率的最佳范围。
分类交叉损失函数采用公式1:
其中,N表示训练样本的数量,M表示分类结果的个数,表示分类结果的真实值,表示预测为该分类结果的概率。Li指的分类信息熵,反映的是分类结果和预测为该分类结果概率的乘积之和,i为自然数,表示是当前是第几个样本,i的取值范围从1~N。
带有热重启的Adam优化器进行学习的方法包括:在每个训练周期后使用验证集验证模型,直到连续训练周期的分类交叉损失函数的数值不再减少。
每计算一次分类交叉损失函数即为一个训练周期。
可以设置分类交叉损失函数的计算频率,如每次在纳入10个或50个样本后再计算分类交叉损失函数。
通过对第一预测模型进行训练,可以知晓显著特征参数与冠心病复发风险的关系,经过训练后的第一预测模型,当输入显著特征参数时,输出的数值为0-1的数值,该数值反应冠心病复发风险,即为复发风险概率。
运动信息数据5获取模块用于获取冠心病患者的运动信息数据。所述运动信息数据包括运动量、运动强度、运动时间、运动类型、运动类型个数、运动次数和运动时机选择。
反应效果获取模块7用于通过第一预测模型4和运动信息数据获取模块5的数据评价康复方案反应效果。基于第一预测模型4,在给定运动周期开始到给定运动周期结束期间内,结合运动信息数据观察冠心病复发风险变化情况,基于该变化情况评价运动康复方案效果反应,运动康复方案效果反应为有反应或无反应。
所述观察冠心病复发风险变化情况包括:
分别在给定运动周期开始时和给定运动周期结束时用第一预测模型评价冠心病的复发风险,第一预测模型的输出结果即为复发风险的概率,如40%、50%。
基于该变化情况评价运动方案效果包括:
给定运动周期结束时的复发风险低于给定运动周期开始时的复发风险,则认为运动康复方案有反应,否则认为无反应。
如果给定运动周期结束时的复发风险低于给定运动周期开始时的复发风险,表示在这个运动干预周期内导致冠心病在预定时间段内复发风险变低,即为运动有反应。
给定运动周期指的是在预定时间段内用户执行运动的干预周期,例如一个月、两个月。类似吃药,吃一个疗程的药可以管一年。
第二预测模型6用于建立显著特征参数、运动信息数据、冠心病复发风险与运动康复方案的关系,以显著特征参数、运动信息数据和冠心病的复发情况作为输入,输出运动康复方案效果反应,运动康复方案效果反应为有反应或无反应。
建立第二预测模型6:以显著特征参数、运动信息数据和第一预测模型的输出作为第二预测模型6的输入,以运动康复方案效果反应作为第二预测模型的输出,对第二预测模型6进行训练。
第二预测模型6的输入为m例训练用的与冠心病复发风险有关的显著特征参数、运动信息数据和第一预测模型4的输出,样本量为m,m为大于或等于1000的自然数。采集上述显著特征参数时,对所有数据进行归一化处理,保证所有数值的范围在0-1之间。为了提高计算效率,m例数据不需要重新采集,可以在第一次采集n例数据中获取。
基于Xgboost算法采用十折交叉验证法对第二预测模型6进行训练,具体为:
采用Xgboost算法生成一个集合,集合内包括多棵回归树。
每棵回归树利用十折交叉验证法进行训练,将对应每棵回归树的叶子节点得分相加,对第i个样本,计算出第二预测估计值,采用如下公式2:
其中,fk表示每一棵回归树的预测结果,K为回归树的数量,xi表示第i个样本的特征,yi表示第i个样本的第二预测估计值,公式1表示给定一个输入值xi,输出值为K棵回归树的第二预测估计值yi;k为1到K的自然数,i为1到m的自然数,m表示训练用的冠心病患者的样本数量。φ表示求和符号,是对每棵树生成的预测结果的求和,在本实施例中,树的个数为K个。
建立第二迭代函数,第t次第二迭代函数采用如下公式3:
其中,l是一个可以微分的损失函数,yi,yi ~(t-1)表示第t-1次迭代中的第i个样本的实际值yi ~(t-1)与第二预测估计值yi之间的差异,ft(xi)是经过第t次迭代后第i个样本的预测结果,i为1到m的自然数,Ω(ft)是复杂度的惩罚函数,t为大于1的自然数。m表示训练样本的个数。m为大于或等于1000的自然数。
以上公式2和公式3采用的训练参数为:将回归树的最大深度的缺省值设置为6,取值范围为1到正无穷大。回归树的深度即为回归树的层数,通常回归树的层数越深,模型越复杂。将学习率设置为0.3,学习率表示每次调整参数需要走的距离,即为模型参数的最优解,适当的学习率可以准确找到损失函数的最小值。将学习目标进行定义,学习目标即为模型的输出,运动康复方案是否有反应为二分类问题,有反应为1,无反应为0。将迭代次数默认设置为50棵,迭代次数表示模型会用到多少棵树。
将训练参数(最大深度、学习率、学习目标)输入到函数2中进行设置:
其中,max_depth表示树的深度,eta表示学习率,slient表示静默模式,objective表示学习目标参数,binary:logistic表示模型中输出二分类的逻辑回归的概率。
接着调用Xgboost中的train函数进行第二次迭代训练,输入以下代码:
其中,bst表示最后train训练后的结果,50表示迭代次数。
设置第二预测模型6的输出,输出的数据是预测运动康复方案反应效果的百分数,如果输出的数值≥0.5则认为运动康复方案有效,如果输出的数值<0.5则认为运动康复方案无效。
本运动康复方案反应性预测模型将输出的数值进行四舍五入,转换成0或1,使得医疗工作人员更明显地看出运动康复方案的效果。
通过医疗数据获取模块1获取医疗指标和冠心病复发情况,通过显著特征参数获取模块2筛选出显著特征参数2,通过预处理模块3得到特征数据,将特征数据和冠心病复发情况训练第一预测模型4,反应效果获取模块7通过第一预测模型4和运动信息数据获取模块5评价康复反应效果,通过显著特征参数、运动信息数据、复发情况和反应效果对第二预测模型6进行训练,最终得到反应性预测模型。
采用该反应性预测模型对冠心病患者的运动康复方案进行预测:
获取测试的冠心病患者的样本,测试样本≥100例。将测试样本输入冠心病运动反应性预测模型中,调用Xgboost的predict函数,输入以下代码:
如果测试结果为0,则表示该运动康复方案无效,如果测试结果为1,则表示该运动康复方案有效。
本发明的第三种实施方式公开了一种电子设备,如图3所示,包括:处理器7、存储器8和总线9;
处理器7和存储器8通过总线9完成相互间的通信,存储器8存储有可被处理器7执行的程序指令,处理器7调用程序指令,能够执行第一种实施方式的冠心病运动反应性预测模型的建模方法。
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节,也不限制该发明仅为所述的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。
Claims (8)
1.一种冠心病运动反应性预测模型的建模方法,其特征在于,包括如下步骤:
S101:采集冠心病患者的医疗数据和运动信息数据,所述医疗数据包括冠心病患者的各项医疗指标和预定时间段内的冠心病的复发情况,通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数,所述冠心病的复发情况为复发或不复发;
S102:对显著特征参数进行预处理得到特征数据;
S103:建立第一预测模型:以特征数据作为第一预测模型的输入,以冠心病的复发情况作为第一预测模型的输出,对第一预测模型进行训练;
S104:基于第一预测模型,在给定运动周期开始到给定运动周期结束期间内,结合运动信息数据观察冠心病复发风险变化情况,基于该变化情况评价运动康复方案效果反应,所述运动康复方案效果反应为有反应或无反应;
S105:建立第二预测模型:以显著特征参数、运动信息数据和第一预测模型的输出作为第二预测模型的输入,以运动康复方案效果反应作为第二预测模型的输出,对第二预测模型进行训练;
步骤S103中,所述对第一预测模型进行训练的方法包括基于递归神经网络算法对第一预测模型进行训练;所述递归神经网络采用分类交叉损失函数和带有热重启的Adam优化器进行学习;所述分类交叉损失函数采用公式1:
所述带有热重启的Adam优化器进行学习的方法包括:在每个训练周期后使用验证集验证模型,直到连续训练周期的分类交叉损失函数的数值不再减少;
步骤S105中,所述对第二预测模型进行训练的方法包括:
采用XGBoost算法生成一个集合,集合内包括多棵回归树,每棵回归树利用十折交叉验证法进行训练,将对应每棵回归树的叶子节点得分相加,对第i个样本,计算出第二预测估计值,采用如下公式2:
其中,fk表示每一棵回归树的预测结果,K为回归树的数量,xi表示第i个样本的特征,yi表示第i个样本的第二预测估计值,公式2表示给定一个输入值xi,输出值为K棵回归树的第二预测估计值yi;k为1到K的自然数,i为1到m的自然数,m表示训练用的冠心病患者的样本数量,φ表示求和符号,是对每棵树生成的预测结果的求和;
建立第二迭代函数,第t次第二迭代函数采用如下公式3:
其中,l是一个可以微分的损失函数,yi,yi ~(t-1)表示第t-1次迭代中的第i个样本的实际值yi ~(t-1)与第二预测估计值yi之间的差异,ft(xi)是经过第t次迭代后第i个样本的预测结果,Ω(ft)是复杂度的惩罚函数,t为大于1的自然数,m表示训练样本的个数。
2.根据权利要求1所述的冠心病运动反应性预测模型的建模方法,其特征在于,步骤S101中,所述运动信息数据包括:运动量、运动强度、运动时间、运动类型、运动类型个数、运动次数和运动时机选择。
3.根据权利要求1所述的冠心病运动反应性预测模型的建模方法,其特征在于,步骤S101中,所述通过分析冠心病患者的各项医疗指标和预定时间段内的复发情况筛选出影响冠心病复发的显著特征参数的方法包括:基于单变量分析筛选法分析冠心病患者的各项医疗指标,对于服从正态分布的数据采用独立T检验,对于不服从正态分布的数据采用Wilconxon秩和检验。
4.根据权利要求1所述的冠心病运动反应性预测模型的建模方法,其特征在于,步骤S102中,所述对显著特征参数进行预处理得到特征数据的方法包括:将显著特征参数输入卷积神经网络,得到特征数据。
5.根据权利要求1所述的冠心病运动反应性预测模型的建模方法,其特征在于,步骤S104中,所述基于该变化情况评价运动康复方案效果的方法包括:
给定运动周期结束时的复发风险低于给定运动周期开始时的复发风险,则认为运动康复方案有反应,否则认为无反应。
6.根据权利要求1所述的冠心病运动反应性预测模型的建模方法,其特征在于,还包括S106:设置第二预测模型的输出,将第二预测模型输出经四舍五入后输出0或1,0代表运动康复方案无反应,1代表运动康复方案有反应。
7.一种冠心病运动反应性预测模型,其特征在于,采用权利要求1-6中任一项所述的冠心病运动反应性预测模型的建模方法进行构建;
所述反应性预测模型包括:医疗数据获取模块、显著特征参数获取模块、预处理模块、第一预测模型、运动信息数据获取模块、第二预测模型和反应效果获取模块;
所述医疗数据获取模块用于获取冠心病患者的医疗数据,所述医疗数据包括冠心病患者的各项医疗指标和预定时间段内的冠心病的复发情况,所述冠心病的复发情况为复发或不复发;
所述显著特征参数获取模块用于通过分析冠心病患者的各项医疗指标和复发情况筛选出影响冠心病复发的显著特征参数;
所述预处理模块用于对显著特征参数进行预处理得到特征数据;
所述第一预测模型用于以特征数据作为输入,输出冠心病的复发情况;
所述运动信息数据获取模块用于获取冠心病患者的运动信息数据;
所述反应效果获取模块用于通过第一预测模型和运动信息数据获取模块的数据评价康复方案反应效果;
所述第二预测模型用于以显著特征参数、运动信息数据和冠心病的复发情况作为输入,输出运动康复方案效果反应,所述运动康复方案效果反应为有反应或无反应。
8.一种电子设备,其特征在于,包括:处理器、存储器和总线,
所述处理器和所述存储器通过所述总线完成相互间的通信,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令,能够执行如权利要求1-6中任一项所述的冠心病运动反应性预测模型的建模方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211330829.9A CN115394448B (zh) | 2022-10-28 | 2022-10-28 | 冠心病运动反应性预测模型的建模方法、模型、设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211330829.9A CN115394448B (zh) | 2022-10-28 | 2022-10-28 | 冠心病运动反应性预测模型的建模方法、模型、设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115394448A CN115394448A (zh) | 2022-11-25 |
CN115394448B true CN115394448B (zh) | 2023-03-28 |
Family
ID=84115035
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211330829.9A Active CN115394448B (zh) | 2022-10-28 | 2022-10-28 | 冠心病运动反应性预测模型的建模方法、模型、设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115394448B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117409930B (zh) * | 2023-12-13 | 2024-02-13 | 江西为易科技有限公司 | 基于ai技术实现的医疗康复数据处理方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4057297A1 (en) * | 2021-03-12 | 2022-09-14 | Dassault Systèmes | Machine-learning with respect to multi-state model of an illness |
CN115240845A (zh) * | 2022-07-19 | 2022-10-25 | 中山大学 | 一种基于体外反搏的冠心病分类方法、系统及装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10580131B2 (en) * | 2017-02-23 | 2020-03-03 | Zebra Medical Vision Ltd. | Convolutional neural network for segmentation of medical anatomical images |
CN110998744B (zh) * | 2017-08-01 | 2024-04-05 | 西门子医疗有限公司 | 针对弥漫性和串联性病变中冠状动脉疾病的非侵入性评估和治疗指导 |
CN109117864B (zh) * | 2018-07-13 | 2020-02-28 | 华南理工大学 | 基于异构特征融合的冠心病风险预测方法、模型及系统 |
CN109754852A (zh) * | 2019-01-08 | 2019-05-14 | 中南大学 | 基于电子病历的心血管疾病风险预测方法 |
CN111292853B (zh) * | 2020-01-15 | 2023-05-26 | 长春理工大学 | 基于多参数的心血管疾病风险预测网络模型及其构建方法 |
CN112270994B (zh) * | 2020-10-14 | 2021-08-17 | 中国医学科学院阜外医院 | 一种风险预测模型的构建方法、设备、终端及存储介质 |
CN113345584A (zh) * | 2021-06-07 | 2021-09-03 | 医渡云(北京)技术有限公司 | 疾病预后风险预测模型训练方法、装置及电子设备 |
-
2022
- 2022-10-28 CN CN202211330829.9A patent/CN115394448B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4057297A1 (en) * | 2021-03-12 | 2022-09-14 | Dassault Systèmes | Machine-learning with respect to multi-state model of an illness |
CN115240845A (zh) * | 2022-07-19 | 2022-10-25 | 中山大学 | 一种基于体外反搏的冠心病分类方法、系统及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN115394448A (zh) | 2022-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106778014B (zh) | 一种基于循环神经网络的患病风险预测建模方法 | |
CN111292853B (zh) | 基于多参数的心血管疾病风险预测网络模型及其构建方法 | |
CN113421652B (zh) | 对医疗数据进行分析的方法、训练模型的方法及分析仪 | |
CN110334843B (zh) | 一种时变注意力改进Bi-LSTM的住院就医行为预测方法及装置 | |
CN111445978B (zh) | 运动方案反应性预测模型及建模方法、电子设备 | |
Li et al. | Identifying informative risk factors and predicting bone disease progression via deep belief networks | |
CN110619322A (zh) | 一种基于多流态卷积循环神经网络的多导联心电异常信号识别方法及系统 | |
CN110111885B (zh) | 属性预测方法、装置、计算机设备及计算机可读存储介质 | |
WO2022166158A1 (zh) | 一种基于卷积生存网络的血透并发症长期风险预测系统 | |
CN111370084A (zh) | 基于BiLSTM的电子健康记录表示学习方法及系统 | |
CN114203295B (zh) | 脑卒中风险预测干预方法及系统 | |
CN115394448B (zh) | 冠心病运动反应性预测模型的建模方法、模型、设备 | |
Kumar et al. | Hybrid Bijective soft set-Neural network for ECG arrhythmia classification | |
CN111297327B (zh) | 一种睡眠分析方法、系统、电子设备及存储介质 | |
CN112768074A (zh) | 一种基于人工智能的重疾风险预测方法及系统 | |
CN109887606B (zh) | 一种基于注意力的双向递归神经网络的诊断预测方法 | |
CN115376638A (zh) | 一种基于多源健康感知数据融合的生理特征数据分析方法 | |
CN112329921B (zh) | 基于深度表征学习和强化学习的利尿剂剂量推理设备 | |
CN111445980B (zh) | 膳食方案反应性预测模型及建模方法、电子设备 | |
CN117010971B (zh) | 一种基于人像识别的智能健康险提供方法及系统 | |
CN111916179A (zh) | 一种基于人工智能自适应个体体征进行“量身定制”饮食滋补模型方法 | |
CN115547502B (zh) | 基于时序数据的血透病人风险预测装置 | |
CN115147768B (zh) | 一种跌倒风险评估方法及系统 | |
An et al. | PARSE: A personalized clinical time-series representation learning framework via abnormal offsets analysis | |
Aswini et al. | Exploratory Data Analysis for Social Big Data Using Regression and Recurrent Neural Networks. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |