The present study introduces new regression formulae that address several challenges of current subadult stature estimation methods by 1) using a large, contemporary, cross-sectional sample of subadult skeletal remains; 2) generating regression models using both lengths and breadths; 3) utilizing both linear and nonlinear regression models to accommodate the nonlinear shape of long bone growth; and 4) providing usable prediction intervals for estimating stature. Eighteen long bone measurements, stature, and age were collected from computed tomography images for a sample of individuals (n = 990) between birth and 20 years from the United States. The bivariate relationship between long bone measurements and stature was modeled using linear and nonlinear methods on an 80% training sample and evaluated on a 20% testing sample. Equations were generated using pooled-sex samples. Goodness of fit was evaluated using Kolmogorov-Smirnov tests and mean absolute deviation (MAD). Accuracy and precision were quantified using percent testing accuracy and Bland-Altman plots. In total, 38 stature estimation equations were created and evaluated, all achieving testing accuracies greater than 90%. Nonlinear models generated better fits compared to linear counterparts and generally produced smaller MAD (3.65 - 15.90cm). Length models generally performed better than breadth models, and a mixture of linear and nonlinear methods resulted in highest testing accuracies. Model performance was not biased by sex, age, or measurement type. A freely available, online graphical user interface is provided for immediate use of the models by practitioners in forensic anthropology and will be expanded to include bioarchaeological contexts in the future.