RPEM: Randomized Monte Carlo parametric expectation maximization algorithm
Open Access
- 15 April 2024
- journal article
- research article
- Published by Wiley in CPT: Pharmacometrics & Systems Pharmacology
- Vol. 13 (5) , 759-780
Abstract
Inspired from quantum Monte Carlo, by sampling discrete and continuous variables at the same time using the Metropolis–Hastings algorithm, we present a novel, fast, and accurate high performance Monte Carlo Parametric Expectation Maximization (MCPEM) algorithm. We named it Randomized Parametric Expectation Maximization (RPEM). We compared RPEM with NONMEM's Importance Sampling Method (IMP), Monolix's Stochastic Approximation Expectation Maximization (SAEM), and Certara's Quasi‐Random Parametric Expectation Maximization (QRPEM) for a realistic two‐compartment voriconazole model with ordinary differential equations using simulated data. We show that RPEM is as fast and as accurate as the algorithms IMP, QRPEM, and SAEM for the voriconazole model in reconstructing the population parameters, for the normal and log‐normal cases.Keywords
All Related Versions
This publication has 29 references indexed in Scilit:
- A Regression Approach to Visual Predictive Checks for Population Pharmacometric ModelsCPT: Pharmacometrics & Systems Pharmacology, 2018
- Model Evaluation of Continuous Data Pharmacometric Models: Metrics and GraphicsCPT: Pharmacometrics & Systems Pharmacology, 2016
- Comparing the performance of FOCE and different expectation-maximization methods in handling complex population physiologically-based pharmacokinetic modelsJournal of Pharmacokinetics and Pharmacodynamics, 2016
- Introduction to Population Pharmacokinetic / Pharmacodynamic Analysis with Nonlinear Mixed Effects Models2014
- An improved SAEM algorithm for maximum likelihood estimation in mixtures of non linear mixed effects modelsStatistics and Computing, 2013
- Monte Carlo Methods2008
- Nonlinear random effects mixture models: Maximum likelihood estimation via the EM algorithmComputational Statistics & Data Analysis, 2007
- Maximum likelihood estimation in nonlinear mixed effects modelsComputational Statistics & Data Analysis, 2005
- Monte Carlo calculations of triton and 4He nuclei with the Reid potentialNuclear Physics A, 1981
- Evaluation of methods for estimating population pharmacokinetic parameters. I. Michaelis-menten model: Routine clinical pharmacokinetic dataJournal of Pharmacokinetics and Biopharmaceutics, 1980