Objective: The present study aims to establish a population pharmacokinetic model of ganciclovir and optimize the dosing regimen in critically ill children suffering from cytomegalovirus related disease.Methods: A total of 104 children were included in the study. The population pharmacokinetic model was developed using the Phoenix NLME program. The final model was validated by diagnostic plots, nonparametric bootstrap, visual predictive check, and normalized prediction distribution errors. To further evaluate and optimize the dosing regimens, Monte Carlo simulations were performed. Moreover, the possible association between systemic exposure and hematological toxicity were also monitored in the assessment of adverse events.Results: The ganciclovir pharmacokinetics could be adequately described by a one-compartment model with first-order elimination along with body weight and estimated glomerular filtration rate as significant covariates. As showed in this study, the typical population parameter estimates of apparent volume of distribution and apparent clearance were 11.35 L and 5.23 L/h, respectively. Simulations indicated that the current regimen at a dosage of 10 mg/kg/d would result in subtherapeutic exposure, and elevated doses might be required to reach the target ganciclovir level. No significant association between neutropenia, the most frequent toxicity reported in our study (19.23%), and ganciclovir exposure was observed.Conclusion: A population pharmacokinetic model of intravenous ganciclovir for critically ill children with cytomegalovirus infection was successfully developed. Results showed that underdosing of ganciclovir was relatively common in critically ill pediatric patients, and model-based approaches should be applied in the optimizing of empiric dosing regimens.