We investigate the personalisation and prediction accuracy of mathematical models for white blood cell (WBC) count dynamics during consolidation treatment using intermediate or high-dose cytarabine (Ara-C) in acute myeloid leukemia (AML). Ara-C is the clinically most relevant cytotoxic agent for AML treatment.We extend the gold-standard model of myelosuppression and a pharmacokinetic model of Ara-C with different hypotheses of Ara-C's pharmacodynamic effects. We cross-validate 12 mathematical models using dense WBC count measurements from 23 AML patients. Surprisingly, the prediction accuracies are similarly good despite different modelling hypotheses. Therefore, we compare average clinical and calculated WBC recovery times for different Ara-C schedules as a successful methodology for model discrimination. As a result, a new hypothesis of a secondary pharmacodynamic effect on the proliferation rate seems plausible. Furthermore, we demonstrate how personalized predictions of the impact of treatment timing on subsequent nadir values could be used for clinical decision support.Preamble on terminology and potentially confusing 1 synonyms 2 Our work is located in the intersection of mathematics, control theory, systems biology, 3 pharmacology, and medicine. Words like "model" or "parameter" have different 4 meanings in these scientific communities, and similar concepts have different names like 5