Age-sex-specific Population forecasts are derived through stochastic population renewal using forecasts of mortality, fertility and net migration. Functional data models with time series coefficients are used to model age-specific mortality and fertility rates. As detailed migration data are lacking, net migration by age and sex is estimated as the diffierence between historic annual population data and successive populations one year ahead derived from a projection using fertility and mortality data. This estimate, which includes error. is also modeled using a functional data model. The three models involve different strengths of the general Box-Cox transformation chosen to minimise out-of-sample forecast error. Uncertainty is estimated from the model, with an adjustment to ensure that the one-step-forecast variances are equal to those obtained with historical data. The three models are then used in a Monte Carlo simulation of future fertility, mortality and net migration, which are combined using the cohort-component method to obtain age-specific forcecasts is of the population by sex. The distribution of the forecasts provides probabilistic prediction intervals. The method is demonstrated by making, 20-year forecasts using Australian data for the period 1921-2004. The advantages of our method are: (I) it is a coherent stochastic model of the three demographic components; (2) it is estimated entirely from historical data with no subjective inputs required: and (3) it provides probabilistic prediction intervals for any demographic variable that is derived from population numbers and vital events, including life expectancies, total fertility rates and dependency ratios. (C) 2008 International Institute of Forecasters. Published by Elsevier B.V. All rights reserved.