This study aimed to suggest an initial pediatric vancomycin dose regimen through population pharmacokinetic-pharmacodynamic modeling. A population pharmacokinetic approach was used to analyze vancomycin concentration-time data from a large pediatric cohort. Pharmacokinetic target attainment for patients with bloodstream isolates was compared with clinical outcome using logistic regression and classification and regression trees. Change in serum creatinine during treatment was used as an indicator of acute nephrotoxicity. Probability of acute kidney injury (50% increase from baseline) or kidney failure (75% increase from baseline) was evaluated using logistic regression. An initial dosing regimen was derived, personalized by age, weight, and serum creatinine, using stochastic simulations. Data from 785 hospitalized pediatric patients (1 day to 21 years of age) with suspected Gram-positive infections were collected. Estimated (relative standard error) typical clearance, volume of distribution 1, intercompartmental clearance, and volume of distribution 2 were (standardized to 70 kg) 4.84 (2.38) liters/h, 39.9 (8.15) liters, 3.85 (17.3) liters/h, and 37.8 (10.2) liters, respectively. While cumulative vancomycin exposure correlated positively with the development of nephrotoxicity (713 patients), no clear relationship between vancomycin area under the plasma concentration-time curve and efficacy was found (102 patients). Predicted probability of acute kidney injury and kidney failure with the optimized dosing regimen at day 5 was 10 to 15% and 5 to 10%, increasing by approximately 50% on day 7 and roughly 100% on day 10 across all age groups. This study presents the first data-driven pediatric dose selection to date accounting for nephrotoxicity, and it indicates that cumulative vancomycin exposure best describes risk of acute kidney injury and acute kidney failure.