We describe population and individual midazolam pharmacokinetic parameter estimates in pediatric intensive care patients by using a population modeling approach. Lower midazolam elimination was observed in comparison to other studies in pediatric intensive care patients, probably as a result of differences in study design and patient differences such as age and disease state. Covariates such as renal failure, hepatic failure, and concomitant administration of CYP3A inhibitors are important predictors of altered midazolam and metabolite pharmacokinetics in pediatric intensive care patients. The derived population model can be useful for future dose optimization and Bayesian individualization.