Simultaneous versus sequential pharmacokinetic‐pharmacodynamic population analysis using an Iterative Two‐Stage Bayesian technique

Abstract

A method for simultaneous pharmacokinetic-pharmacodynamic (PK-PD) population analysis using an Iterative Two-Stage Bayesian (ITSB) algorithm was developed. The method was evaluated using clinical data and Monte Carlo simulations. Data from a clinical study with rocuronium in nine anesthetized patients and data generated by Monte Carlo simulation using a similar study design were analysed by sequential PK-PD analysis, PD analysis with nonparametric PK data and simultaneous PK-PD analysis. Both PK and PD data sets were 'rich' with respect to the number of measurements per individual. The accuracy and precision of the estimated population parameters were evaluated by comparing their mean error (ME) and root mean squared error (RMSE), respectively. The influence of PD model misspecification on the results was also investigated. The simultaneous PK-PD analysis resulted in slightly more precise population parameter estimates than the sequential PK-PD analysis and the nonparametric PK method. In the presence of PD model misspecification, however, simultaneous analysis resulted in poor PK parameter estimates, while sequential PK-PD analysis performed well. In conclusion, ITSB is a valuable technique for PK-PD population analysis of rich data sets. The sequential PK-PD method is better suited for the analysis of rich data than the simultaneous analysis.

A pharmacokinetic–pharmacodynamic model for intrathecal baclofen in patients with severe spasticity

“…Our PKPD model on ITB is the first one based on human data. The observed concentration data were used instead of the PK model‐predicted data, because of the high inter‐individual and intra‐individual differences in the model parameters, as well as to prevent PK modelling errors to be carried over to the PKPD model . The PKPD model is based on a delayed effect model and shows an adequate fit to the effect data as shown in Figure .…”

Multiple intravenous doses of paracetamol result in a predictable pharmacokinetic profile in very preterm infants

“…The software package takes into account individual patient characteristics (gestational age, birth weight, sex en creatinine) for prediction of the time–concentration profile. The KinPop module of the program uses an iterative two‐stage Bayesian fitting and calculates means, medians and standard deviation (SD) of the pharmacokinetic parameters . Before the actual modelling starts, rough estimates of the model parameters and their SD are entered from basic pharmacokinetics of paracetamol in adults.…”

Pharmacokinetic/pharmacodynamic data extrapolation models for improved pediatric efficacy and toxicity estimation, with application to secondary hyperparathyroidism

“…In modeling PK‐PD data, a choice must be made to analyze the PK and PD data simultaneously or sequentially 32‐36 . Simultaneous modeling may be considered optimal from a hierarchical Bayesian point of view, but PD models are often highly complex (involving systems of differential equations), so even computing joint maximum likelihood estimates can be numerically difficult.…”