Two general methods for population pharmacokinetic modeling: non-parametric adaptive grid and non-parametric Bayesian

Population pharmacokinetic (PK) modeling methods can be statistically classified as either parametric or nonparametric (NP). Each classification can be divided into maximum likelihood (ML) or Bayesian (B) approaches. In this paper we discuss the nonparametric case using both maximum likelihood and Bayesian approaches. We present two nonparametric methods for estimating the unknown joint population distribution of model parameter values in a pharmacokinetic/pharmacodynamic (PK/PD) dataset. The first method is the NP Adaptive Grid (NPAG). The second is the NP Bayesian (NPB) algorithm with a stick-breaking process to construct a Dirichlet prior. Our objective is to compare the performance of these two methods using a simulated PK/PD dataset. Our results showed excellent performance of NPAG and NPB in a realistically simulated PK study. This simulation allowed us to have benchmarks in the form of the true population parameters to compare with the estimates produced by the two methods, while incorporating challenges like unbalanced sample times and sample numbers as well as the ability to include the covariate of patient weight. We conclude that both NPML and NPB can be used in realistic PK/PD population analysis problems. The advantages of one versus the other are discussed in the paper. NPAG and NPB are implemented in R and freely available for download within the Pmetrics package from www.lapk.org.

[1]  B. Lindsay.  The Geometry of Mixture Likelihoods: A General Theory , 1983 .

[2]  A. Mallet.  A maximum likelihood estimation method for random coefficient regression models , 1986 .

[3]  D J Spiegelhalter,et al.  Flexible random‐effects models using Bayesian semi‐parametric models: applications to institutional comparisons , 2007, Statistics in medicine.

[4]  Bennett L. Fox,et al.  Algorithm 647: Implementation and Relative Efficiency of Quasirandom Sequence Generators , 1986, TOMS.

[5]  Stephen G. Walker,et al.  Sampling the Dirichlet Mixture Model with Slices , 2006, Commun. Stat. Simul. Comput..

[6]  David B. Dunson,et al.  Bayesian Nonparametrics: Nonparametric Bayes applications to biostatistics , 2010 .

[7]  W. Karush.  Minima of Functions of Several Variables with Inequalities as Side Conditions , 2014 .

[8]  L B Sheiner,et al.  Estimating population kinetics. , 1982, Critical reviews in biomedical engineering.

[9]  Alan Schumitzky,et al.  Accurate Detection of Outliers and Subpopulations With Pmetrics, a Nonparametric and Parametric Pharmacometric Modeling and Simulation Package for R , 2012, Therapeutic drug monitoring.

[10]  G. Roberts,et al.  Retrospective Markov chain Monte Carlo methods for Dirichlet process hierarchical models , 2007, 0710.4228.

[11]  K. Schittkowski,et al.  NONLINEAR PROGRAMMING , 2022 .

[12]  Hoon Kim,et al.  Monte Carlo Statistical Methods , 2000, Technometrics.

[13]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[14]  J. Kiefer,et al.  CONSISTENCY OF THE MAXIMUM LIKELIHOOD ESTIMATOR IN THE PRESENCE OF INFINITELY MANY INCIDENTAL PARAMETERS , 1956 .

[15]  Stephen G. Walker,et al.  Slice sampling mixture models , 2011, Stat. Comput..

[16]  Alan Schumitzky,et al.  Kullback-Leibler Markov Chain Monte Carlo - a New Algorithm for Finite Mixture Analysis and its Application to Gene Expression Data , 2008, J. Bioinform. Comput. Biol..

[17]  Bayesian Nonparametric Population Models: Formulation and Comparison with Likelihood Approaches , 2004, Journal of Pharmacokinetics and Biopharmaceutics.

[18]  Peter Müller,et al.  A Bayesian Population Model with Hierarchical Mixture Priors Applied to Blood Count Data , 1997 .

[19]  Xiaoning Wang,et al.  Population pharmacokinetic/pharmacodynamic mixture models via maximum a posteriori estimation , 2009, Comput. Stat. Data Anal..

[20]  Mats O. Karlsson,et al.  Two bootstrapping routines for obtaining imprecision estimates for nonparametric parameter distributions in nonlinear mixed effects models , 2011, Journal of Pharmacokinetics and Pharmacodynamics.

[21]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[22]  Jing Wang,et al.  Dirichlet Processes in Nonlinear Mixed Effects Models , 2010, Commun. Stat. Simul. Comput..

[23]  Roger W. Jelliffe,et al.  The Nonparametric Adaptive Grid Algorithm for Population Pharmacokinetic Modeling , 2013 .

[24]  Christian P. Robert,et al.  The Bayesian choice , 1994 .

[25]  Roger W. Jelliffe,et al.  An adaptive grid non-parametric approach to pharmacokinetic and dynamic (PK/PD) population models , 2001, Proceedings 14th IEEE Symposium on Computer-Based Medical Systems. CBMS 2001.

[26]  L Aarons,et al.  Population pharmacokinetics. , 1992, International journal of clinical pharmacology, therapy, and toxicology.

[27]  Sylvia Frühwirth-Schnatter,et al.  Finite Mixture and Markov Switching Models , 2006 .

[28]  A. Schumitzky.  Nonparametric EM Algorithms for estimating prior distributions , 1991 .

[29]  Peter Müller,et al.  Bayesian Population Pharmacokinetic and Pharmacodynamic Analyses Using Mixture Models , 1997, Journal of Pharmacokinetics and Biopharmaceutics.

[30]  Marc Lavielle,et al.  Estimation of Population Pharmacokinetic Parameters of Saquinavir in HIV Patients with the MONOLIX Software , 2007, Journal of Pharmacokinetics and Pharmacodynamics.

[31]  Xiaoning Wang,et al.  Nonlinear random effects mixture models: Maximum likelihood estimation via the EM algorithm , 2007, Comput. Stat. Data Anal..

[32]  Lancelot F. James,et al.  Gibbs Sampling Methods for Stick-Breaking Priors , 2001 .

[33]  D Bayard,et al.  Achieving target goals most precisely using nonparametric compartmental models and "multiple model" design of dosage regimens. , 2000, Therapeutic drug monitoring.

[34]  J. Sethuraman.  A CONSTRUCTIVE DEFINITION OF DIRICHLET PRIORS , 1991 .

[35]  Pulak Ghosh,et al.  A semi‐parametric Bayesian approach to average bioequivalence , 2007, Statistics in medicine.

[36]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.