An adaptive grid non-parametric approach to pharmacokinetic and dynamic (PK/PD) population models

摘要

Our NPEM software for non-parametric PK/PD population modeling employs the classical EM Optimization algorithm to compute a maximum likelihood distribution on a large multidimensional grid. In order to achieve good resolution, a large number of grid points must be chosen, which can lead to high computational demands requiring a large-scale parallel supercomputer. Here we describe an improved method NPAG that uses a sequence of adaptively refitted grids, as well as a new, state-of-the-art interior point algorithm for solving the associated maximum likelihood problem on each successive grid. The combination of the adaptive grid strategy with the interior point algorithm is far foster than the original NPEM method. Also, NPAG requires much less memory, thus making many computations feasible on a PC or workstation that previously required supercomputer resources. Finally, the new algorithm easily and naturally accommodates the simultaneous maximum likelihood estimation of both intra- individual and inter-individual variability, thus improving usability and removing a major limitation of the original NPEM pro grain.