The objective of this study was to determine whether any pretreatment parameters were associated with pharmacological effect or toxicity parameters after vincristine administration and to describe a mathematical model, which explains the interpatient pharmacodynamic variability. The relationship between patient characteristics and vincristine dose and hematological toxicity were evaluated. 68 pediatric and adolescence patients and 107 adults with acute lymphoblastic leukemia were treated with vincristine $1.5mg/m^2/day$ IV and other anticancer drugs as scheduled. Complete blood counts and other blood test results were obtained. The input variables were age, gender, weight, lean body weight (LBW), height, body surface area, vincristine dose and total vincristine dose. The outcome measures were nadir values (white blood cells, absolute neutrophil counts, hemoglobin, and platelets); the absolute decrease, relative decrease, and survival fraction of blood cells. Polynomial regression analysis was carried out to determine the other significant covariates. The variability of $WBC_{nadir}$ was modeled with good precision and accuracy with a two-covariate model. This model should be validated and improved on with further clinical data. We believe that such pharmacodynamic modeling should be explored further to determine its performance and clinical relevance compared with modeling using pharmacokinetic parameter.