一种基于nlmixr包的药动学-药效学模型超参数自动学习方法及装置与流程

一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习方法及装置
技术领域
1.本发明涉及一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习方法及装置，可用于人工智能药物设计技术领域。


背景技术：

2.药动学
‑
药效学结合模型反映的是药物与机体之间的双向相互作用，其中，机体对药物的作用可用药动学模型表述，包括吸收、分布、代谢和排泄四个环节，在模型中用药物浓度随时间的变化进行表述。药物对机体的作用反映在药效学模型中，描述了效应随着浓度而变化的动力学过程。据文献报道，进入临床试验约有40％的候选化合物是由于药动学
‑
药效学方面的原因而淘汰，这足以说明药动学
‑
药效学研究在创新药开发研究中的作用。
3.nlmixr是一个可用来构建药动学
‑
药效学方式模型、传统的隔间药动学模型以及其他更复杂模型的软件包。由于其免费开源，操作简单，功能强大，已经逐渐成为国内外药理学研究中使用最广泛的药动
‑
药效学软件之一。
4.使用nlmixr软件包对药动
‑
药效学模型深入研究一方面可解决研究人员手工实验的繁琐性，加速新药研发的进程，提高药物开发决策效率，另一方面为临床用药的安全性和有效性提供了更为科学的理论依据，但在构建药动学
‑
药效学模型时，涉及到模型超参数多且复杂，通过人工调参的方法门槛高，而且很难训练出好的模型。
5.自动机器学习(automl)旨在通过让一些通用步骤(如数据预处理、模型选择和调整超参数)自动化，来简化机器学习中生成模型的过程。automl是指尽量不通过人来设定超参数，而是使用某种学习机制，来调节这些超参数。超参数与一般模型参数不同，超参数是在训练前提前设置的。超参数优化最常见的类型是黑盒优化(black
‑
box function optimization)，就是将决策网络当作是一个黑盒来进行优化，仅关心输入和输出，而忽略其内部机制。找到一组超参数，这些超参数返回一个优化模型，该模型减少了预定义的损失函数，进而提高了给定独立数据的预测或者分类精度。
6.目前需要本领域技术人员迫切解决的一个技术问题就是：如何能够有效的设计出一种新的基于nlmixr软件包的药动学
‑
药效学模型超参数自动学习方法及装置。


技术实现要素：

7.本发明的目的就是为了解决现有技术中存在的上述问题，提出一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习方法及装置。
8.本发明的目的将通过以下技术方案得以实现：一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习方法，该方法包括以下步骤：
9.s1：基于nlmixr软件包构建药动学
‑
药效学模型；
10.s2：确定超参数空间；
11.s3：结合机器学习算法来获得候选超参数集合；
12.s4：交叉验证机制；
13.s5：包括已取得最佳结果的超参数的组合评分及提供优化过程期间观察到的最好的评分。
14.优选地，所述s1步骤又包括以下步骤：
15.s10：指定特定算法：在通过nlmixr包构建药动学
‑
药效学模型时可使用odes模型或solved system模型；
16.s11：构建模型：药动学
‑
药效学模型包括ini块和model块，其中ini块指定初始条件，包括初始估计，以及支持它们的算法的边界；model模型块用来指定模型。
17.优选地，所述s2步骤又包括以下步骤：
18.s20：输入超参数包括一个超参数空间上的网格数(n)、超参数个数(p)、每个超参数的下界值和上界值；
19.s21：根据每个超参数的网格数(n)和超参数个数(p)确定网格点总数(n)；
20.s22：将超参数空间划分为若干个网格点，在所有网格点中，每个网格点可以与下一个网格点按步长分开；
21.s23：根据步长确定网格点的坐标。
22.优选地，在所述s20步骤中，获取基于药动
‑
药效学的超参数，输入超参数包括一个超参数空间上的网格数(n)、超参数个数(p)、每个超参数的下界值和上界值，网格的数量n等于4，超参数的数量p等于2；
23.s21：根据每个超参数的网格数(n)和超参数个数(p)确定网格点总数(n)；
24.超参数1的范围为[0，5]，即下界为“0”，上界为“5”；超参数2的输入范围为[0，10]，即下界为“0”，下界为“10”，网格点总数(n)确定为“n＝n
p”，n＝42＝16；
[0025]
s23：计算每个网格点的坐标值及步长，包括以下过程：
[0026]
网格搜索技术将超参数空间划分为若干个网格点，在所有网格点中，每个网格点可以与下一个网格点按步长分开，如公式(2)所示：
[0027][0028]
式中，ub
i
和lb
i
分别为超参数“i”的上界和下界；
[0029]
此外，网格点总数中的每个网格点可以用一组坐标表示，
[0030][0031]
其中，
″
r
i
″
＝0，1，2，
…
，(n
‑
1)，
″
i
″
＝0，1，2，
…
，(p
‑
1)，网格的坐标形式为(x，y)；
[0032]
为每个网格点计算的坐标如下式所示，(2，4)是坐标为(1，1)的格点：
[0033][0034]
优选地，所述s3步骤又包括以下步骤：
[0035]
s31：将目标函数值相互比较，以识别具有最小目标函数值的网格点，目标函数值为当前目标函数值；
[0036]
s32：通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合。
[0037]
优选地，在所述s31步骤中，将目标函数值(即当前目标函数值)相互比较，以识别具有最小目标函数值的网格点；
[0038][0039]
为因变量的观测值，为因变量的预测值，所识别的目标函数值最小的网格点即为药动学
‑
药效学超参数。
[0040]
优选地，所述s4步骤又包括以下步骤：
[0041]
s41：交叉验证对训练集等分成n份，n为用户指定的值；
[0042]
s42：将其中一份作为验证集，其余n
‑
1份作为训练集，经过n次测试，每次都更换不同的验证集，得到n个模型结果，取最优结果；
[0043]
s43：使用最优超参数重新训练模型，实现自动调节超参数的过程。
[0044]
优选地，在所述s41步骤中，指定n为10，即为10折交叉验证。
[0045]
本发明还揭示了一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习装置，该装置包括：基于nlmixr包的药动学
‑
药效学模型构建模块，用来生成药动
‑
药效学模型及提供超参数优化的数据集；
[0046]
超参数空间生成模块，用来接收生成超参数空间中的各个超参数及构建超参数空间；
[0047]
超参数自动优化模块，用来实现候选超参数的自动优化；
[0048]
药动学
‑
药效学模型超参数自动学习的性能评价模块，用来表示所选超参数的得分情况。
[0049]
优选地，所述基于nlmixr包的药动学
‑
药效学模型构建模块包括：ini模块：指定初始条件，包括初始估计，以及支持它们的算法的边界；model模块：用来构建模型，model模块选择使用残差模型，加性残差模型或者比例残差模型；
[0050]
超参数空间生成模块包括：初始模块，用于接收药动学
‑
药效学超参数的超参数上下界值，超参数个数，每个超参数的网格数及网格总数信息；构建模块，根据输入超参数将超参数空间划分为若干个网格点，每个网格点可以与下一个网格点按步长分开；根据步长确定网格点的坐标；
[0051]
超参数自动优化模块包括：搜索模块，通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合；重新训练模型，使用最优超参数重新训练模型；
[0052]
药动学
‑
药效学模型超参数自动学习的性能评价模块包括：best_parameters模块：描述了已取得最佳结果的超参数的组合；best_score模块：提供优化过程期间观察到的最好的评分。
[0053]
本发明采用以上技术方案与现有技术相比，具有以下技术效果：本技术方案首先通过nlmixr软件包进行药动
‑
药效学模型的构建，接着确定模型的超参数空间，并进行药动学
‑
药效学超参数的初始估计，再结合机器学习相关算法实现初始超参数的自动调优，降低了人工调参的门槛，有助于构建更优模型，以加快药物研发进程。
附图说明
[0054]
图1为本发明的一种基于nlmixr包的药动学
‑
药效学结合模型的示意图。
[0055]
图2为本发明的一种基于nlmixr包的药动学
‑
药效学结合模型的示意图。
[0056]
图3为本发明的一种基于nlmixr包的药动学
‑
药效学结合模型的示意图。
[0057]
图4为本发明的一种基于nlmixr包的药动学
‑
药效学结合模型的示意图。
[0058]
图5为本发明的一种药动学
‑
药效学模型超参数自动学习方法的结构示意图。
[0059]
图6为本发明中一种药动学
‑
药效学模型超参数自动学习装置的实现示意图。
[0060]
图7是本发明中超参数空间中网格点以坐标表示的格点示意图。
[0061]
图8是本发明中超参数自动优化模块303的结构示意图。
[0062]
图9是本发明中药动学
‑
药效学模型超参数自动学习的性能评价模块304的结构示意图。
具体实施方式
[0063]
本发明的目的、优点和特点，将通过下面优选实施例的非限制性说明进行图示和解释。这些实施例仅是应用本发明技术方案的典型范例，凡采取等同替换或者等效变换而形成的技术方案，均落在本发明要求保护的范围之内。
[0064]
本发明揭示了一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习方法及装置，如图1、图2和图3所示，首先通过nlmixr软件包进行药动
‑
药效学模型的构建，接着确定模型的超参数空间，并进行药动学
‑
药效学超参数的初始估计，再结合机器学习相关算法实现初始超参数的自动调优。
[0065]
一种基于nlmixr软件包的药动学
‑
药效学模型超参数自动学习方法，该方法包括以下步骤：
[0066]
s1：基于nlmixr软件包构建药动学
‑
药效学模型；
[0067]
所述s1步骤又包括以下步骤：
[0068]
s10：指定特定算法：在通过nlmixr包构建药动学
‑
药效学模型时可使用odes模型或者solved system模型；
[0069]
s11：构建模型：药动学
‑
药效学模型包括ini块和model块，其中ini块指定初始条件，包括初始估计，以及支持它们的算法的边界；model模型块用来指定模型，类似于nonmem中的$pk、$pred和$error块。
[0070]
s2：确定超参数空间；
[0071]
所述s2步骤又包括以下步骤：
[0072]
s20：输入超参数包括一个超参数空间e的网格数(n)、超参数个数(p)、每个超参数的下界值和上界值；
[0073]
s21：根据每个超参数的网格数(n)和超参数个数(p)确定网格点总数(n)；
[0074]
s22：将超参数空间划分为若干个网格点，在所有网格点中，每个网格点可以与下一个网格点按步长分开；
[0075]
s23：根据步长确定网格点的坐标。
[0076]
s3：结合机器学习算法来获得候选超参数集合，
[0077]
所述s3步骤又包括以下步骤：
[0078]
s31：将目标函数值相互比较，以识别具有最小目标函数值的网格点，目标函数值为当前目标函数值；
[0079]
s32：通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合。
[0080]
s4：交叉验证机制；
[0081]
所述s4步骤又包括以下步骤：
[0082]
s41：交叉验证对训练集等分成n份，n为用户指定的值，比如可指定为10，即为10折交叉验证；
[0083]
s42：将其中一份作为验证集，其余n
‑
1份作为训练集，经过n次测试，每次都更换不同的验证集，得到n个模型结果，取最优结果；
[0084]
s43：使用最优超参数重新训练模型，实现自动调节超参数的过程。
[0085]
s5：包括已取得最佳结果的超参数的组合评分及提供优化过程期间观察到的最好的评分。
[0086]
一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习装置，该装置包括：基于nlmixr包的药动学
‑
药效学模型构建模块301，用来生成药动
‑
药效学模型及提供超参数优化的数据集；超参数空间生成模块302，用来接收生成超参数空间中的各个超参数及构建超参数空间。
[0087]
超参数自动优化模块303，用来实现候选超参数的自动优化；药动学
‑
药效学模型超参数自动学习的性能评价模块304，用来表示所选超参数的得分情况。
[0088]
其中，基于nlmixr包的药动学
‑
药效学模型构建模块301具体包括：ini模块：指定初始条件，包括初始估计，以及支持它们的算法的边界；model模块：用来构建模型，此步骤可选择使用残差模型，比如加性残差模型或者比例残差模型。
[0089]
其中，超参数空间生成模块302具体包括：初始模块，用于接收药动学
‑
药效学超参数的超参数上下界值，超参数个数，每个超参数的网格数及网格总数信息。
[0090]
构建模块，如图4所示，根据输入超参数将超参数空间划分为若干个网格点，每个网格点可以与下一个网格点按步长分开；根据步长确定网格点的坐标。其中，超参数自动优化模块具体包括：搜索模块，通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合；重新训练模型，使用最优超参数重新训练模型。
[0091]
如图5所示，具体包括：交叉验证对训练集等分成n份，n为用户指定的值，比如可指定为10，即为10折交叉验证；将其中一份作为验证集，其余n
‑
1份作为训练集，经过n次测试，每次都更换不同的验证集，得到n个模型结果，取最优结果；使用最优超参数重新训练模型。
[0092]
其中，药动学
‑
药效学模型超参数自动学习的性能评价模块，如图6所示，具体包括：best_parameters模块601：描述了已取得最佳结果的超参数的组合；best_score模块602：提供优化过程期间观察到的最好的评分。
[0093]
实施例：
[0094]
如图1和图2所示，该药动学
‑
药效学模型超参数自动学习的方法，包括以下步骤：
[0095]
s1：构建基于nlmixr包的药动学
‑
药效学模型，
[0096]
药动学
‑
药效学模型可选择房室模型或单房室模型或多房室模型、代数模型和微分方程模型，如图1所示是一个二房室的通用结构，对应方程(1)列出了常见的二房室模型
结构。
[0097]
c＝ae
‑
αt
be
‑
βt
ꢀꢀꢀ
(1)
[0098]
其中，c为药物在人/动物体内的血药浓度，t为时间，α、β在二房室模型中分别为分布速率常数和消除速率常数，a、b为α、β相延伸线在纵轴的截距，不同模型可进行相应修改。
[0099]
图2、图3和图4描述了药动
‑
药效学结合模型以及药动学模型，药效学模型的曲线图。药动学
‑
药效学结合模型反映的是药物与机体之间的双向相互作用。其中，机体对药物的作用可用药动学模型表述，在模型中用药物浓度随时间的变化进行表述；药物对机体的作用反映在药效学模型中，描述了效应随着浓度而变化的动力学过程。
[0100]
所述s1步骤又包括以下步骤：
[0101]
s10：指定特定算法：在使用nlmixr包构建药动学
‑
药效学模型时可使用odes模型或者solved system模型；
[0102]
s11：构建模型：
[0103]
如图2所示，使用nlmixr软件包药动学
‑
药效学模型包括ini块和model块，其中ini块指定初始条件，包括初始估计，以及支持它们的算法的边界；model模型块用来指定模型，类似于nonmem中的$pk、$pred和$error块；
[0104]
s2：确定超参数空间；
[0105]
所述s2步骤又包括以下步骤：
[0106]
s20：输入超参数包括一个超参数空间上的网格数(n)、超参数个数(p)、每个超参数的下界值和上界值；
[0107]
获取基于药动
‑
药效学的超参数，输入超参数包括一一个超参数空间上的网格数(n)、超参数个数(p)、每个超参数的下界值和上界值。本实施例中网格的数量即n等于4，超参数的数量即p等于2。
[0108]
s21：根据每个超参数的网格数(n)和超参数个数(p)确定网格点总数(n)；
[0109]
此外，超参数1的范围为[0，5]，即下界为“0”，上界为“5”；超参数2的输入范围为[0，10]，即下界为“0”，下界为“10”，可将网格点总数(n)确定为“n＝n
p”，即n＝42＝16。
[0110]
s23：计算每个网格点的坐标值及步长，包括以下过程：
[0111]
网格搜索技术将超参数空间划分为若干个网格点，在所有网格点中，每个网格点可以与下一个网格点按步长分开，如公式(2)所示：
[0112][0113]
式中，ub
i
和lb
i
分别为超参数“i”的上界和下界。
[0114]
此外，网格点总数中的每个网格点可以用一组坐标表示，如图2所示，
[0115][0116]
其中，
″
r
i
″
＝0，1，2，
…
，(n
‑
1)，
″
i
″
＝0，1，2，
…
，(p
‑
1)，网格的坐标形式为(x，y)。
[0117]
为每个网格点计算的坐标如下所示，如图7所示，(2，4)是坐标为(1，1)的格点：
[0118]
[0119]
s3：结合机器学习算法来获得候选超参数集合，
[0120]
所述s3步骤又包括以下步骤：
[0121]
s31：将目标函数值相互比较，以识别具有最小目标函数值的网格点，目标函数值为当前目标函数值；
[0122]
将目标函数值(即当前目标函数值)相互比较，以识别具有最小目标函数值的网格点。
[0123][0124]
为因变量的观测值，为因变量的预测值，所识别的目标函数值最小的网格点即为药动学
‑
药效学超参数。
[0125]
s32：通过网格搜索技术从所有网格点中选取目标函数值最小的网格点作为候选超参数集合。
[0126]
s4：交叉验证机制；
[0127]
所述s4步骤又包括以下步骤：
[0128]
s41：交叉验证对训练集等分成n份，n为用户指定的值，本技术方案中，n指定为5，即为5折交叉验证；
[0129]
s42：将其中一份作为验证集，其余4份作为训练集，经过5次测试，每次都更换不同的验证集，得到5个模型结果，取最优结果；
[0130]
s43：使用最优超参数重新训练模型。
[0131]
s5：评分函数，包括以下过程：best_parameters模块：描述了已取得最佳结果的超参数的组合；best_score模块：提供优化过程期间观察到的最好的评分。
[0132]
在本技术方案中，基于nlmixr包构建的模型性能与超参数直接相关，这说明准确设定模型的超参数是非常必要的。该方法包括构建药动学
‑
药效学模型；确定超参数空间；通过搜索的方式来获得候选参数集合；交叉验证机制；评分函数。该装置包括基于nlmixr包的药动学
‑
药效学模型构建模块；参数空间生成模块；超参数自动优化模块；药动学
‑
药效学模型超参自动学习的性能评价模块。根据本技术方案提出的一种基于nlmixr包的药动学
‑
药效学模型超参数自动学习的方法，可以简单实现模型超参数自动调优问题，帮助构建更优的药动学
‑
药效学模型，以便加快药物研发进程和节约研发成本。
[0133]
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神和基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内，不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
[0134]
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。本发明尚有多种实施方式，凡采用等同变换或者等效变换而形成的所有技术方案，均落在本发明的保护范围之内。