PKPDindex: an R/CRAN package for PK/PD index analysis

Abstract

PKPDindex() is a free and open-source R package for analysing pharmacokinetic/pharmacodynamic (PK/PD) indices, including AUC/MIC, Cmax/MIC, and T>MIC. Development was motivated by previously identified inconsistencies in PK/PD index modelling and reporting practices. The package fits eight variations of the Emax model to data and compares model performance using the Akaike information criterion (AIC) andR² values. As input, the package requires a dataset containing PK/PD indices and a response variable. As output, it generates tables summarising all fitted models, identifies the best-fitting model for each index, and produces customisable plots of model fits with parameter estimates. PKPDindex() is freely available on CRAN (https://CRAN.R-project.org/package=PKPDindex).

Keywords

pharmacokinetic/pharmacodynamic (PK/PD) index / Emax model / R package / antimicrobial pharmacology / model fitting / dose fractionation studies / PKPDindex

References

[1] Mouton JW, Dudley MN, Cars O, Derendorf H, Drusano GL. Standardization of pharmacokinetic/pharmacodynamic (PK/PD) terminology for anti—infective drugs: an update. J Antimicrob Chemother. 2005; 55: 601-7.

[2] Nielsen EI, Cars O, Friberg LE. Pharmacokinetic/Pharmacodynamic (PK/PD) Indices of Antibiotics Predicted by a Semimechanistic PKPD Model: a Step toward Model—Based Dose Optimization. Antimicrob Agents Chemother. 2011; 55: 4619—30.

[3] Nielsen EI, Friberg LE. Pharmacokinetic—Pharmacodynamic Modeling of Antibacterial Drugs. Pharmacol Rev. 2013; 65: 1053-90.

[4] Rotschafer JC, Andes DR, Rodvold KA. Antibiotic Pharmacodynamics. New York: Humana Press; 2016.

[5] Goutelle S, Maurin M, Rougier F, Barbaut X, Bourguignon L, Ducher M, et al. The Hill equation: a review of its capabilities in pharmacological modelling. Fundam Clin Pharmacol. 2008; 22: 633-48.

[6] Alabdulkarim N, Readman J, Mead A, Oakley R, Wenker S, Kembou J, et al. Systematic review of antimicrobial pharmacokinetic/pharmacodynamic indices in murine thigh and hollow fibre dose fractionation studies analysed with a standard method. J Antimicrob Chemother. 2025; 81: e81.

[7] Spiess AN, Neumeyer N. An evaluation of R2 as an inadequate measure for nonlinear models in pharmacological and biochemical research: a Monte Carlo approach. BMC Pharmacol. 2010; 10: 6.

[8] Burnham KP, Anderson DR. Multimodel Inference. Sociol Methods Res. 2004; 33: 261-304.

[9] Akaike H. Information theory and an extension of the maximum likelihood principle. In: Petrov BN, Csaki F, editors. Budapest (Hungary): Akadémiai Kiadó; 1973. pp. 267-81.

[10] Kuha J. AIC and BIC. Sociol Methods Res. 2004; 33: 188-229.

[11] Ting N. Dose finding in drug development. Springer Science & Business Media; 2006.

[12] Schwarz G. Estimating the Dimension of a Model. Ann Stat. 1978; 6: e6.