One major task in clinical pharmacology is to determine the pharmacokinetic-pharmacodynamic (PK-PD) parameters of a drug in a patient population. NONMEM is a program commonly used to build population PK-PD models, that is, models that characterize the relationship between a patient's PK-PD parameters and other patient specific covariates such as the patient's (patho)physiological condition, concomitant drug therapy, etc. This paper extends a previously described approach to efficiently find the relationships between the PK-PD parameters and covariates. In a first step, individual estimates of the PK-PD parameters are obtained as empirical Bayes estimates, based on a prior NONMEM fit using no covariates. In a second step, the individual PK-PD parameter estimates are regressed on the covariates using a generalized additive model. In a third and final step, NONMEM is used to optimize and finalize the population model. Four real-data examples are used to demonstrate the effectiveness of the approach. The examples show that the generalized additive model for the individual parameter estimates is a good initial guess for the NONMEM population model. In all four examples, the approach successfully selects the most important covariates and their functional representation. The great advantage of this approach is speed. The time required to derive a population model is markedly reduced because the number of necessary NONMEM runs is reduced. Furthermore, the approach provides a nice graphical representation of the relationships between the PK-PD parameters and covariates.
[1]
L B Sheiner,et al.
Estimating population kinetics.
,
1982,
Critical reviews in biomedical engineering.
[2]
H. Akaike.
A new look at the statistical model identification
,
1974
.
[3]
Jerome H. Friedman.
Multivariate adaptive regression splines (with discussion)
,
1991
.
[4]
Carl de Boor,et al.
A Practical Guide to Splines
,
1978,
Applied Mathematical Sciences.
[5]
D. Stanski,et al.
A three-step approach combining bayesian regression and NONMEM population analysis: Application to midazolam
,
1991,
Journal of Pharmacokinetics and Biopharmaceutics.
[6]
Leo Breiman,et al.
[The ∏ Method for Estimating Multivariate Functions from Noisy Data]: Response
,
1991
.
[7]
R. Bellman,et al.
V. Adaptive Control Processes
,
1964
.
[8]
L B Sheiner,et al.
Forecasting individual pharmacokinetics
,
1979,
Clinical pharmacology and therapeutics.
[9]
J. Freidman,et al.
Multivariate adaptive regression splines
,
1991
.
[10]
R. Tibshirani,et al.
Generalized Additive Models
,
1991
.
[11]
L B Sheiner,et al.
Bayesian individualization of pharmacokinetics: simple implementation and comparison with non-Bayesian methods.
,
1982,
Journal of pharmaceutical sciences.