IV-29 Sungwoo Goo

Explainable covariate screening method using Artificial Neural Network (ANN) for Non-linear Mixed Effect Modeling

Sungwoo Goo(1), Heeyoon Jung(1), Woojin Jung(2), Jung-woo Chae*(1,2), Hwi-yeol Yun*(1,2)

(1)Department of Bio-AI convergence, Chungnam National University, Daejeon, Republic of Korea (2)College of Pharmacy, Chungnam National University, Daejeon, Republic of Korea

Objectives: 

As the fast growing of deep learning methodologies, there has been various attempt to incorporate deep learning method into pharmacometrics. Nevertheless, it is still difficult to use the deep learning library into the NLME model library at the same time in the Python circumstance, because of language differences. To overcome those limitation, we developed a non-linear mixed model library based on pytorch and efficient covariate screening method using an artificial neural network(ANN).

Methods:

The NLME model could be implemented with ANN, since linear algebra operated in both of them. The Objective function valeues (OFV) was not only calculated by referring to the paper of Bae at al.(2016)[1], but the prediction function was also conducted using the symbolic math library, ‘sympy’. In addition, ‘torchdiffeq’ was used as a module for solving Ordinary differential equation (ODE) what it was based on Runge–Kutta methods and Dormand–Prince methods. and Limited-memory- Broyden–Fletcher–Goldfarb–Shanno algorithm was conducted as well for optimization algorithm. Theophylline dataset and PK base model, reported by Boeckmann at al.(1994)[2], was evaluated for verification of the develped library.

Then, covariate screening method using ANN was attempted to estimate the relationship among covariates and PK parameters (Ka; absorption rate constant, Kel; elimination rate constant, and Vd: volume of distribution). The real covariate (Body weight, BWT) and dummy covariates having with or without relationship on BWT were tested as the covariates. Covariate screening was performed 100 times with randomly generated 100 datasets and counted the number of covariate selection.

Results: 

It was confirmed that the NLME model implemented in the Pytoch environment showed similar performance compared to NONMEM. (Ex., difference of OFV and PK parameters < 3%) The results of the VPC using the Pytoch NLME model and NONMEM looks similar as well. For the covariate screening method using ANN, the BWT was not only selected, similar to the results of the original model using the existing NONMEM, but it also was confirmed that dummy covariates having a correlation with BWT were also selected depending on the correlation of BWT. And covariate searching was attempted as a backward elimination method using ANN.
(Counts for selected covariates in 100 times covariate searching. ‘covariate name’: counts for selection / ‘BWT_Norm’: 82 /  ‘BWT_Cor_0.5(dummy covariates having correlation 0.5 with BWT)’: 71 / ‘BWT_Cor_0.1(dummy covariates having correlation 0.1 with BWT)’: 63 / ‘rand-1+1(random number between -1 and ~ 1)’: 62 / ‘norm(0,1)( random variables of standard normal distribution’: 55 / ‘fixed(zero numbers)’ : 0)

NLME model optimization was performed using torchpm. And for each ID, simulation was performed 300 times to draw a graph. The graph was plotted by displaying the 5 percentile, median, mean, and 95 percentile.

Conclusions:

NLME model was successfully developed by implementation in the PyTorch environment and efficient covariate screening method using ANN in PyTorch was well conducted. Considering the cons of traditional covariate searching methods like stepwise covariate model(SCM), this methods could be used for alternative method for replace traditional covariate screening method for its strength on less time consuming and high performance even correlated covariates. In addition, it might contribute to the integration of deep learning into pharmacometrics. The entire code is shared on github. (https://github.com/yeoun9/torchpm) Those who are interested, please contribute to the code.

References:
[1] Bae, K.-S. and D.-S. Yim (2016). “R-based reproduction of the estimation process hidden behind
NONMEM® Part 2: First-order conditional estimation”. In: Translational and Clinical Pharma-
cology 24.4, pp. 161–168.
[2] Boeckmann, A. J., L. B. Sheiner, and S. L. Beal (1994). NONMEM users guide: Part V.

Reference: PAGE 30 (2022) Abstr 10000 [www.page-meeting.org/?abstract=10000]

Poster: Methodology – AI/Machine Learning

PDF poster / presentation (click to open)