Optimal sampling schedule design for populations of patients

Generation of pharmacodynamic relationships in the clinical arena requires estimation of pharmacokinetic parameter values for individual patients. When the target population is severely ill, the ability to obtain traditional intensive blood sampling schedules is curtailed. Population modeling guided by optimal sampling theory has provided robust estimates of individual patient pharmacokinetic parameter values. Because of the wide range of parameter values seen in this circumstance, it is important to know how the range of parameter values in the population affects the timing of the optimal samples. We describe a new, simple technique to obtain optimal samples for a population of patients. This technique uses the nonparametric distribution associated with a nonparametric adaptive grid population pharmacokinetic analysis. We used the distribution from an analysis of 58 patients receiving levofloxacin for nosocomial pneumonia at a dose of 750 mg. The collection of parameter vectors and their associated probabilities were entered into a D-optimal design evaluation by using ADAPT II. The sampling times, weighted for their probabilities, were displayed in a frequency histogram (an expression of how system information varies with time for the population). Such an explicit expression of the time distribution of information allows rational sampling design that is robust not only for the population mean vector, as in traditional D-optimal design theory, but also for large portions of the total population. For levofloxacin, one reasonable six-sample design would be 1.5, 2, 2.25, 4, 4.75, and 24 h after starting a 90-min infusion. Such sampling designs allow informative population pharmacokinetic analysis with precise and unbiased estimates after the maximal a posteriori probability Bayesian step. This allows the highest probability of delineating a pharmacodynamic relationship.