Comparing non-hierarchical models: Application to non-linear mixed effects modeling

There is no method available to compare the fit of two non-hierarchical non-linear mixed effects models, although the common practice is to select the model with the lower objective function. Bootstrapping the log-likelihood differences (LLDs) of non-hierarchical models and constructing a bootstrap confidence interval on the LLDs is proposed for comparing the goodness-of-fit of such models. This is illustrated with different parameterizations of clearance models for an anti-infective agent in a longitudinal pharmacokinetic study which are compared. Additive and exponential models of creatinine clearance as a predictor of clearance are used as examples.