Introduction: Nonparametric population modeling algorithms have a theoretical superiority over parametric methods to detect pharmacokinetic and pharmacodynamic subgroups and outliers within a study population. Methods: The authors created “Pmetrics,” a new Windows and Unix R software package that updates the older MM-USCPACK software for nonparametric and parametric population modeling and simulation of pharmacokinetic and pharmacodynamic systems. The parametric iterative 2-stage Bayesian and the nonparametric adaptive grid (NPAG) approaches in Pmetrics were used to fit a simulated population with bimodal elimination (K_el) and unimodal volume of distribution (V_d), plus an extreme outlier, for a 1-compartment model of an intravenous drug. Results: The true means (SD) for K_el and V_d in the population sample were 0.19 (0.17) and 102 (22.3), respectively. Those found by NPAG were 0.19 (0.16) and 104 (22.6). The iterative 2-stage Bayesian estimated them to be 0.18 (0.16) and 104 (24.4). However, given the bimodality of K_el, no subject had a value near the mean for the population. Only NPAG was able to accurately detect the bimodal distribution for K_el and to find the outlier in both the population model and in the Bayesian posterior parameter estimates. Conclusions: Built on over 3 decades of work, Pmetrics adopts a robust, reliable, and mature nonparametric approach to population modeling, which was better than the parametric method at discovering true pharmacokinetic subgroups and an outlier.

We hypothesized that dosing vancomycin to achieve trough concentrations of >15 mg/liter overdoses many adults compared to area under the concentration-time curve (AUC)-guided dosing. We conducted a 3-year, prospective study of vancomycin dosing, plasma concentrations, and outcomes. In year 1, nonstudy clinicians targeted trough concentrations of 10 to 20 mg/liter (infection dependent) and controlled dosing. In years 2 and 3, the study team controlled vancomycin dosing with BestDose Bayesian software to achieve a daily, steady-state AUC/MIC ratio of ≥400, with a maximum AUC value of 800 mg · h/liter, regardless of trough concentration. For Bayesian estimation of AUCs, we used trough samples in years 1 and 2 and optimally timed samples in year 3. We enrolled 252 adults who were ≥18 years old with ≥1 available vancomycin concentration. Only 19% of all trough concentrations were therapeutic versus 70% of AUCs ( P < 0.0001). After enrollment, median trough concentrations by year were 14.4, 9.7, and 10.9 mg/liter ( P = 0.005), with 36%, 7%, and 6% over 15 mg/liter ( P < 0.0001). Bayesian AUC-guided dosing in years 2 and 3 was associated with fewer additional blood samples per subject (3.6, 2.0, and 2.4; P = 0.003), shorter therapy durations (8.2, 5.4, and 4.7 days; P = 0.03), and reduced nephrotoxicity (8%, 0%, and 2%; P = 0.01). The median inpatient stay was 20 days among nephrotoxic patients versus 6 days ( P = 0.002). There was no difference in efficacy by year, with 42% of patients having microbiologically proven infections. Compared to trough concentration targets, AUC-guided, Bayesian estimation-assisted vancomycin dosing was associated with decreased nephrotoxicity, reduced per-patient blood sampling, and shorter length of therapy, without compromising efficacy. These benefits have the potential for substantial cost savings. (This study has been registered at ClinicalTrials.gov under registration no. NCT01932034.)