CN113806923B - 一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置 - Google Patents
一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置 Download PDFInfo
- Publication number
- CN113806923B CN113806923B CN202110991468.1A CN202110991468A CN113806923B CN 113806923 B CN113806923 B CN 113806923B CN 202110991468 A CN202110991468 A CN 202110991468A CN 113806923 B CN113806923 B CN 113806923B
- Authority
- CN
- China
- Prior art keywords
- super
- model
- parameter
- pharmacokinetic
- grid
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005457 optimization Methods 0.000 claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000010801 machine learning Methods 0.000 claims abstract description 10
- 238000002790 cross-validation Methods 0.000 claims abstract description 8
- 230000007246 mechanism Effects 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 42
- 238000012795 verification Methods 0.000 claims description 18
- 238000012549 training Methods 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 5
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 239000003814 drug Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 229940079593 drug Drugs 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 230000003285 pharmacodynamic effect Effects 0.000 description 3
- 238000012827 research and development Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 235000000332 black box Nutrition 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009509 drug development Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 229940121657 clinical drug Drugs 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009510 drug design Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000029142 excretion Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229940042040 innovative drug Drugs 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 239000002547 new drug Substances 0.000 description 1
- 230000000144 pharmacologic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/10—ICT 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Primary Health Care (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Medicinal Chemistry (AREA)
- Public Health (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Epidemiology (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明揭示了一种基于nlmixr包的药动学‑药效学模型超参数自动学习方法及装置,该方法包括以下步骤:S1:基于nlmixr软件包构建药动学‑药效学模型;S2:确定超参数空间;S3:结合机器学习算法来获得候选超参数集合,S4:交叉验证机制;S5:包括已取得最佳结果的超参数的组合评分及提供优化过程期间观察到的最好的评分。本技术方案首先通过nlmixr软件包进行药动‑药效学模型的构建,接着确定模型的超参数空间,并进行药动学‑药效学超参数的初始估计,再结合机器学习相关算法实现初始超参数的自动调优。
Description
技术领域
本发明涉及一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置,可用于人工智能药物设计技术领域。
背景技术
药动学-药效学结合模型反映的是药物与机体之间的双向相互作用,其中,机体对药物的作用可用药动学模型表述,包括吸收、分布、代谢和***四个环节,在模型中用药物浓度随时间的变化进行表述。药物对机体的作用反映在药效学模型中,描述了效应随着浓度而变化的动力学过程。据文献报道,进入临床试验约有40%的候选化合物是由于药动学-药效学方面的原因而淘汰,这足以说明药动学-药效学研究在创新药开发研究中的作用。
nlmixr是一个可用来构建药动学-药效学方式模型、传统的隔间药动学模型以及其他更复杂模型的软件包。由于其免费开源,操作简单,功能强大,已经逐渐成为国内外药理学研究中使用最广泛的药动-药效学软件之一。
使用nlmixr软件包对药动-药效学模型深入研究一方面可解决研究人员手工实验的繁琐性,加速新药研发的进程,提高药物开发决策效率,另一方面为临床用药的安全性和有效性提供了更为科学的理论依据,但在构建药动学-药效学模型时,涉及到模型超参数多且复杂,通过人工调参的方法门槛高,而且很难训练出好的模型。
自动机器学习(AutoML)旨在通过让一些通用步骤(如数据预处理、模型选择和调整超参数)自动化,来简化机器学习中生成模型的过程。AutoML是指尽量不通过人来设定超参数,而是使用某种学习机制,来调节这些超参数。超参数与一般模型参数不同,超参数是在训练前提前设置的。超参数优化最常见的类型是黑盒优化(black-box functionoptimization),就是将决策网络当作是一个黑盒来进行优化,仅关心输入和输出,而忽略其内部机制。找到一组超参数,这些超参数返回一个优化模型,该模型减少了预定义的损失函数,进而提高了给定独立数据的预测或者分类精度。
目前需要本领域技术人员迫切解决的一个技术问题就是:如何能够有效的设计出一种新的基于nlmixr软件包的药动学-药效学模型超参数自动学习方法及装置。
发明内容
本发明的目的就是为了解决现有技术中存在的上述问题,提出一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置。
本发明的目的将通过以下技术方案得以实现:一种基于nlmixr包的药动学-药效学模型超参数自动学习方法,该方法包括以下步骤:
S1:基于nlmixr软件包构建药动学-药效学模型;
S2:确定超参数空间;
S3:结合机器学习算法来获得候选超参数集合;
S4:交叉验证机制;
S5:包括已取得最佳结果的超参数的组合评分及提供优化过程期间观察到的最好的评分。
优选地,所述S1步骤又包括以下步骤:
S10:指定特定算法:在通过nlmixr包构建药动学-药效学模型时可使用ODEs模型或solved system模型;
S11:构建模型:药动学-药效学模型包括ini块和model块,其中ini块指定初始条件,包括初始估计,以及支持它们的算法的边界;model模型块用来指定模型。
优选地,所述S2步骤又包括以下步骤:
S20:输入超参数包括一个超参数空间上的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值;
S21:根据每个超参数的网格数(N)和超参数个数(p)确定网格点总数(n);
S22:将超参数空间划分为若干个网格点,在所有网格点中,每个网格点可以与下一个网格点按步长分开;
S23:根据步长确定网格点的坐标。
优选地,在所述S20步骤中,获取基于药动-药效学的超参数,输入超参数包括一个超参数空间上的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值,网格的数量N等于4,超参数的数量p等于2;
S21:根据每个超参数的网格数(N)和超参数个数(p)确定网格点总数(n);
超参数1的范围为[0,5],即下界为“0”,上界为“5”;超参数2的输入范围为[0,10],即下界为“0”,下界为“10”,网格点总数(n)确定为“n=Np”,n=42=16;
S23:计算每个网格点的坐标值及步长,包括以下过程:
网格搜索技术将超参数空间划分为若干个网格点,在所有网格点中,每个网格点可以与下一个网格点按步长分开,如公式(2)所示:
式中,UBi和LBi分别为超参数“i”的上界和下界;
此外,网格点总数中的每个网格点可以用一组坐标表示,
其中,″ri″=0,1,2,…,(N-1),″i″=0,1,2,…,(p-1),网格的坐标形式为(x,y);
为每个网格点计算的坐标如下式所示,(2,4)是坐标为(1,1)的格点:
优选地,所述S3步骤又包括以下步骤:
S31:将目标函数值相互比较,以识别具有最小目标函数值的网格点,目标函数值为当前目标函数值;
S32:通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合。
优选地,在所述S31步骤中,将目标函数值(即当前目标函数值)相互比较,以识别具有最小目标函数值的网格点;
为因变量的观测值,/>为因变量的预测值,所识别的目标函数值最小的网格点即为药动学-药效学超参数。
优选地,所述S4步骤又包括以下步骤:
S41:交叉验证对训练集等分成N份,N为用户指定的值;
S42:将其中一份作为验证集,其余N-1份作为训练集,经过N次测试,每次都更换不同的验证集,得到N个模型结果,取最优结果;
S43:使用最优超参数重新训练模型,实现自动调节超参数的过程。
优选地,在所述S41步骤中,指定N为10,即为10折交叉验证。
本发明还揭示了一种基于nlmixr包的药动学-药效学模型超参数自动学习装置,该装置包括:基于nlmixr包的药动学-药效学模型构建模块,用来生成药动-药效学模型及提供超参数优化的数据集;
超参数空间生成模块,用来接收生成超参数空间中的各个超参数及构建超参数空间;
超参数自动优化模块,用来实现候选超参数的自动优化;
药动学-药效学模型超参数自动学习的性能评价模块,用来表示所选超参数的得分情况。
优选地,所述基于nlmixr包的药动学-药效学模型构建模块包括:ini模块:指定初始条件,包括初始估计,以及支持它们的算法的边界;model模块:用来构建模型,model模块选择使用残差模型,加性残差模型或者比例残差模型;
超参数空间生成模块包括:初始模块,用于接收药动学-药效学超参数的超参数上下界值,超参数个数,每个超参数的网格数及网格总数信息;构建模块,根据输入超参数将超参数空间划分为若干个网格点,每个网格点可以与下一个网格点按步长分开;根据步长确定网格点的坐标;
超参数自动优化模块包括:搜索模块,通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合;重新训练模型,使用最优超参数重新训练模型;
药动学-药效学模型超参数自动学习的性能评价模块包括:best_parameters模块:描述了已取得最佳结果的超参数的组合;best_score模块:提供优化过程期间观察到的最好的评分。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:本技术方案首先通过nlmixr软件包进行药动-药效学模型的构建,接着确定模型的超参数空间,并进行药动学-药效学超参数的初始估计,再结合机器学习相关算法实现初始超参数的自动调优,降低了人工调参的门槛,有助于构建更优模型,以加快药物研发进程。
附图说明
图1为本发明的一种基于nlmixr包的药动学-药效学结合模型的示意图。
图2为本发明的一种基于nlmixr包的药动学-药效学结合模型的示意图。
图3为本发明的一种基于nlmixr包的药动学-药效学结合模型的示意图。
图4为本发明的一种基于nlmixr包的药动学-药效学结合模型的示意图。
图5为本发明的一种药动学-药效学模型超参数自动学习方法的结构示意图。
图6为本发明中一种药动学-药效学模型超参数自动学习装置的实现示意图。
图7是本发明中超参数空间中网格点以坐标表示的格点示意图。
图8是本发明中超参数自动优化模块303的结构示意图。
图9是本发明中药动学-药效学模型超参数自动学习的性能评价模块304的结构示意图。
具体实施方式
本发明的目的、优点和特点,将通过下面优选实施例的非限制性说明进行图示和解释。这些实施例仅是应用本发明技术方案的典型范例,凡采取等同替换或者等效变换而形成的技术方案,均落在本发明要求保护的范围之内。
本发明揭示了一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置,如图1、图2和图3所示,首先通过nlmixr软件包进行药动-药效学模型的构建,接着确定模型的超参数空间,并进行药动学-药效学超参数的初始估计,再结合机器学习相关算法实现初始超参数的自动调优。
一种基于nlmixr软件包的药动学-药效学模型超参数自动学习方法,该方法包括以下步骤:
S1:基于nlmixr软件包构建药动学-药效学模型;
所述S1步骤又包括以下步骤:
S10:指定特定算法:在通过nlmixr包构建药动学-药效学模型时可使用ODEs模型或者solved system模型;
S11:构建模型:药动学-药效学模型包括ini块和model块,其中ini块指定初始条件,包括初始估计,以及支持它们的算法的边界;model模型块用来指定模型,类似于NONMEM中的$PK、$PRED和$ERROR块。
S2:确定超参数空间;
所述S2步骤又包括以下步骤:
S20:输入超参数包括一个超参数空间E的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值;
S21:根据每个超参数的网格数(N)和超参数个数(p)确定网格点总数(n);
S22:将超参数空间划分为若干个网格点,在所有网格点中,每个网格点可以与下一个网格点按步长分开;
S23:根据步长确定网格点的坐标。
S3:结合机器学习算法来获得候选超参数集合,
所述S3步骤又包括以下步骤:
S31:将目标函数值相互比较,以识别具有最小目标函数值的网格点,目标函数值为当前目标函数值;
S32:通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合。
S4:交叉验证机制;
所述S4步骤又包括以下步骤:
S41:交叉验证对训练集等分成N份,N为用户指定的值,比如可指定为10,即为10折交叉验证;
S42:将其中一份作为验证集,其余N-1份作为训练集,经过N次测试,每次都更换不同的验证集,得到N个模型结果,取最优结果;
S43:使用最优超参数重新训练模型,实现自动调节超参数的过程。
S5:包括已取得最佳结果的超参数的组合评分及提供优化过程期间观察到的最好的评分。
一种基于nlmixr包的药动学-药效学模型超参数自动学习装置,该装置包括:基于nlmixr包的药动学-药效学模型构建模块301,用来生成药动-药效学模型及提供超参数优化的数据集;超参数空间生成模块302,用来接收生成超参数空间中的各个超参数及构建超参数空间。
超参数自动优化模块303,用来实现候选超参数的自动优化;药动学-药效学模型超参数自动学习的性能评价模块304,用来表示所选超参数的得分情况。
其中,基于nlmixr包的药动学-药效学模型构建模块301具体包括:ini模块:指定初始条件,包括初始估计,以及支持它们的算法的边界;model模块:用来构建模型,此步骤可选择使用残差模型,比如加性残差模型或者比例残差模型。
其中,超参数空间生成模块302具体包括:初始模块,用于接收药动学-药效学超参数的超参数上下界值,超参数个数,每个超参数的网格数及网格总数信息。
构建模块,如图4所示,根据输入超参数将超参数空间划分为若干个网格点,每个网格点可以与下一个网格点按步长分开;根据步长确定网格点的坐标。其中,超参数自动优化模块具体包括:搜索模块,通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合;重新训练模型,使用最优超参数重新训练模型。
如图5所示,具体包括:交叉验证对训练集等分成N份,N为用户指定的值,比如可指定为10,即为10折交叉验证;将其中一份作为验证集,其余N-1份作为训练集,经过N次测试,每次都更换不同的验证集,得到N个模型结果,取最优结果;使用最优超参数重新训练模型。
其中,药动学-药效学模型超参数自动学习的性能评价模块,如图6所示,具体包括:best_parameters模块601:描述了已取得最佳结果的超参数的组合;best_score模块602:提供优化过程期间观察到的最好的评分。
实施例:
如图1和图2所示,该药动学-药效学模型超参数自动学习的方法,包括以下步骤:
S1:构建基于nlmixr包的药动学-药效学模型,
药动学-药效学模型可选择房室模型或单房室模型或多房室模型、代数模型和微分方程模型,如图1所示是一个二房室的通用结构,对应方程(1)列出了常见的二房室模型结构。
C=Ae-αt+Be-βt (1)
其中,C为药物在人/动物体内的血药浓度,t为时间,α、β在二房室模型中分别为分布速率常数和消除速率常数,A、B为α、β相延伸线在纵轴的截距,不同模型可进行相应修改。
图2、图3和图4描述了药动-药效学结合模型以及药动学模型,药效学模型的曲线图。药动学-药效学结合模型反映的是药物与机体之间的双向相互作用。其中,机体对药物的作用可用药动学模型表述,在模型中用药物浓度随时间的变化进行表述;药物对机体的作用反映在药效学模型中,描述了效应随着浓度而变化的动力学过程。
所述S1步骤又包括以下步骤:
S10:指定特定算法:在使用nlmixr包构建药动学-药效学模型时可使用ODEs模型或者solved system模型;
S11:构建模型:
如图2所示,使用nlmixr软件包药动学-药效学模型包括ini块和model块,其中ini块指定初始条件,包括初始估计,以及支持它们的算法的边界;model模型块用来指定模型,类似于NONMEM中的$PK、$PRED和$ERROR块;
S2:确定超参数空间;
所述S2步骤又包括以下步骤:
S20:输入超参数包括一个超参数空间上的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值;
获取基于药动-药效学的超参数,输入超参数包括一一个超参数空间上的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值。本实施例中网格的数量即N等于4,超参数的数量即p等于2。
S21:根据每个超参数的网格数(N)和超参数个数(p)确定网格点总数(n);
此外,超参数1的范围为[0,5],即下界为“0”,上界为“5”;超参数2的输入范围为[0,10],即下界为“0”,下界为“10”,可将网格点总数(n)确定为“n=Np”,即n=42=16。
S23:计算每个网格点的坐标值及步长,包括以下过程:
网格搜索技术将超参数空间划分为若干个网格点,在所有网格点中,每个网格点可以与下一个网格点按步长分开,如公式(2)所示:
式中,UBi和LBi分别为超参数“i”的上界和下界。
此外,网格点总数中的每个网格点可以用一组坐标表示,如图2所示,
其中,″ri″=0,1,2,…,(N-1),″i″=0,1,2,…,(p-1),网格的坐标形式为(x,y)。
为每个网格点计算的坐标如下所示,如图7所示,(2,4)是坐标为(1,1)的格点:
S3:结合机器学习算法来获得候选超参数集合,
所述S3步骤又包括以下步骤:
S31:将目标函数值相互比较,以识别具有最小目标函数值的网格点,目标函数值为当前目标函数值;
将目标函数值(即当前目标函数值)相互比较,以识别具有最小目标函数值的网格点。
为因变量的观测值,/>为因变量的预测值,所识别的目标函数值最小的网格点即为药动学-药效学超参数。
S32:通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合。
S4:交叉验证机制;
所述S4步骤又包括以下步骤:
S41:交叉验证对训练集等分成N份,N为用户指定的值,本技术方案中,N指定为5,即为5折交叉验证;
S42:将其中一份作为验证集,其余4份作为训练集,经过5次测试,每次都更换不同的验证集,得到5个模型结果,取最优结果;
S43:使用最优超参数重新训练模型。
S5:评分函数,包括以下过程:best_parameters模块:描述了已取得最佳结果的超参数的组合;best_score模块:提供优化过程期间观察到的最好的评分。
在本技术方案中,基于nlmixr包构建的模型性能与超参数直接相关,这说明准确设定模型的超参数是非常必要的。该方法包括构建药动学-药效学模型;确定超参数空间;通过搜索的方式来获得候选参数集合;交叉验证机制;评分函数。该装置包括基于nlmixr包的药动学-药效学模型构建模块;参数空间生成模块;超参数自动优化模块;药动学-药效学模型超参自动学习的性能评价模块。根据本技术方案提出的一种基于nlmixr包的药动学-药效学模型超参数自动学习的方法,可以简单实现模型超参数自动调优问题,帮助构建更优的药动学-药效学模型,以便加快药物研发进程和节约研发成本。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神和基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内,不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。本发明尚有多种实施方式,凡采用等同变换或者等效变换而形成的所有技术方案,均落在本发明的保护范围之内。
Claims (7)
1.一种基于nlmixr包的药动学-药效学模型超参数自动学习方法,其特征在于:该方法包括以下步骤:
S1:基于nlmixr软件包构建药动学-药效学模型;
S2:确定超参数空间,包括以下步骤:
S20:输入超参数包括一个超参数空间上的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值;其中,获取基于药动-药效学的超参数,输入超参数包括一个超参数空间上的网格数(N)、超参数个数(p)、每个超参数的下界值和上界值,网格的数量N等于4,超参数的数量p等于2;
S21:根据每个超参数的网格数(N)和超参数个数(p)确定网格点总数(n);超参数1的范围为[0,5],即下界为“0”,上界为“5”;超参数2的输入范围为[0,10],即下界为“0”,下界为“10”,网格点总数(n)确定为“n=Np”,n=42=16;
S22:将超参数空间划分为若干个网格点,在所有网格点中,每个网格点可以与下一个网格点按步长分开;
S23:根据步长确定网格点的坐标,包括以下过程:
网格搜索技术将超参数空间划分为若干个网格点,在所有网格点中,每个网格点可以与下一个网格点按步长分开,如公式(2)所示:
式中,UBi和LBi分别为超参数“i”的上界和下界;
此外,网格点总数中的每个网格点可以用一组坐标表示,
其中,″ri″=0,1,2,…,(N-1),″i″=0,1,2,…,(p-1),网格的坐标形式为(x,y);
为每个网格点计算的坐标如下式所示,(2,4)是坐标为(1,1)的格点:
S3:结合机器学习算法来获得候选超参数集合,包括以下步骤:
S31:将目标函数值相互比较,以识别具有最小目标函数值的网格点,目标函数值为当前目标函数值;
S32:通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合;
S4:交叉验证机制;
S5:包括已取得最佳结果的超参数的组合评分及提供优化过程期间观察到的最好的评分。
2.根据权利要求1所述的一种基于nlmixr包的药动学-药效学模型超参数自动学习方法,其特征在于:所述S1步骤又包括以下步骤:
S10:指定特定算法:在通过nlmixr包构建药动学-药效学模型时可使用ODEs模型或solved system模型;
S11:构建模型:药动学-药效学模型包括ini块和model块,其中ini块指定初始条件,包括初始估计,以及支持它们的算法的边界;model模型块用来指定模型。
3.根据权利要求1所述的一种基于nlmixr包的药动学-药效学模型超参数自动学习方法,其特征在于:
在所述S31步骤中,将目标函数值(即当前目标函数值)相互比较,以识别具有最小目标函数值的网格点;
为因变量的观测值,/>为因变量的预测值,所识别的目标函数值最小的网格点即为药动学-药效学超参数。
4.根据权利要求1所述的一种基于nlmixr包的药动学-药效学模型超参数自动学习方法,其特征在于:所述S4步骤又包括以下步骤:
S41:交叉验证对训练集等分成N份,N为用户指定的值;
S42:将其中一份作为验证集,其余N-1份作为训练集,经过N次测试,每次都更换不同的验证集,得到N个模型结果,取最优结果;
S43:使用最优超参数重新训练模型,实现自动调节超参数的过程。
5.根据权利要求4所述的一种基于nlmixr包的药动学-药效学模型超参数自动学习方法,其特征在于:在所述S41步骤中,指定N为10,即为10折交叉验证。
6.一种基于nlmixr包的药动学-药效学模型超参数自动学习装置,用于实现如权利要求1-5任意一项基于nlmixr包的药动学-药效学模型超参数自动学习方法,其特征在于:该装置包括:基于nlmixr包的药动学-药效学模型构建模块(301),用来生成药动-药效学模型及提供超参数优化的数据集;
超参数空间生成模块(302),用来接收生成超参数空间中的各个超参数及构建超参数空间;
超参数自动优化模块(303),用来实现候选超参数的自动优化;
药动学-药效学模型超参数自动学习的性能评价模块(304),用来表示所选超参数的得分情况。
7.根据权利要求6所述的一种基于nlmixr包的药动学-药效学模型超参数自动学习装置,其特征在于:所述基于nlmixr包的药动学-药效学模型构建模块(301)包括:ini模块:指定初始条件,包括初始估计,以及支持它们的算法的边界;model模块:用来构建模型,model模块选择使用残差模型,加性残差模型或者比例残差模型;
超参数空间生成模块(302)包括:初始模块,用于接收药动学-药效学超参数的超参数上下界值,超参数个数,每个超参数的网格数及网格总数信息;构建模块,根据输入超参数将超参数空间划分为若干个网格点,每个网格点可以与下一个网格点按步长分开;根据步长确定网格点的坐标;
超参数自动优化模块(303)包括:搜索模块,通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合;重新训练模型,使用最优超参数重新训练模型;
药动学-药效学模型超参数自动学习的性能评价模块(304)包括:best_parameters模块(601):描述了已取得最佳结果的超参数的组合;best_score模块(602):提供优化过程期间观察到的最好的评分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110991468.1A CN113806923B (zh) | 2021-08-26 | 2021-08-26 | 一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110991468.1A CN113806923B (zh) | 2021-08-26 | 2021-08-26 | 一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113806923A CN113806923A (zh) | 2021-12-17 |
CN113806923B true CN113806923B (zh) | 2023-12-22 |
Family
ID=78941853
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110991468.1A Active CN113806923B (zh) | 2021-08-26 | 2021-08-26 | 一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113806923B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114331184B (zh) * | 2021-12-31 | 2022-11-18 | 北京森科医药有限公司 | 一种核素药物实验规划方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109657805A (zh) * | 2018-12-07 | 2019-04-19 | 泰康保险集团股份有限公司 | 超参数确定方法、装置、电子设备及计算机可读介质 |
WO2020090826A1 (ja) * | 2018-10-30 | 2020-05-07 | 日本電気株式会社 | 解析装置、機械学習装置、解析システム、解析方法および記録媒体 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11705226B2 (en) * | 2019-09-19 | 2023-07-18 | Tempus Labs, Inc. | Data based cancer research and treatment systems and methods |
-
2021
- 2021-08-26 CN CN202110991468.1A patent/CN113806923B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020090826A1 (ja) * | 2018-10-30 | 2020-05-07 | 日本電気株式会社 | 解析装置、機械学習装置、解析システム、解析方法および記録媒体 |
CN109657805A (zh) * | 2018-12-07 | 2019-04-19 | 泰康保险集团股份有限公司 | 超参数确定方法、装置、电子设备及计算机可读介质 |
Non-Patent Citations (2)
Title |
---|
基于强化学习的模型选择和超参数优化;吴佳;陈森朋;陈修云;周瑞;;电子科技大学学报(02);全文 * |
基于贝叶斯网的抗肺结核诊疗数据分析;陆维嘉;;计算机与数字工程(12);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113806923A (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104200087B (zh) | 用于机器学习的参数寻优及特征调优的方法及*** | |
Liu et al. | Fast and adaptive sparse precision matrix estimation in high dimensions | |
JP2019535057A (ja) | タンパク質結合部位予測の方法、装置、設備及び記憶媒体 | |
CN103473482B (zh) | 基于差分进化和构象空间退火的蛋白质三维结构预测方法 | |
KR101544457B1 (ko) | 최적 설계 파라미터 탐색을 위한 최적화 방법 | |
JP6694447B2 (ja) | ビッグデータの計算方法及びシステム、プログラムならびに記録媒体 | |
CN113806923B (zh) | 一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置 | |
CN106796533A (zh) | 自适应地选择执行模式的***和方法 | |
Chen et al. | Interactive multiobjective evolutionary algorithm based on decomposition and compression | |
CN109360596B (zh) | 一种基于差分进化局部扰动的蛋白质构象空间优化方法 | |
CN108491686B (zh) | 一种基于双向XGBoost的基因调控网络构建方法 | |
Urbanowicz et al. | An extended michigan-style learning classifier system for flexible supervised learning, classification, and data mining | |
Soto et al. | Boosting autonomous search for CSPs via skylines | |
Mandal et al. | Learning fitness functions for machine programming | |
CN114627980A (zh) | 一种化学逆合成分析方法及*** | |
Stenger et al. | Benchmark of bayesian optimization and metaheuristics for control engineering tuning problems with crash constraints | |
CN111820922B (zh) | 一种面向幼儿计算思维的测评方法 | |
Stober et al. | Learning geometry from sensorimotor experience | |
van Kesteren et al. | Structural equation models as computation graphs | |
CN106815858A (zh) | 一种运动目标提取方法及装置 | |
Chen et al. | Efficient model evaluation in the search-based approach to latent structure discovery | |
Benmouna et al. | New method for Bayesian network learning | |
Hironaka et al. | Analysis and improvements of the Pareto optimal solution visualization method using the self-organizing maps | |
Botzheim et al. | Fuzzy rule base model identification by bacterial memetic algorithms | |
CN105022798A (zh) | 一种基于预测关系的离散贝叶斯网络分类数据挖掘方法 |
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 |