demography (version 1.22)

pop.sim: Population simulation

Description

Simulate future sample paths of a population using functional models for mortality, fertility and migration.

Usage

pop.sim(mort, fert = NULL, mig = NULL, firstyearpop, N = 100,
  mfratio = 1.05, bootstrap = FALSE)

Value

A list of two arrays containing male and female future simulated population values. The arrays are of dimension (p,h,N) where p is the number of age groups, h is the forecast horizon and N is the number of simulated sample paths.

Arguments

mort

Forecasts of class fmforecast2 for mortality.

fert

Forecasts of class fmforecast for female fertility.

mig

Forecasts of class fmforecast2 for net migration.

firstyearpop

Population for first year of simulation.

N

Number of sample paths to simulate.

mfratio

Male-female ratio used in distributing births.

bootstrap

If TRUE, simulation uses resampled errors rather than normally distributed errors.

Author

Rob J Hyndman

See Also

simulate.fmforecast, simulate.fmforecast2.

Examples

Run this code
if (FALSE) {
require(addb)
# Construct data objects
mort.sm <- smooth.demogdata(set.upperage(extract.years(australia,1950:2002),100))
fert.sm <- smooth.demogdata(extract.years(aus.fertility,1950:2002))
aus.mig <- netmigration(set.upperage(australia,100),aus.fertility,mfratio=1.0545)
# Fit models
mort.fit <- coherentfdm(mort.sm)
fert.fit <- fdm(fert.sm)
mig.fit <- coherentfdm(aus.mig)
# Produce forecasts
mort.fcast <- forecast(mort.fit)
fert.fcast <- forecast(fert.fit)
mig.fcast <- forecast(mig.fit)
# Simulate
aus.sim <- pop.sim(mort.fcast,fert.fcast,mig.fcast,australia)}

Run the code above in your browser using DataCamp Workspace