Mixed Effects Modeling Using Stochastic Differential Equations: Illustrated by Pharmacokinetic Data of Nicotinic Acid in Obese Zucker Rats
Journal Title: The AAPS Journal - Year 2015, Vol 17, Issue 3
Abstract
Inclusion of stochastic differential equations in mixed effects models provides means to quantify and distinguish three sources of variability in data. In addition to the two commonly encountered sources, measurement error and interindividual variability, we also consider uncertainty in the dynamical model itself. To this end, we extend the ordinary differential equation setting used in nonlinear mixed effects models to include stochastic differential equations. The approximate population likelihood is derived using the first-order conditional estimation with interaction method and extended Kalman filtering. To illustrate the application of the stochastic differential mixed effects model, two pharmacokinetic models are considered. First, we use a stochastic one-compartmental model with first-order input and nonlinear elimination to generate synthetic data in a simulated study. We show that by using the proposed method, the three sources of variability can be successfully separated. If the stochastic part is neglected, the parameter estimates become biased, and the measurement error variance is significantly overestimated. Second, we consider an extension to a stochastic pharmacokinetic model in a preclinical study of nicotinic acid kinetics in obese Zucker rats. The parameter estimates are compared between a deterministic and a stochastic NiAc disposition model, respectively. Discrepancies between model predictions and observations, previously described as measurement noise only, are now separated into a comparatively lower level of measurement noise and a significant uncertainty in model dynamics. These examples demonstrate that stochastic differential mixed effects models are useful tools for identifying incomplete or inaccurate model dynamics and for reducing potential bias in parameter estimates due to such model deficiencies.
Authors and Affiliations
Jacob Leander, Joachim Almquist, Christine Ahlström, Johan Gabrielsson, Mats Jirstrand
Keywords
Related Articles
- EP ID EP680945
- DOI 10.1208/s12248-015-9718-8
- Views 44
- Downloads 0
An Electrically Tight In Vitro Blood–Brain Barrier Model Displays Net Brain-to-Blood Efflux of Substrates for the ABC Transporters, P-gp, Bcrp and Mrp-1
Efflux transporters of the ATP-binding cassette superfamily including breast cancer resistance protein (Bcrp/Abcg2), P-glycoprotein (P-gp/Abcb1) and multidrug resistance-associated proteins (Mrp’s/Abcc’s)...
Development of an In Vivo Ovine Dry Powder Inhalation Model for the Evaluation of Conventional and Controlled Release Microparticles
Microparticles of DSCG alone or with 90% w/w PVA where prepared and characterised as described previously (5,6). Briefly, DSCG was spray dried under optimised conditions to produce a yield of microparticles with desired...
Identification of Key Licorice Constituents Which Interact with Cytochrome P450: Evaluation by LC/MS/MS Cocktail Assay and Metabolic Profiling
The online version of this article (doi:10.1208/s12248-013-9544-9) contains supplementary material, which is available to authorized users.
Conformational Analysis of Therapeutic Proteins by Hydroxyl Radical Protein Footprinting
The online version of this article (doi:10.1208/s12248-012-9336-7) contains supplementary material, which is available to authorized users.
Simultaneous population optimal design for pharmacokinetic-pharmacodynamic experiments
Multiple outputs or measurement types are commonly gathered in biological experiments. Often, these experiments are expensive (such as clinical drug trials) or require careful design to achieve the desired information co...