CN115050451A - 败血症临床用药方案自动生成系统 - Google Patents

败血症临床用药方案自动生成系统 Download PDF

Info

Publication number
CN115050451A
CN115050451A CN202210984717.9A CN202210984717A CN115050451A CN 115050451 A CN115050451 A CN 115050451A CN 202210984717 A CN202210984717 A CN 202210984717A CN 115050451 A CN115050451 A CN 115050451A
Authority
CN
China
Prior art keywords
patient
data
clinical
dense
drug
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.)
Granted
Application number
CN202210984717.9A
Other languages
English (en)
Other versions
CN115050451B (zh
Inventor
傅天童
欧阳波
杨善林
乔灵杰
马志敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202210984717.9A priority Critical patent/CN115050451B/zh
Publication of CN115050451A publication Critical patent/CN115050451A/zh
Application granted granted Critical
Publication of CN115050451B publication Critical patent/CN115050451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Medicinal Chemistry (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Accommodation For Nursing Or Treatment Tables (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Abstract

本发明提供一种败血症临床用药方案自动生成系统,涉及医疗技术领域。本发明中,训练模块根据每个分类簇中的历史临床数据,采用Decision‑Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;验证模块将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision‑Transformer深度学习网络进行迭代。在每个子空间中使用Decision Transformer架构进行训练,基于过去真实的用药数据,选取最好治疗效果的药物集,效果更优。

Description

败血症临床用药方案自动生成系统
技术领域
本发明涉及医疗技术领域,具体涉及一种败血症临床用药方案自动生成系统。
背景技术
随着医疗信息化的快速发展,电子病历和健康档案的实行,产生了大量的文档、表格、图像、语音等多媒体信息。利用人工智能技术辅助开展医疗过程,对数据进行整合分析,为提升医疗卫生服务能力,解决医疗资源紧缺带来了新契机。
目前,部分回顾性研究利用过往数据调查研究患者死亡率、统计患者死亡原因等,为医务工作者提供相对准确的医学经验,但由于每位患者自身状态各不相同,不能以经验主义一概而论,其实际意义相对低微。还有部分预测性研究通过EMR或公开数据集提取相应数据,对患者进行个性化分析和诊断,提高患者的生存率。由于如下原因:(1)患者病情较为复杂,用药数量、剂量不确定;(2)临床医生由于医学责任、既往经验等顾虑而更相信自己的判断。大多数研究止步于病情的预测、预后,而未有根据诊断结果进行辅助用药预测。
此外,采用的机器学习训练方案各不相同,算法设计较为复杂或不合适,例如败血症临床用药方案自动生成方案中,对高维数据采用k-means聚类算法,可能导致聚类效果不理想的问题,进而导致药物使用集合估计结果不准确。
发明内容
(一)解决的技术问题
针对现有技术的不足,本发明提供了一种败血症临床用药方案自动生成系统,解决了药物使用集合估计结果不准确的技术问题。
(二)技术方案
为实现以上目的,本发明通过以下技术方案予以实现:
一种败血症临床用药方案自动生成系统,包括:
预处理模块,用于获取并预处理败血症患者的历史临床数据,并将预处理后的历史临床数据划分为训练集和验证集;
聚类模块,用于根据所述训练集,采用子空间聚类算法进行聚类,获取若干分类簇;
训练模块,用于根据每个分类簇中的历史临床数据,采用Decision-Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;
验证模块,用于将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision-Transformer深度学习网络进行迭代;
生成模块,根据败血症患者的当前临床数据,选定患者的存活率最高的网络模型自动生成临床用药方案。
优选的,所述预处理模块具体用于:
S11、用于在重症监护医疗数据集MIMIC-IV中依据国际Spesis-3标准,提取入院后24h内,被诊断为败血症的成年患者,在医院接受诊断及治疗前72h的历史临床数据;
S12、用于平均划分每位患者前72h的历史临床数据为M段,将其编码为具有∆t=72/M时间步长的多维离散时间序列,由此对每位患者得到M*L维数据表,每一个患者的数据记作:
Figure 270840DEST_PATH_IMAGE002
,其中Pn im表示第n个患者的第m条时间序列数据向量,L表示历史临床数据的总类数。
优选的,所述历史临床数据包括47种:
(1)基础数据:患者id、性别,年龄,身高,首次测量体重,身体质量指数;
(2)生理数据:心率,呼吸频率,平均压,收缩压,舒张压;
(3)血液元素值:白蛋白,过量碱,碳酸氢盐,血尿素氮,钙,钾,氯,钠,二氧化碳,肌酐,游离钙,葡萄糖,血红蛋白,国际标准化比值,乳酸,镁,酸碱度,血小板,凝血酶原时间,凝血激活酶时间,血氧饱和度,总胆红素,白细胞计数;
(4)医学诊断标准:急性生理功能评分系统首次评分值,序贯器官衰竭的检测首次评分值,序贯器官衰竭的检测最大评分值,患者于重症监护病房停留日长,患者于重症监护病房死亡率——布尔值,患者于医院内停留日长,患者于医院内死亡率——布尔值;
(5)用药情况:多巴胺,肾上腺素,去甲肾上腺素,苯肾上腺素,氢吗啡酮,超过4小时的静脉输液摄入。
优选的,所述预处理模块执行步骤S12前,还用于
若某一患者的历史临床数据符合以下情形,剔除该患者的数据;
(1)入院时未成年;
(2)最终死亡或存活未被记录;
(3)未有静脉输液和血管升压药使用记录;
(4)缺失部分较多;
(5)超出临床数据正常范围过大;
(6)由于患者病情严重,采取大剂量血液升压药的急救治疗策略,且于12h内死亡。
优选的,所述聚类模块具体用于:
S21、将训练集中所有患者的历史临床数据定义为L*1*N维,记作去id化、去时间序列的N条特征向量,其中N表示训练集对应的患者总数;
S22、初始化特征空间,把空间划分若不重叠的矩形单元,并计算每个网格的密度,根据给定的阈值,识别密集网格和非密集网格,将这N条特征向量的数据集放入初始化特征空间;
S23、遍历训练集,得到1维密集网格单元之后,进行迭代,由k-1维单元产生候选的k维密集网格单元,之后再遍历一次训练集确定真正的k维密集网格单元,直至无法产生新的候选密集k维密集网格单元,迭代结束;
S24、标记相同的密集网格组成分类簇;
S25、重复执行S23和S24,直至产生t个分类簇,其中同一分类簇对应的各个患者属于同一病理特征。
优选的,所述S24中标记相同的密集网格组成分类簇,具体包括:
当发现第一个密集网格时,以该网格开始扩展,若一个网格与已知密集区域内的网格邻接并且其自身也是密集的,则合并网格到该密集区域,直到不再有这样的网格被发现为止。
优选的,所述聚类模块在执行步骤25后,还具体用于:
S26、统计每t个分类簇中患者的死亡率,若大于预设的死亡率阈值,则剔除该部分分类簇,否则保留用于训练模块学习Decision-Transformer深度学习网络。
优选的,所述训练模块具体用于:
对每个分类簇中抽取序列长度为M的历史临床数据,训练与输入标记St对应的预测头,预测离散动作At的交叉熵损失、以及连续动作的均方误差,并对每个时间步的损失进行平均,获取所述离线轨迹数据集,
其中患者数据轨迹记为τ
Figure 860084DEST_PATH_IMAGE004
其中
Figure 905401DEST_PATH_IMAGE006
、st、at分别为t时刻的Reward,state和action,t=1,2,⋯,T;
Figure 707134DEST_PATH_IMAGE006
设置为 t 时刻之后的所有奖励之和,最终患者存活为奖励,患者死亡为惩罚;
state设置为患者的N*L维数据表;
action设置为药物治疗剂量组合。
优选的,所述药物治疗剂量组合具体设置为静脉输液剂量和和血管升压药输入剂量的若干种剂量组合。
优选的,所述若干种剂量组合包括25种:
接受血管升压药剂量和为0、1~50、50~200、200~500、>500mg,与接受静脉输液剂量和为0、0.001~0.08、0.08~0.22、0.22~0.45、>0.45ml之间的两两组合。
(三)有益效果
本发明提供了一种败血症临床用药方案自动生成系统。与现有技术相比,具备以下有益效果:
本发明中,训练模块根据每个分类簇中的历史临床数据,采用Decision-Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;验证模块将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision-Transformer深度学习网络进行迭代。在每个子空间中使用DecisionTransformer架构进行训练,基于过去真实的用药数据,选取最好治疗效果的药物集,效果更优。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种败血症临床用药方案自动生成系统的结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本申请实施例通过提供一种败血症临床用药方案自动生成系统,解决了药物使用集合估计结果不准确的技术问题。
本申请实施例中的技术方案为解决上述技术问题,总体思路如下:
之前的研究在对数据集选取时未考虑部分生理数据值在发病前并未测量,数据缺失量极大。本实施例将每位患者第一次测量心率数据的时间作为判定该患者接受诊断及治疗的开始时间,并提取此后其他生理检测的全部数据,最大程度的确保数据完整性及准确性。
之前的研究使用基于距离进行聚类的k-means++算法。由于本申请中数据维度较高,数据较低维空间中数据分布要稀疏,在高维空间中基于距离来构建簇误差较大。因此本发明实施例采用子空间聚类算法对患者数据进行聚类,更好的对患者状态作出区分。
之前的研究采用传统强化学习中的TD差分学习,本发明实施例在每个子空间中使用Decision Transformer架构进行训练,基于过去真实的用药数据,选取最好治疗效果的药物集,效果更优。
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
实施例:
如图1所示,本发明实施例提供了一种败血症临床用药方案自动生成系统,包括:
预处理模块,用于获取并预处理败血症患者的历史临床数据,并将预处理后的历史临床数据划分为训练集和验证集;
聚类模块,用于根据所述训练集,采用子空间聚类算法进行聚类,获取若干分类簇;
训练模块,用于根据每个分类簇中的历史临床数据,采用Decision-Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;
验证模块,用于将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision-Transformer深度学习网络进行迭代;
生成模块,根据败血症患者的当前临床数据,选定患者的存活率最高的网络模型自动生成临床用药方案。
本发明实施例在每个子空间中使用Decision Transformer架构进行训练,基于过去真实的用药数据,选取最好治疗效果的药物集,效果更优。
接下来将结合具体内容详细介绍上述技术方案的各个模块:
在预处理模块中,其用于获取并预处理败血症患者的历史临床数据,并将预处理后的历史临床数据划分为训练集和验证集;具体用于:
S11、用于在重症监护医疗数据集MIMIC-IV贝斯以色列迪康医学中心2001年6月——2012年10月重症监护室病人的临床数据记录)中依据国际Spesis-3标准,提取入院后24h内,被诊断为败血症的成年患者,在医院接受诊断及治疗前72h的历史临床数据;
S12、用于平均划分每位患者前72h的历史临床数据为M段,将其编码为具有∆t=72/M时间步长的多维离散时间序列,由此对每位患者得到M*L维数据表,每一个患者的数据记作:
Figure 787086DEST_PATH_IMAGE002
,其中Pn im表示第n个患者的第m条时间序列数据向量,L表示历史临床数据的总类数,在此本发明实施取M=18,则∆t=4h。
所述历史临床数据包括L=47种:
(1)基础数据:患者id、性别,年龄,身高,首次测量体重,身体质量指数;
(2)生理数据:心率,呼吸频率,平均压,收缩压,舒张压;
(3)血液元素值:白蛋白,过量碱,碳酸氢盐,血尿素氮,钙,钾,氯,钠,二氧化碳,肌酐,游离钙,葡萄糖,血红蛋白,国际标准化比值,乳酸,镁,酸碱度,血小板,凝血酶原时间,凝血激活酶时间,血氧饱和度,总胆红素,白细胞计数;
(4)医学诊断标准:急性生理功能评分系统首次评分值,序贯器官衰竭的检测首次评分值,序贯器官衰竭的检测最大评分值,患者于重症监护病房停留日长,患者于重症监护病房死亡率——布尔值,患者于医院内停留日长,患者于医院内死亡率——布尔值;
(5)用药情况:多巴胺,肾上腺素,去甲肾上腺素,苯肾上腺素,氢吗啡酮,超过4小时的静脉输液摄入。
特别的,所述预处理模块执行步骤S12前,还用于
若某一患者的历史临床数据符合以下情形,剔除该患者的数据;
(1)入院时未成年;
(2)最终死亡或存活未被记录;
(3)未有静脉输液和血管升压药使用记录;
(4)缺失部分较多;
(5)超出临床数据正常范围过大;
(6)由于患者病情严重,采取大剂量血液升压药的急救治疗策略,且于12h内死亡。
由预处理模块获取的M*L维数据表,能够较为清楚的展示了患者生理状态随着用药剂量而变化的趋势。选取其中80%作为样本训练,其余20%留作模型验证。
在聚类模块中,其用于根据所述训练集,采用子空间聚类算法进行聚类,获取若干分类簇;
由于高维空间中数据较低维空间中数据分布要稀疏,其中数据间距离几乎相等是普遍现象,而传统聚类方法是基于距离进行聚类的,因此在高维空间中无法基于距离来构建簇。所以本发明实施例采用子空间聚类(Subspace Clustering)算法对患者数据进行聚类,它先将原始特征空间分为若干个网格,再以落到某网格中样本点的概率表示该子空间的密度情况。对于密度超过一定阈值的子空间作为密集单元进行保留,而对非密集子空间进行舍弃。
算法流程如下:
(1)将数据划分成m个子空间,识别每个子空间包含的t个簇;
(2)识别这t个簇;
(3)生成描述这些簇的集群CT
具体用于:
S21、将训练集中所有患者的历史临床数据定义为L*1*N维,记作去id化、去时间序列的N条特征向量,其中N表示训练集对应的患者总数,根据上文的介绍可知L=47。
S22、初始化特征空间,把空间划分若不重叠的矩形单元,并计算每个网格的密度,根据给定的阈值,识别密集网格和非密集网格,将这N条特征向量的数据集放入初始化特征空间。
如果一个向量点集(簇)S是47维空间的一个簇,那么S是该空间的任意k−1维子空间中某个簇的一部分。因此,单次迭代过程包括如下S23和S24:
S23、遍历训练集,得到1维密集网格单元之后,进行迭代,由k-1维单元产生候选的k维密集网格单元,之后再遍历一次训练集确定真正的k维密集网格单元,直至无法产生新的候选密集k维密集网格单元,迭代结束。
S24、标记相同的密集网格组成分类簇;具体包括:
当发现第一个密集网格时,以该网格开始扩展,若一个网格与已知密集区域内的网格邻接并且其自身也是密集的,则合并网格到该密集区域,直到不再有这样的网格被发现为止。
S25、重复执行S23和S24,直至产生t个分类簇,其中同一分类簇对应的各个患者属于同一病理特征。
S26、统计每t个分类簇中患者的死亡率,若大于预设的死亡率阈值,则剔除该部分分类簇,否则保留用于训练模块学习Decision-Transformer深度学习网络。
在训练模块中,其用于根据每个分类簇中的历史临床数据,采用Decision-Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;具体用于:
对每个分类簇中抽取序列长度为M(M=18)的历史临床数据,训练与输入标记St对应的预测头,预测离散动作At的交叉熵损失、以及连续动作的均方误差,并对每个时间步的损失进行平均,获取所述离线轨迹数据集,
其中患者数据轨迹记为τ
Figure 16073DEST_PATH_IMAGE004
其中
Figure 814265DEST_PATH_IMAGE007
、st、at分别为t时刻的Reward,state和action,t=1,2,⋯,T;
Figure 419690DEST_PATH_IMAGE007
设置为 t 时刻之后的所有奖励之和,最终患者存活为奖励+100,患者死亡为惩罚-100;
state设置为患者的N*L维数据表;
action设置为药物治疗剂量组合。
如表1所示,所述药物治疗剂量组合具体设置为静脉输液剂量和和血管升压药输入剂量的若干种剂量组合,即包括以下25种:
接受血管升压药剂量和为0、1~50、50~200、200~500、>500mg,与接受静脉输液剂量和为0、0.001~0.08、0.08~0.22、0.22~0.45、>0.45ml之间的两两组合。
表1
Figure 619727DEST_PATH_IMAGE009
上文所言的Decision Transformer对Transformer架构进行较小修改,它对轨迹进行自回归建模——即对状态、动作和奖励序列的联合分布建模,以替代传统RL算法。它根据过往已经收集到的经验——即患者过往的生理状态数据表及治疗过程,使用序列建模目标来训练Transformer模型。
特别的,患者第一条被完整记录的4h时间步长序列被当作初始状态,作为初始生成的条件信息。将每个P n i m *47数据向量经过一个嵌入,再加上位置编码得到tokens,经过Transformer处理得到输出tokens,再经过多层感知机的decoder 层得到输出的action(用药组合),训练网络通过执行每一个当前的给药状态生成操作后,通过患者最终生存或死亡的结果而设置的奖惩设置,以取得的奖励来减少目标回报,并一直重复到这一患者记录的数据状态终止。
本发明实施例将每位患者第一次测量心率数据的时间作为判定该患者接受诊断及治疗的开始时间,并提取此后其他生理检测的全部数据,最大程度的确保数据完整性及准确性。
在验证模块中,其用于将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision-Transformer深度学习网络进行迭代;
在生成模块中,其根据败血症患者的当前临床数据,选定患者的存活率最高的网络模型自动生成临床用药方案,为临床医生提供辅助治疗建议。
综上所述,与现有技术相比,具备以下有益效果:
本发明实施例中,训练模块根据每个分类簇中的历史临床数据,采用Decision-Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;验证模块将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision-Transformer深度学习网络进行迭代。在每个子空间中使用DecisionTransformer架构进行训练,基于过去真实的用药数据,选取最好治疗效果的药物集,效果更优。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种败血症临床用药方案自动生成系统,其特征在于,包括:
预处理模块,用于获取并预处理败血症患者的历史临床数据,并将预处理后的历史临床数据划分为训练集和验证集;
聚类模块,用于根据所述训练集,采用子空间聚类算法进行聚类,获取若干分类簇;
训练模块,用于根据每个分类簇中的历史临床数据,采用Decision-Transformer深度学习网络获取离线轨迹数据集,所述离线轨迹数据集包括预测得到的患者生理状态、药物治疗剂量组合和奖励大小;
验证模块,用于将预测得到的所述药物治疗剂量组合返回验证集进行验证,比较所述验证集中的实际药物治疗剂量与前述预测得到的药物治疗剂量组合比较类似程度,若患者的存活率状况满足预设条件,则将对应的参数返回Decision-Transformer深度学习网络进行迭代;
生成模块,根据败血症患者的当前临床数据,选定患者的存活率最高的网络模型自动生成临床用药方案。
2.如权利要求1所述的败血症临床用药方案自动生成系统,其特征在于,所述预处理模块具体用于:
S11、用于在重症监护医疗数据集MIMIC-IV中依据国际Spesis-3标准,提取入院后24h内,被诊断为败血症的成年患者,在医院接受诊断及治疗前72h的历史临床数据;
S12、用于平均划分每位患者前72h的历史临床数据为M段,将其编码为具有∆t=72/M时间步长的多维离散时间序列,由此对每位患者得到M*L维数据表,每一个患者的数据记作:
Figure 628658DEST_PATH_IMAGE002
,其中Pnim表示第n个患者的第m条时间序列数据向量,L表示历史临床数据的总类数。
3.如权利要求2所述的败血症临床用药方案自动生成系统,其特征在于,所述历史临床数据包括47种:
(1)基础数据:患者id、性别,年龄,身高,首次测量体重,身体质量指数;
(2)生理数据:心率,呼吸频率,平均压,收缩压,舒张压;
(3)血液元素值:白蛋白,过量碱,碳酸氢盐,血尿素氮,钙,钾,氯,钠,二氧化碳,肌酐,游离钙,葡萄糖,血红蛋白,国际标准化比值,乳酸,镁,酸碱度,血小板,凝血酶原时间,凝血激活酶时间,血氧饱和度,总胆红素,白细胞计数;
(4)医学诊断标准:急性生理功能评分系统首次评分值,序贯器官衰竭的检测首次评分值,序贯器官衰竭的检测最大评分值,患者于重症监护病房停留日长,患者于重症监护病房死亡率——布尔值,患者于医院内停留日长,患者于医院内死亡率——布尔值;
(5)用药情况:多巴胺,肾上腺素,去甲肾上腺素,苯肾上腺素,氢吗啡酮,超过4小时的静脉输液摄入。
4.如权利要求2所述的败血症临床用药方案自动生成系统,其特征在于,所述预处理模块执行步骤S12前,还用于
若某一患者的历史临床数据符合以下情形,剔除该患者的数据;
入院时未成年;
最终死亡或存活未被记录;
未有静脉输液和血管升压药使用记录;
缺失部分较多;
超出临床数据正常范围过大;
由于患者病情严重,采取大剂量血液升压药的急救治疗策略,且于12h内死亡。
5.如权利要求2所述的败血症临床用药方案自动生成系统,其特征在于,所述聚类模块具体用于:
S21、将训练集中所有患者的历史临床数据定义为L*1*N维,记作去id化、去时间序列的N条特征向量,其中N表示训练集对应的患者总数;
S22、初始化特征空间,把空间划分若不重叠的矩形单元,并计算每个网格的密度,根据给定的阈值,识别密集网格和非密集网格,将这N条特征向量的数据集放入初始化特征空间;
S23、遍历训练集,得到1维密集网格单元之后,进行迭代,由k-1维单元产生候选的k维密集网格单元,之后再遍历一次训练集确定真正的k维密集网格单元,直至无法产生新的候选密集k维密集网格单元,迭代结束;
S24、标记相同的密集网格组成分类簇;
S25、重复执行S23和S24,直至产生t个分类簇,其中同一分类簇对应的各个患者属于同一病理特征。
6.如权利要求5所述的败血症临床用药方案自动生成系统,其特征在于,所述S24中标记相同的密集网格组成分类簇,具体包括:
当发现第一个密集网格时,以该网格开始扩展,若一个网格与已知密集区域内的网格邻接并且其自身也是密集的,则合并网格到该密集区域,直到不再有这样的网格被发现为止。
7.如权利要求5所述的败血症临床用药方案自动生成系统,其特征在于,所述聚类模块在执行步骤25后,还具体用于:
S26、统计每t个分类簇中患者的死亡率,若大于预设的死亡率阈值,则剔除该部分分类簇,否则保留用于训练模块学习Decision-Transformer深度学习网络。
8.如权利要求5~7任一项所述的败血症临床用药方案自动生成系统,其特征在于,所述训练模块具体用于:
对每个分类簇中抽取序列长度为M的历史临床数据,训练与输入标记St对应的预测头,预测离散动At作的交叉熵损失、以及连续动作的均方误差,并对每个时间步的损失进行平均,获取所述离线轨迹数据集,
其中患者数据轨迹记为τ
Figure 280220DEST_PATH_IMAGE004
其中
Figure 404165DEST_PATH_IMAGE006
、st、at分别为t时刻的Reward,state和action,t=1,2,⋯,T;
Figure 268215DEST_PATH_IMAGE006
设置为t时刻之后的所有奖励之和,最终患者存活为奖励,患者死亡为惩罚;
state设置为患者的N*L维数据表;
action设置为药物治疗剂量组合。
9.如权利要求8所述的败血症临床用药方案自动生成系统,其特征在于,所述药物治疗剂量组合具体设置为静脉输液剂量和和血管升压药输入剂量的若干种剂量组合。
10.如权利要求9所述的败血症临床用药方案自动生成系统,其特征在于,所述若干种剂量组合包括25种:
接受血管升压药剂量和为0、1~50、50~200、200~500、>500mg,与接受静脉输液剂量和为0、0.001~0.08、0.08~0.22、0.22~0.45、>0.45ml之间的两两组合。
CN202210984717.9A 2022-08-17 2022-08-17 败血症临床用药方案自动生成系统 Active CN115050451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210984717.9A CN115050451B (zh) 2022-08-17 2022-08-17 败血症临床用药方案自动生成系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210984717.9A CN115050451B (zh) 2022-08-17 2022-08-17 败血症临床用药方案自动生成系统

Publications (2)

Publication Number Publication Date
CN115050451A true CN115050451A (zh) 2022-09-13
CN115050451B CN115050451B (zh) 2022-11-04

Family

ID=83168380

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210984717.9A Active CN115050451B (zh) 2022-08-17 2022-08-17 败血症临床用药方案自动生成系统

Country Status (1)

Country Link
CN (1) CN115050451B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116453706A (zh) * 2023-06-14 2023-07-18 之江实验室 一种基于强化学习的血液透析方案制定方法及系统
CN116779186A (zh) * 2023-08-23 2023-09-19 寿光市人民医院 一种用于心律失常患者的环磷腺苷用药数据信息处理系统
CN117275661A (zh) * 2023-11-23 2023-12-22 太原理工大学 一种基于深度强化学习的肺癌患者用药预测方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170193176A1 (en) * 2015-12-30 2017-07-06 Insilico Medicine, Inc. System, Method, and Software for Improved Drug Efficacy and Safety in a Patient
KR101886374B1 (ko) * 2017-08-16 2018-08-07 재단법인 아산사회복지재단 딥러닝 기반의 패혈증 조기 감지방법 및 프로그램
WO2020043922A1 (en) * 2018-08-31 2020-03-05 Novo Nordisk A/S Retrospective horizon based insulin dose prediction
CN111696678A (zh) * 2020-06-15 2020-09-22 中南大学 一种基于深度学习的用药决策方法和系统
WO2020187987A1 (en) * 2019-03-19 2020-09-24 Koninklijke Philips N.V. Population-level gaussian processes for clinical time series forecasting
WO2020245727A1 (en) * 2019-06-02 2020-12-10 Predicta Med Analytics Ltd. A method of evaluating autoimmune disease risk and treatment selection
CN114330566A (zh) * 2021-12-30 2022-04-12 中山大学 一种脓毒症治疗策略的学习方法及装置
CN114388095A (zh) * 2021-12-22 2022-04-22 中山大学 脓毒症治疗策略优化方法、系统、计算机设备和存储介质
WO2022104412A1 (en) * 2020-11-18 2022-05-27 Epi-Minder Pty Ltd Methods and systems for determination of treatment therapeutic window, detection, prediction, and classification of neuroelectrical, cardiac and/or pulmonary events, and optimization of treatment according to the same

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170193176A1 (en) * 2015-12-30 2017-07-06 Insilico Medicine, Inc. System, Method, and Software for Improved Drug Efficacy and Safety in a Patient
KR101886374B1 (ko) * 2017-08-16 2018-08-07 재단법인 아산사회복지재단 딥러닝 기반의 패혈증 조기 감지방법 및 프로그램
WO2020043922A1 (en) * 2018-08-31 2020-03-05 Novo Nordisk A/S Retrospective horizon based insulin dose prediction
WO2020187987A1 (en) * 2019-03-19 2020-09-24 Koninklijke Philips N.V. Population-level gaussian processes for clinical time series forecasting
WO2020245727A1 (en) * 2019-06-02 2020-12-10 Predicta Med Analytics Ltd. A method of evaluating autoimmune disease risk and treatment selection
CN111696678A (zh) * 2020-06-15 2020-09-22 中南大学 一种基于深度学习的用药决策方法和系统
WO2022104412A1 (en) * 2020-11-18 2022-05-27 Epi-Minder Pty Ltd Methods and systems for determination of treatment therapeutic window, detection, prediction, and classification of neuroelectrical, cardiac and/or pulmonary events, and optimization of treatment according to the same
CN114388095A (zh) * 2021-12-22 2022-04-22 中山大学 脓毒症治疗策略优化方法、系统、计算机设备和存储介质
CN114330566A (zh) * 2021-12-30 2022-04-12 中山大学 一种脓毒症治疗策略的学习方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KUBICEK, J.1等: "《Modeling and objectification of blood vessel calcification with using of multiregional segmentation》", 《VIETNAM JOURNAL OF COMPUTER SCIENCE》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116453706A (zh) * 2023-06-14 2023-07-18 之江实验室 一种基于强化学习的血液透析方案制定方法及系统
CN116453706B (zh) * 2023-06-14 2023-09-08 之江实验室 一种基于强化学习的血液透析方案制定方法及系统
CN116779186A (zh) * 2023-08-23 2023-09-19 寿光市人民医院 一种用于心律失常患者的环磷腺苷用药数据信息处理系统
CN116779186B (zh) * 2023-08-23 2023-12-12 寿光市人民医院 一种用于心律失常患者的环磷腺苷用药数据信息处理系统
CN117275661A (zh) * 2023-11-23 2023-12-22 太原理工大学 一种基于深度强化学习的肺癌患者用药预测方法及装置
CN117275661B (zh) * 2023-11-23 2024-02-09 太原理工大学 一种基于深度强化学习的肺癌患者用药预测方法及装置

Also Published As

Publication number Publication date
CN115050451B (zh) 2022-11-04

Similar Documents

Publication Publication Date Title
CN115050451B (zh) 败血症临床用药方案自动生成系统
Ma et al. Risk prediction on electronic health records with prior medical knowledge
CN106934235B (zh) 一种基于迁移学习的疾病领域间病人相似性度量迁移系统
CN106339593B (zh) 基于医疗数据建模的川崎病分类预测方法
US20190378619A1 (en) Using machine learning to predict health conditions
Huddar et al. Predicting complications in critical care using heterogeneous clinical data
CN107194138B (zh) 一种基于体检数据建模的空腹血糖预测方法
García-Gallo et al. A machine learning-based model for 1-year mortality prediction in patients admitted to an Intensive Care Unit with a diagnosis of sepsis
CN110246577B (zh) 一种基于人工智能辅助妊娠期糖尿病遗传风险预测的方法
CN116631562B (zh) 一种生成电子病历出院记录的方法、装置及电子设备
CN114023449A (zh) 基于深度自编码器的糖尿病风险预警方法与系统
WO2021073255A1 (zh) 基于时序聚类的用药提醒方法及相关设备
CN113128654B (zh) 一种用于冠心病预诊断中的改进型随机森林模型及其预诊断系统
CN113057587A (zh) 一种病症预警方法、装置、电子设备及存储介质
Yang Multi-series Time-aware Sequence Partitioning for Disease Progression Modeling.
Duan et al. Early prediction of sepsis using double fusion of deep features and handcrafted features
CN109119155A (zh) 基于深度学习的icu死亡危险性评估系统
CN114098638B (zh) 一种可解释的动态疾病严重程度预测方法
CN115602317A (zh) 一种基于深度学习的七点血糖预测模型构建方法及装置
CN114613465A (zh) 一种脑卒中患病风险预测和个性化治疗推荐方法及系统
CN113057588A (zh) 一种病症预警方法、装置、设备及介质
Pappy et al. Predicting High-Flow Nasal Cannula Failure in an ICU Using a Recurrent Neural Network with Transfer Learning and Input Data Perseveration: A Retrospective Analysis
CN115762698B (zh) 一种医疗慢病检查报告数据提取方法及系统
US20100121655A1 (en) Patient data record and user interface
Chen et al. Preoperative risk prediction of heart failure with numerical and textual attributes

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