CN109243604B - 一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统 - Google Patents
一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统 Download PDFInfo
- Publication number
- CN109243604B CN109243604B CN201811076751.6A CN201811076751A CN109243604B CN 109243604 B CN109243604 B CN 109243604B CN 201811076751 A CN201811076751 A CN 201811076751A CN 109243604 B CN109243604 B CN 109243604B
- Authority
- CN
- China
- Prior art keywords
- model
- kawasaki disease
- data
- neural network
- risk assessment
- 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
- 208000011200 Kawasaki disease Diseases 0.000 title claims abstract description 92
- 208000001725 mucocutaneous lymph node syndrome Diseases 0.000 title claims abstract description 92
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 57
- 238000012502 risk assessment Methods 0.000 title claims abstract description 41
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 33
- 238000010276 construction Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 63
- 238000003745 diagnosis Methods 0.000 claims abstract description 37
- 238000012795 verification Methods 0.000 claims abstract description 35
- 238000012549 training Methods 0.000 claims abstract description 27
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims abstract description 11
- 201000010099 disease Diseases 0.000 claims abstract description 10
- 238000011156 evaluation Methods 0.000 claims abstract description 10
- 238000012216 screening Methods 0.000 claims abstract description 10
- 238000012360 testing method Methods 0.000 claims description 19
- 230000035945 sensitivity Effects 0.000 claims description 13
- 238000013210 evaluation model Methods 0.000 claims description 10
- 206010037660 Pyrexia Diseases 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 102100032752 C-reactive protein Human genes 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 108010088751 Albumins Proteins 0.000 claims description 2
- 102000009027 Albumins Human genes 0.000 claims description 2
- 108010074051 C-Reactive Protein Proteins 0.000 claims description 2
- 108010028780 Complement C3 Proteins 0.000 claims description 2
- 102000016918 Complement C3 Human genes 0.000 claims description 2
- 108010049003 Fibrinogen Proteins 0.000 claims description 2
- 102000008946 Fibrinogen Human genes 0.000 claims description 2
- 102000006395 Globulins Human genes 0.000 claims description 2
- 108010044091 Globulins Proteins 0.000 claims description 2
- 102000007584 Prealbumin Human genes 0.000 claims description 2
- 108010071690 Prealbumin Proteins 0.000 claims description 2
- 239000006185 dispersion Substances 0.000 claims description 2
- 229940012952 fibrinogen Drugs 0.000 claims description 2
- 229940027941 immunoglobulin g Drugs 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims 1
- 238000011282 treatment Methods 0.000 abstract description 12
- 230000002265 prevention Effects 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 9
- 208000024891 symptom Diseases 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 3
- 238000013145 classification model Methods 0.000 description 3
- 210000004351 coronary vessel Anatomy 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 208000019622 heart disease Diseases 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 238000011369 optimal treatment Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 206010015150 Erythema Diseases 0.000 description 2
- 208000032912 Local swelling Diseases 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000031225 myocardial ischemia Diseases 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 102100027211 Albumin Human genes 0.000 description 1
- 208000023275 Autoimmune disease Diseases 0.000 description 1
- 206010011071 Coronary artery aneurysm Diseases 0.000 description 1
- 206010065420 Coronary artery dilatation Diseases 0.000 description 1
- 206010013700 Drug hypersensitivity Diseases 0.000 description 1
- 206010015218 Erythema multiforme Diseases 0.000 description 1
- 208000010201 Exanthema Diseases 0.000 description 1
- 208000015220 Febrile disease Diseases 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 206010020565 Hyperaemia Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 244000132436 Myrica rubra Species 0.000 description 1
- 235000013559 Schnittsellerie Nutrition 0.000 description 1
- 244000169997 Schnittsellerie Species 0.000 description 1
- 206010040047 Sepsis Diseases 0.000 description 1
- 206010040844 Skin exfoliation Diseases 0.000 description 1
- 206010042434 Sudden death Diseases 0.000 description 1
- 208000004732 Systemic Vasculitis Diseases 0.000 description 1
- 206010047115 Vasculitis Diseases 0.000 description 1
- 206010001093 acute tonsillitis Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 210000000436 anus Anatomy 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 229960000190 bacillus calmette–guérin vaccine Drugs 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 210000000795 conjunctiva Anatomy 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 201000005311 drug allergy Diseases 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 201000005884 exanthem Diseases 0.000 description 1
- 210000000416 exudates and transudate Anatomy 0.000 description 1
- 210000004904 fingernail bed Anatomy 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000011081 inoculation Methods 0.000 description 1
- 238000003771 laboratory diagnosis Methods 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 201000003265 lymphadenitis Diseases 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 210000004379 membrane Anatomy 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 210000000214 mouth Anatomy 0.000 description 1
- 210000004400 mucous membrane Anatomy 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 210000003800 pharynx Anatomy 0.000 description 1
- 206010037844 rash Diseases 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000013058 risk prediction model Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000001507 sample dispersion Methods 0.000 description 1
- 208000013223 septicemia Diseases 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 208000011580 syndromic disease 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/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- 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/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Abstract
本发明公开了一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统。所述构建方法包括：从样本数据集中提取可用于建模评估的有效样本；从有效样本的特征集合中筛选出符合现场医疗辅助诊断应用的10项特征；将有效样本的不完整数据集随机分割为训练集和验证集；使用神经网络的方法拟合训练集进行模型构建，采用十折交叉验证法，记录最优模型参数；根据ROC曲线使用验证集计算模型分类阈值t，从而构建得到川崎病风险评估模型。本发明还构建相应的川崎病风险评估系统应用于对待评估数据进行评估，得到KDx评分。本发明有助于降低川崎病的误诊率和漏诊率，使患者在发病早期可以获得有效的预防、干预和治疗。
Description
技术领域
本发明涉及一种模型的构建方法，具体来说涉及一种基于神经网络算法的预测川崎病风险的评估模型的构建方法、构建系统，以及评估系统，属于风险评估模型构建技术领域。
背景技术
川崎病又称皮肤黏膜淋巴结综合征，是一种以全身血管炎为主要病变的自身免疫性疾病，目前已累及全球60多个国家。其中冠状动脉是较易受累部位，是原因不明的发热性出疹性疾病，川崎病主要表现为持续性发热5天以上，还包括：(1)两眼结膜出现充血症状，但未出现渗出物；(2)口唇发红，出现杨梅舌，口腔和咽部黏膜存在弥漫性充血症状；(3)皮肤出现多形性红斑和皮疹；部分患儿可出现卡介苗接种处红肿，是一种特异性表现；(4)四肢末端发生变化；若手足存在硬性肿胀，掌跖和指端充血，则为急性期；若指端甲床皮肤移行部位膜状蜕皮，则为恢复期；肛门周围也多见脱皮症状；(5)急性期表现为非化脓性颈部淋巴结肿大，普遍为单侧，直径在1.5cm以上等临床症状。2017年美国心脏病协会(AHA)制定的川崎病诊断标准：若患者发热≥5天，且以上主要条件中≥4项者确诊为川崎病。含有上述若发热≥5天，主要临床表现不足4项，但在超声心动图或者血管造影发现有冠状动脉病变者，也诊断为川崎病。
川崎病高发人群为5岁以下儿童，最主要且严重的并发症是冠状动脉病变，如果不能进行及时诊断和治疗，会对心血管系统造成严重损伤，冠状动脉扩张和动脉瘤均是该病临床上发生率较高的并发症，还可能直接导致患者发生缺血性心脏病和猝死，目前已成为小儿后天获得性心脏病最常见的病因之一，也是成年后缺血性心脏病发生的危险因素。由此可见，KD的早期诊断具有重要地位。
目前的诊断依据必须发热≥5天，且需要等待临床症状出现，辅以实验室诊断和超心电图检查，容易使患儿错过最佳治疗时间。同时川崎病病因病理目前尚不明确，且发病后可引起多种症状，在一定程度上增加了小儿川崎病的诊断难度。由于患儿本身年龄较小，在没有确诊的情况下，治疗风险较大。小儿川崎病治疗的前提就是要明确诊断，这样才能够及时地对患儿实施治疗。目前尚没有特异性的诊断方法，容易造成患儿临床治疗延误。此外，川崎病的临床症状表现复杂多样，发病初期临床症状不明显，而且与临床上的败血症、淋巴结炎、急性扁挑体炎、药物过敏综合征等疾病症状极为相似，早期误诊率较高。出现误诊的患儿很容易延误病情，进而造成更大的危害。
综上所述，难确诊、易误诊是川崎病患者在诊断过程中的两大难题，是川崎病诊断过程中的临床痛点。因此，研发灵敏度高，特异性强的诊断方式成为川崎病诊疗的中急需满足的需求。
基于医疗数据建模的川崎病患病预测模型可以辅助诊断，有助于降低其漏诊率和误诊率，进一步指导其后续治疗过程。目前存在的基于数据的川崎病分类模型多采用线性方法，典型代表为逻辑回归分析方法。因其敏感性、特异性不足而造成川崎病患者漏诊、误诊情况，从而延误患者治疗。
因此，如何对现有的川崎病患病预测模型进行优化，构建一种具有高敏感性、特异性的风险评估模型，已然成为业界研究人员长期以来一直努力的方向。
发明内容
本发明的主要目的在于提供一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统，以克服现有技术中的不足。
本发明的另一目的还在于提供一种基于神经网络算法的川崎病风险评估系统。
为实现前述发明目的，本发明采用的技术方案包括：
本发明实施例提供了一种基于神经网络算法的川崎病风险评估模型的构建方法，其包括：
从样本数据集中提取可用于建模评估模型的有效样本；
从所述有效样本的特征集中筛选出符合现场医疗辅助诊断应用的10项特征；
将所述有效样本的不完整数据集随机分割为训练集和验证集；
使用神经网络的方法拟合训练集进行模型构建，采用十折交叉验证法，记录最优模型参数；同时，根据ROC曲线使用验证集计算模型分类阈值t，从而构建得到川崎病风险评估模型。
本发明实施例还提供了一种基于神经网络算法的川崎病风险评估模型的构建系统，应用于前述的构建方法，其包括：
数据采集模块，至少用于数据采集，获取样本数据集；
数据处理模块，至少用于从样本数据集中提取可用于构建评估模型的有效样本；
模型构建模块，至少用于将所述有效样本的不完整数据集随机分割为训练集和验证集，并使用神经网络的方法拟合训练集，采用十折交叉验证法，记录最优模型参数；
阈值计算模块，至少用于根据ROC曲线使用验证集计算模型分类阈值。
本发明实施例还提供了由前述方法构建得到的基于神经网络算法的川崎病风险评估模型。
本发明实施例还提供了一种基于神经网络算法的川崎病风险评估系统，其包括：
输入模块，至少用于输入待评估数据；
由前述方法构建得到的基于神经网络算法的川崎病风险评估模型，至少用于对该待评估数据进行评估；
显示模块，至少用于显示评估结果，即KDx评分。
1)与现有技术相比，本发明提供的基于神经网络算法的川崎病风险评估模型构建方法及系统，使用与川崎病相关的医疗数据进行系统的统计分析、建模，并给出模型评价方法，神经网络克服了大多数分类器产生的过拟合问题，是一种表现极好的集成分类器，通过该模型能够基于已有的川崎病医疗数据，对疑似川崎病的患者进行科学有效的辅助评估，有助于降低其误诊率和漏诊率，使患者在发病早期可以获得有效的预防、干预，并科学可靠地指导后续治疗过程，为达到最佳治疗效果提供依据，有效地避免了现有诊断方式中因没有高敏感性和特异性的评估模型而造成川崎病患者漏诊、误诊情况，防止延误患者治疗情况的发生；
2)出于对诊断用时的考虑，本发明所选特征项的检测用时较短，大大缩短医生诊断所用时间。并且，特征项选取较少，降低检测所用成本。
3)本发明数据样本量庞大，优势突出。
附图说明
为了更清楚地说明本发明的技术方案，下面将对实施例或现有技术描述中所需要使用的附图进行简单的介绍，显而易见地，下面描述的附图仅仅作为本文发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。
图1是本发明的一典型实施例中一种基于神经网络算法的川崎病风险评估模型的构建方法的流程示意图。
图2是本发明实施例1中基于神经网络算法的川崎病风险评估模型的ROC曲线图。
具体实施方式
如前所述，鉴于现有技术的不足，本案发明人经长期研究和大量实践，得以提出本发明的技术方案。下面结合附图以及本发明的实施例对一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统等作进一步详细的说明。本发明的保护内容包含但不局限于以下实施案例。在不背离发明构思的精神和范围下，本领域技术人员能够想到的变化和优点都被包括在本发明中，并且以所附的权利要求书为保护范围。
人工神经网络是由大量处理单元互联组成的非线性、自适应信息处理系统。它是在现代神经科学研究成果的基础上提出的，试图通过模拟大脑神经网络处理、记忆信息的方式进行信息处理。它的构筑理念是受到生物(人或其他动物)神经网络功能的运作启发而产生的。人工神经网络通常是通过一个基于数学统计学类型的学习方法(LearningMethod)得以优化，所以人工神经网络也是数学统计学方法的一种实际应用，通过统计学的标准数学方法我们能够得到大量的可以用函数来表达的局部结构空间，另一方面在人工智能学的人工感知领域，我们通过数学统计学的应用可以来做人工感知方面的决定问题(也就是说通过统计学的方法，人工神经网络能够类似人一样具有简单的决定能力和简单的判断能力)，这种方法比起正式的逻辑学推理演算更具有优势。
本发明主要基于电子病例中的医疗数据进行建模，使用数据中蕴含的信息对病人患有川崎病的风险进行评估，并将评估结果进行数字化描述，即得到KDx评分。本发明包括针对医疗数据进行建模的数据处理流程和进行川崎病分类预测、分析、数字化等重要方法和结果。本发明结合了医疗数据和数据挖掘方法，是医疗数据与大数据分析方法结合的一种创新，本发明在一定程度上填补了国内医疗数据研究的空白，在利用医疗数据进行川崎病辅助检测分析方面具有创新性。
本发明实施例的一个方面提供的一种基于神经网络算法的川崎病风险评估模型的构建方法，其包括：
从样本数据集中提取可用于建模评估模型的有效样本；
从所述有效样本的特征集中筛选出符合现场医疗辅助诊断应用的10项特征；
将所述有效样本的不完整数据集随机分割为训练集和验证集；
使用神经网络的方法拟合训练集进行模型构建，采用十折交叉验证法，记录最优模型参数；同时，根据ROC曲线使用验证集计算模型分类阈值t，从而构建得到川崎病风险评估模型。
在一些实施例中，所述构建方法包括：
第一步：数据样本选择；从样本数据集中提取可用于建模及模型评估的有效样本，其具体步骤如下：
1.对重复数据进行删除处理；
2.对数据量不足80％的指标进行删除处理；
3.对残缺、错误数据进行中位数填充；
4.对整体样本离差标准化，将样本数据压缩到[0,1]区间内，并消除量纲：
其中，xi为第i个特征向量，maxi、mini分别为第i个特征向量的最大值、最小值，xi *代表经过变换后的特征向量i；
第二步：特征筛选；从构建样本数据的特征集中筛选出符合现场医疗辅助诊断应用的10项特征；其具体步骤如下：
1.对训练集D，使用神经网络默认参数建立模型；
2.将其中一个特征变量Vi的数值加上或减少10％并进行替换，得到新数据集D1,D2；
3.将D1,D2作为仿真数据集，利用步骤1)中的模型进行仿真预测，记录其正确率的差值MIVi＝D1-D2；
4.重复步骤2)～3)，得到每一个变量的MIV值并进行排序；
5.根据步骤4)所得结果，并结合在实际应用于现场医疗辅助诊断中，取得各项特征值所用时间较短者，进行综合比较得到。
第三步：川崎病患病风险预测模型构建；采用神经网络的方法进行模型构建，其步骤如下：
(1)现有不完整数据集和完整数据集：将不完整数据集随机分割为训练集Xrain、验证集Xderivation，比例为1：1～10：1，并以完整数据集作为测试集Xtest；其具体步骤如下：
1.将训练集数据平均分为十部分；
2.取其中九折数据，使用神经网络的方法进行拟合，得到模型；
3.利用步骤2所得模型，对剩余一折的数据集进行预测，并计算其预测误差；
4.改变模型参数，重复步骤2～3；
5.比较预测误差，记录使得预测误差最小的模型所对应的参数，作为最优模型参数。
(2)使用神经网络的方法拟合Xtrain数据集进行模型构建，采用十折交叉验证法，记录最优模型参数；
(3)根据ROC曲线使用验证集计算模型分类阈值t，阈值t计算具体步骤如下：
1.利用最优参数模型，在训练集上建立最优模型；
2.使用验证集观测值，在模型上进行预测，得到分类得分；
3.在[0,1]范围内，选取不同数值作为分类阀阈值，对步骤2所得分类得分进行划分；
4.计算不同分类阀域下，预测的敏感度、特异度和正确率，并绘制ROC曲线图；
5.根据ROC曲线图选取使得同时满足预测的敏感度、特异度和正确率较优的分类阀域。
在一些实施例中，所述10项特征分别为：
A.性别；
B.年龄；
C.C-反应蛋白浓度(CRP g/L)；
D.纤维蛋白原浓度(FG g/L)；
E.白蛋白浓度(ALB g/L)；
F.球蛋白浓度(GLB g/L)；
G.补体C3浓度(C3 g/L)；
H.免疫球蛋白G浓度(IgG g/L)；
I.前白蛋白PAB浓度(PAB g/L)；
J.白球比例(A/G)。
在一些实施例中，训练集(Xrain)与验证集(Xderivation)的分割比例为1：1～10：1。
在一些实施例中，所述构建方法包括：根据ROC曲线使用验证集计算得到模型分类阈值t，KDx评分高于此分类阈值t预测为川崎病高风险，数值越高，代表川崎病患病概率越大；低于此分类阈值t预测为川崎病低风险，数值越低，代表川崎病患病概率越小。
进一步地，所述构建方法还包括：以完整数据集作为测试集(Xtest)，对构建得到的川崎病风险评估模型进行测试。根据计算所得分类阀域t，进行测试集样本的预测分析。
例如，更具体的，根据训练集构建预测模型并对测试集数据进行预测的步骤包括：
1)使用拟合训练集得到的最优神经网络预测模型，对测试集中每个病人预测其分类得分，即KDx评分。分类得分大于t为川崎病患病高风险患者，分类得分小于t为川崎病患病低风险患者；
2)根据测试集的分类得分计算此模型在辅助川崎病诊断中的敏感性、特异性和准确性。
例如，在一些更为具体的实施方案中，获得可用于构建评估模型的有效样本的过程包括：
(a)根据2017年美国心脏病协会(AHA)制定的川崎病诊断标准将样本数据分为川崎病和普通发热疾病两组，对不能明确诊断结果的样本数据进行删除处理；
(b)对重复数据进行删除处理；
(c)对数据量不足80％的指标进行删除处理；
(d)对残缺、错误数据进行中位数填充，从而获得可用于构建评估模型的有效样本。
本发明使用的医疗数据即样本数据集，来源于医院EDC在线电子病例录入系统，包括医嘱、检验、检查、病程、门诊病历数据、院外随访数据、多中心样本数据、标本分子检测数据等多维数据。
在一些更为具体的实施方案中，参见图1所示，一种基于神经网络算法的川崎病风险评估模型的构建方法，具体步骤如下：
1、样本选择
原始数据集为dataset1，不具有明确诊断结果、重复数据、数据量不足80％的病人被从数据集中移除，此时数据集为dataset2。
2、对整体数据离差标准化，将样本数据压缩到[0,1]区间内，并消除量纲，此时数据集为dataset3。
3、特征筛选
对于dataset2进行特征筛选，通过信息增益计算浏览各特征变量的重要性，删去信息增益接近0的特征变量，同时考虑到特征项数值获得时间长短，取获得时间较短的特征项，此时数据集为dataset4。
4、川崎病分类模型构建
1)现有不完整数据集和完整数据集：将不完整数据集随机分割为训练集Xrain、验证集Xderivation，比例为1：1～10：1，并以完整数据集作为测试集Xtest；
2)使用神经网络的方法拟合Xtrain数据集进行模型构建，采用十折交叉验证法，记录最优模型参数；
3)根据ROC曲线使用验证集计算模型分类阈值t。
本发明实施例的另一个方面还提供了一种基于神经网络算法的川崎病风险评估模型的构建系统，应用于前述的构建方法，其包括：
数据采集模块，至少用于数据采集，获取样本数据集；
数据处理模块，至少用于从样本数据集中提取可用于构建评估模型的有效样本；
模型构建模块，至少用于将所述有效样本的不完整数据集随机分割为训练集和验证集，并使用神经网络的方法拟合训练集，采用十折交叉验证法，记录最优模型参数；
阈值计算模块，至少用于根据ROC曲线使用验证集计算模型分类阈值。
本发明实施例的另一个方面还提供了由前述方法构建得到的基于神经网络算法的川崎病风险评估模型。
相应的，本发明实施例的另一个方面还提供了一种基于神经网络算法的川崎病风险评估系统，其包括：
输入模块，至少用于输入待评估数据；
由前述方法构建得到的基于神经网络算法的川崎病风险评估模型，至少用于对该待评估数据进行评估；
显示模块，至少用于显示评估结果，即KDx评分。
综上所述，本发明的模型构建方法及系统，使用与川崎病相关的医疗数据进行系统的统计分析、建模，并给出模型评价方法，通过该模型能够基于已有的川崎病医疗数据，对疑似川崎病的患者进行科学有效的辅助评估，有助于降低其误诊率和漏诊率，使患者在发病早期可以获得有效的预防、干预，并科学可靠地指导后续治疗过程，为达到最佳治疗效果提供依据，有效地避免了现有诊断方式中因没有高敏感性和特异性的评估模型而造成川崎病患者漏诊、误诊情况，防止延误患者治疗情况的发生。
为使本发明的目的、技术方案和优点更加清楚，下面结合若干优选实施例对本发明的技术方案进行进一步具体描述，但本发明并不仅仅局限于下述实施例，该领域技术人员在本发明核心指导思想下做出的非本质改进和调整，仍然属于本发明的保护范围。
实施例1：
为了验证本发明一种基于神经网络算法的川崎病风险评估模型的构建系统的有效性，本实施例选取时间范围为2008.7-2018.3电子病例中42498个病人数据。本实施例采用神经网络方法。
1、数据处理：
原始数据集经过删除处理之后不完整数据集包括8204个样本，完全数据集包含471个样本。根据本发明采用数据集具有形式为：每行表示为一个病人的信息，每列表示为其一特征信息，如ID，组别，性别，年龄，CRP，FG等，数据集格式如表格1。
通过数据样本选择和特征筛选，最终生成数据集包含的8675行，11列特征，如表1所示。
表1
2、最优模型数据
将不完整数据集随机分为训练集(5742)，验证集(2462)，比例为7:3，完整数据集作为测试集(471)，得到最优模型参数如表2所示：
表2
3、选择分类阀域t
用最优参数模型预测验证集，在[0,1]范围内自动随机生成1259个分类阀域，计算可得对应敏感度、特异度和正确率，并绘制ROC曲线图，如图2所示。
选取靠近曲线左上角并使得敏感度、特异度和正确率较优的分类阀域t＝0.5。
4、对预测结果进行数字化打分
以上模型将作为一种川崎病患病风险评估系统，将测试集中的观测值应用到该系统中进行预测。
测试集结果如表3-1和表3-2所示，本实验中，测试集包括471人。
表3-1
表3-2
附注：关于分类问题一些指标解释，对于二分类问题，定义两个分类分别为正类和负类，正类中的每一个对象成为正实例，负类中的每一个对象成为负实例。通常，在预测川崎病时，川崎病样本为正类，其他发烧患者为负类。使用分类模型对测试样本进行预测，会有四种情况，如果一个实例是正类并被预测为真正类(true positive，TP)，如果实例是负类被预测为正类，称之为假正类(false positive，FP)。相应的，如果实例是负类被预测为负类，称之为真负类(true negative，TN)，正实例被预测为负类则为假负类(falsenegative，FN)。
TP：正实例预测为正类数目；
FN：正实例预测为负类数目；
FP：负实例预测为正类数目；
TN：负实例预测为负类数目；
敏感性(sensitivity)：正类中正确预测为正类的实例比例，即TP/(TP+FN)；
特异性(specificity)：负类中被正确预测为负类的实例比例，即TN/(TN+FP)；
阳性预测值(positive predictive value，PPV)：预测为正类的实例中，正实例占得比例，即TP/(TP+FP)。
正确性：在全部实例中被正确预测的实例比例，即(TP+TN)/(TP+FN+TN+FP)。
实验结果
由测试集数据的真实分类情况可知：278人患有川崎病，193为普通发热。将测试集数据运用到最优神经网络模型中，以其观测值预测其响应值的分类概率KDx(如表3-1所示)，并根据分类阀域t＝0.5对该结果进行划分，得到结果：277人被预测为患有川崎病，194人被预测为普通发热。与测试集中的真实分类比较可得：真正类(TP)为259人，真负类(TN)为166人，假正类(FP)为27人，假负类(FN)为19人(如表3-2所示)。
由测试分类结果可得：敏感度(sensitivity)为93.17％，特异性(specificity)为86.01％，阳性预测值(PPV)为90.56％，正确性为90.23％。
综上所述，由以上数据，本发明一种川崎病患病风险评估系统，通过该模型能够基于已有的川崎病医疗数据，对疑似川崎病的患者进行科学有效的辅助评估，有助于降低其误诊率和漏诊率，使患者在发病早期可以获得有效的预防、干预，并科学可靠地指导后续治疗过程，为达到最佳治疗效果提供依据。出于对诊断用时的考虑，本发明所选特征项的检测用时较短，大大缩短医生诊断所用时间。并且，特征项选取较少，降低检测所用成本。本发明数据样本量庞大，优势突出，原始数据集经过删除处理之后不完整数据集包括8204个样本，完全数据集包含471个样本。
以上所述的实施例对本发明的技术方案进行了详细说明，应理解的是以上所述仅为本发明的具体实施例，并不用于限制本发明，凡在本发明的原则范围内所做的任何修改、补充或类似方式替代等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于包括：
从样本数据集中提取可用于建模评估模型的有效样本；
从所述有效样本的特征集中筛选出符合现场医疗辅助诊断应用的10项特征；
将所述有效样本的不完整数据集随机分割为训练集和验证集；
使用神经网络的方法拟合训练集进行模型构建，采用十折交叉验证法，记录最优模型参数；同时，根据ROC曲线使用验证集计算模型分类阈值t，从而构建得到川崎病风险评估模型；
其中，从所述有效样本的特征集中筛选出符合现场医疗辅助诊断应用的10项特征，具体包括：
步骤11.对训练集D，使用神经网络默认参数建立模型；
步骤12.将其中一个特征变量Vi的数值加上或减少10％并进行替换，得到新数据集D1,D2；
步骤13.将D1,D2作为仿真数据集，利用步骤11中的模型进行仿真预测，记录其正确率的差值MIVi＝D1-D2；
步骤14.重复步骤12～13，得到每一个变量的MIV值并进行排序；
步骤15.根据步骤14所得结果，并结合在实际应用于现场医疗辅助诊断中，取得各项特征值所用时间较短者，进行综合比较得到；
采用十折交叉验证法，记录最优模型参数，具体包括：
步骤21.将训练集数据平均分为十部分；
步骤22.取其中九折数据，使用神经网络的方法进行拟合，得到模型：
步骤23.利用步骤22所得模型，对剩余一折的数据集进行预测，并计算其预测误差；
步骤24.改变参数，重复步骤22～23；
步骤25.比较预测误差，记录使得预测误差最小的模型所对应的参数，作为最优模型参数；
根据ROC曲线使用验证集计算模型分类阈值t，具体包括：
步骤31.利用最优参数模型，在训练集上建立最优模型；
步骤32.使用验证集观测值，在模型上进行预测，得到分类得分；
步骤33.在[0,1]范围内，选取不同数值作为分类阀阈值，对步骤32所得分类得分进行划分；
步骤34.计算不同分类阀域下，预测的敏感度、特异度和正确率，并绘制ROC曲线图；
步骤35.根据ROC曲线图选取使得同时满足预测的敏感度、特异度和正确率较优的分类阀域。
2.根据权利要求1所述的基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于：所述10项特征分别为性别、年龄、C-反应蛋白浓度、纤维蛋白原浓度、白蛋白浓度、球蛋白浓度、补体C3浓度、免疫球蛋白G浓度、前白蛋白浓度和白球比例。
3.根据权利要求1所述的基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于：训练集与验证集的分割比例为1：1～10：1。
4.根据权利要求1所述的基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于包括：根据ROC曲线使用验证集计算模型分类阈值t，KDx评分高于分类阈值t预测为川崎病高风险，低于分类阈值t预测为川崎病低风险。
5.根据权利要求1-4中任一项所述的基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于还包括：以完整的数据集作为测试集，对构建得到的川崎病风险评估模型进行预测。
6.根据权利要求1所述的基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于包括：
根据川崎病诊断标准将样本数据集分为川崎病和普通发热疾病两组，对不能明确诊断结果的样本进行删除处理；
对重复数据进行删除处理；
对数据量不足80％的指标进行删除处理；
对残缺、错误数据进行中位数填充；
对整体数据离差标准化，从而获得可用于构建评估模型的有效样本。
7.根据权利要求6所述的基于神经网络算法的川崎病风险评估模型的构建方法，其特征在于：所述样本数据集来源于医院在线电子病例录入系统，包括医嘱、检验、检查、病程、门诊病历数据、院外随访数据、多中心样本数据和标本分子检测数据。
8.一种基于神经网络算法的川崎病风险评估模型的构建系统，其应用于权利要求1-7中任一项所述的构建方法，其包括：
数据采集模块，至少用于数据采集，获取样本数据集；
数据处理模块，至少用于从样本数据集中提取可用于构建评估模型的有效样本；
模型构建模块，至少用于将所述有效样本的不完整数据集随机分割为训练集和验证集，并使用神经网络的方法拟合训练集，采用十折交叉验证法，记录最优模型参数；
阈值计算模块，至少用于根据ROC曲线使用验证集计算模型分类阈值。
9.一种基于神经网络算法的川崎病风险评估系统，其特征在于包括：
输入模块，至少用于输入待评估数据；
由权利要求1-7中任一项所述方法构建得到的基于神经网络算法的川崎病风险评估模型，至少用于对该待评估数据进行评估；
显示模块，至少用于显示评估结果，即KDx评分。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811076751.6A CN109243604B (zh) | 2018-09-14 | 2018-09-14 | 一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811076751.6A CN109243604B (zh) | 2018-09-14 | 2018-09-14 | 一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109243604A CN109243604A (zh) | 2019-01-18 |
| CN109243604B true CN109243604B (zh) | 2021-11-12 |
Family
ID=65058393
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811076751.6A Active CN109243604B (zh) | 2018-09-14 | 2018-09-14 | 一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109243604B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109949942A (zh) * | 2019-01-30 | 2019-06-28 | 深圳市橙月生物科技有限公司 | 基于铁代谢指标的结核病风险预测模型的构建方法和系统 |
| CN111243736B (zh) * | 2019-10-24 | 2023-09-01 | 中国人民解放军海军军医大学第三附属医院 | 一种生存风险评估方法及系统 |
| CN111462042B (zh) * | 2020-03-03 | 2023-06-13 | 西北工业大学 | 癌症预后分析方法及系统 |
| CN112037919B (zh) * | 2020-09-15 | 2024-02-23 | 南京鼓楼医院 | 一种用于甲状腺结节患者乳头状癌的风险评估模型 |
| CN113936804B (zh) * | 2021-08-23 | 2023-03-28 | 四川大学华西医院 | 一种肺癌切除术后持续漏气风险预测模型构建系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106295229A (zh) * | 2016-08-30 | 2017-01-04 | 青岛大学 | 一种基于医疗数据建模的川崎病分级预测方法 |
| CN106339593A (zh) * | 2016-08-31 | 2017-01-18 | 青岛睿帮信息技术有限公司 | 基于医疗数据建模的川崎病分类预测方法 |
| CN107230108A (zh) * | 2017-06-13 | 2017-10-03 | 北京百分点信息科技有限公司 | 业务数据的处理方法及装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050228692A1 (en) * | 2004-04-08 | 2005-10-13 | Hodgdon Darren W | Incentive based health care insurance program |
| CN107590082A (zh) * | 2017-09-13 | 2018-01-16 | 北京润乾信息系统技术有限公司 | 一种事实数据集在外存中而维数据集在内存中的序号计算多对一join的方法 |
-
2018
- 2018-09-14 CN CN201811076751.6A patent/CN109243604B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106295229A (zh) * | 2016-08-30 | 2017-01-04 | 青岛大学 | 一种基于医疗数据建模的川崎病分级预测方法 |
| CN106339593A (zh) * | 2016-08-31 | 2017-01-18 | 青岛睿帮信息技术有限公司 | 基于医疗数据建模的川崎病分类预测方法 |
| CN107230108A (zh) * | 2017-06-13 | 2017-10-03 | 北京百分点信息科技有限公司 | 业务数据的处理方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 基于数据挖掘技术建立的ＢＰ 神经网络模型;樊楚 等;《中国循症儿科杂志》;20170228;第22-26页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109243604A (zh) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109243604B (zh) | 一种基于神经网络算法的川崎病风险评估模型的构建方法及构建系统 | |
| CN109273093B (zh) | 一种川崎病风险评估模型的构建方法及构建系统 | |
| Gjoreski et al. | Machine learning and end-to-end deep learning for the detection of chronic heart failure from heart sounds | |
| CN109065171B (zh) | 基于集成学习的川崎病风险评估模型的构建方法及系统 | |
| Kam et al. | Learning representations for the early detection of sepsis with deep neural networks | |
| CN109273094B (zh) | 一种基于Boosting算法的川崎病风险评估模型的构建方法及构建系统 | |
| Danso et al. | Developing an explainable machine learning-based personalised dementia risk prediction model: A transfer learning approach with ensemble learning algorithms | |
| Zhang et al. | Lstm for septic shock: Adding unreliable labels to reliable predictions | |
| CN109215781B (zh) | 一种基于logistic算法的川崎病风险评估模型的构建方法及构建系统 | |
| CN106339593A (zh) | 基于医疗数据建模的川崎病分类预测方法 | |
| WO2023087917A1 (zh) | 基于多维分层漂移扩散模型的认知决策评估方法及系统 | |
| Eyvazlou et al. | Prediction of metabolic syndrome based on sleep and work-related risk factors using an artificial neural network | |
| Duan et al. | Early prediction of sepsis using double fusion of deep features and handcrafted features | |
| Kim et al. | Visual field prediction using a deep bidirectional gated recurrent unit network model | |
| Hiroyasu et al. | Gender classification of subjects from cerebral blood flow changes using Deep Learning | |
| CN110473631B (zh) | 基于真实世界研究的智能睡眠监测方法和系统 | |
| Nissimagoudar et al. | AlertNet: Deep convolutional-recurrent neural network model for driving alertness detection | |
| Wang et al. | Method of non-invasive parameters for predicting the probability of early in-hospital death of patients in intensive care unit | |
| CN113128654B (zh) | 一种用于冠心病预诊断中的改进型随机森林模型及其预诊断系统 | |
| Shahul et al. | Machine Learning Based Analysis of Sepsis | |
| CN112329921B (zh) | 基于深度表征学习和强化学习的利尿剂剂量推理设备 | |
| Suneetha et al. | Fine tuning bert based approach for cardiovascular disease diagnosis | |
| Xu et al. | Prediction of Diabetes with its Symptoms Based on Machine Learning | |
| Roversi et al. | A Dynamic Bayesian Network model for simulating the progression to diabetes onset in the ageing population | |
| Chatchumni et al. | Performance of the Simple Clinical Score (SCS) and the Rapid Emergency Medicine Score (REMS) to predict severity level and mortality rate among patients with sepsis in the emergency department |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190118 Assignee: Shanghai Qianbei Medical Technology Co.,Ltd. Assignor: BASEPAIR BIOTECHNOLOGY Co.,Ltd. Contract record no.: X2020980002296 Denomination of invention: Construction method and construction system for Kawasaki disease risk assessment model based on neural network algorithm License type: Common License Record date: 20200518 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210706 Address after: 201600 room 406, no.6, Lane 1015, Longteng Road, Songjiang District, Shanghai Applicant after: Daozhi precision medicine technology (Shanghai) Co.,Ltd. Address before: Unit 426, A2 Floor, 218 Xinghu Street, Suzhou Industrial Park, Jiangsu Province Applicant before: BASEPAIR BIOTECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EC01 | Cancellation of recordation of patent licensing contract |
Assignee: Shanghai Qianbei Medical Technology Co.,Ltd. Assignor: BASEPAIR BIOTECHNOLOGY Co.,Ltd. Contract record no.: X2020980002296 Date of cancellation: 20231218 |
|
| EC01 | Cancellation of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190118 Assignee: Shanghai Haoen Medical Technology Co.,Ltd. Assignor: Daozhi precision medicine technology (Shanghai) Co.,Ltd. Contract record no.: X2024310000027 Denomination of invention: A construction method and system for a risk assessment model of Kawasaki disease based on neural network algorithms Granted publication date: 20211112 License type: Common License Record date: 20240306 |
|
| EE01 | Entry into force of recordation of patent licensing contract |