1. Demography Dynamic microsimulation
Download icon

Stay on Target: Population Projections and Microsimulation Design

  1. Rachel J Bacon  Is a corresponding author
  2. George Hodulik
  3. David Voas
  4. Ivan Puga-Gonzalez
  5. Wesley J Wildman
  1. Center for Mind and Culture, United States
  2. University College London, United Kingdom
  3. NORCE Center for Modelling Social Systems, Norway
Research article
Cite this article as: R. J Bacon, G. Hodulik, D. Voas, I. Puga-Gonzalez, W. J Wildman; 2023; Stay on Target: Population Projections and Microsimulation Design; International Journal of Microsimulation; 16(3); 77-99. doi: 10.34196/ijm.00289
7 figures and 5 tables

Figures

Results of Norway Models, scaled up to actual population size
Results of USA Models, scaled up to actual population size
Results of India Models, scaled up to actual population size
Population Dynamics Among Females of Child-bearing Ages (15-49), UN Data WPP2019
Results of Norway Models with Additional Stochasticity Options, scaled up to actual population size
Results of USA Models with Additional Stochasticity Options, scaled up to actual population size
Results of India Models with Additional Stochasticity Options, scaled up to actual population size

Tables

Table 1
Assignment of Demographic Event Risk by Age for 1-year Split Fertility Design
t=0Rnd. 2t=1Rnd. 2t=2Rnd. 2t=3Rnd. 2t=4Rnd. 2
Agent Age Characteristics
Current Age43444445454646474748
Starting Age4343434343
Ending Age4848484848
Corresponding Age Labels for UN Statistics
Fertility Rate40-4440-4440-4445-4945-4945-4945-4945-4945-4945-49
Survival Ratio40-4440-4440-4440-4440-44
Migration Count45-4945-4945-4945-4945-49
Table 2
Years of Risk among Newborn Cohorts
CohortYear 0Year 1Year 2Year 3Year 4Total Years of RiskPercent of Total
A11111533
B1111427
C111320
D11213
E117
Total15100
Table 3
Model Design Characteristics and Formulas
Intuitive ModelsSplit Fertility Models
Time5-year step and interval
t =T= [1950,1955,…2100]
1-year step t= [0,1,2,3,4],
5-year interval T=[1950,1955,…2100]
5-year step and interval
t =T= [1950,1955,…2100]
1-year step t= (0,1,2,3,4),
5-year interval T=[1950,1955,…2100]
Event ordering0) Initialization,
1) fertility,
2) infant initialization,
3) mortality,
4) aging,
5) migration,
6) time step/interval advancement, T+5.
0) Initialization,
1) fertility,
2) infant initialization,
3) mortality,
4) aging,
5) migration,
6) time step advancement,
7) repeat events 1-6 for t<5 until interval advancement, T+5.
0) Initialization,
1) fertility (round one),
2) infant initialization (round one),
3) mortality,
4) aging,
5) migration
6) fertility round two,
7) infant initialization (round two),
8) second round infant mortality,
9) second round infant aging
10) time step/interval advancement, T+5.
0) Initialization,
1) fertility (round one),
2) infant initialization (round one),
3) mortality,
4) aging,
5) migration
6) fertility round two,
7) infant initialization (round two),
8) second round infant mortality,
9) second round infant aging
10) time step advancement,
11) repeat events 1-10 for t<5 until interval advancement, T+5.
Current age (A)5-year age groups, A=[0,5…105]1-year age groups A=[0,1,…105]5-year age groups, A=[0,5…105]1-year age groups A=[0,1,…105]
For immigrants: A=eA-t
Initial age (iA)Not ApplicableNot ApplicableNot ApplicableiA=A-t
End age (eA)Not ApplicableNot ApplicableNot ApplicableeA=A+(5-t)
Model Age-Specific Fertility Rates (ASFR)(UNASFR)10005(UNASFR)1000;
Apply to current age (A)
(UNASFR)100052(UNASFR)100012; Apply to current age (A)
Infant initializationAssigned age -5, and sex based on UN sex ratioAssigned age -1, and sex based on UN sex ratioAssigned age -5, and sex based on UN sex ratioAssigned age -1, and sex based on UN sex ratio.
Model Mortality Rate (MR)1 - (UNSR)1(UNSR)5;
Apply to current age (A)
1 - (UNSR)1(UNSR)5; Apply to initial age (iA).
For newborn cohorts: 1UNSR[0]3 ;
For agents age 110+: 1.
For immigrants: no mortality
Model Emigration count (EC)(UNEC) * (Sample%)UNEC*Sample%25;
Apply to current age (A)
(UNEC) * (Sample%)(UNEC)(Sample\% )25; Apply to end age (eA).
For newborn cohorts:(UNEC)(Sample\% )15
Model Immigration count (IC)(UNIC) * (Sample%)UNIC*Sample%25;
Apply to current age (A)
(UNIC) * (Sample%)(UNIC)(Sample\% )25; Apply to end age (eA).
For newborn cohorts:(UNIC)(Sample\% )15
Note: All immigrants receive a tag to provide immunity to mortality and emigration.
Immigrant initializationAssign current age (A)Assign current age (A)Assign current age (A)Assign end age (eA), then calculate A and iA by subtracting remaining time steps and 5, respectively.
Period advanceNAt + 1NAt + 1
Interval advanceT + 5, remove any tags that identify newborn cohorts or immigrant agents
Table 4
Average Percent Divergence from UN Target Values in 5-year Models
BirthsDeathsTotal Population
CountryPeriodIntuitiveSF*IntuitiveSF*IntuitiveSF*
Norway1950-20992.70.00.80.01.50.0
1950-19992.60.00.10.00.50.0
2000-20492.70.00.20.01.50.0
2050-20992.70.01.50.02.00.0
USA1950-20993.70.01.80.02.70.0
1950-19993.60.00.20.11.30.0
2000-20494.20.01.00.03.00.0
2050-20993.40.03.10.03.20.0
India1950-20999.30.07.00.08.90.0
1950-19997.30.03.20.04.00.0
2000-204911.20.05.00.09.30.0
2050-20999.40.010.10.010.70.0
  1. *

    Split Fertility Design

Table 5
Average Percent Divergence from UN Target Values in 1-year Models
BirthsDeathsTotal Population
CountryPeriodIntuitiveSF*IntuitiveSF*IntuitiveSF*
Norway1950-20991.60.33.30.52.10.3
1950-19991.00.33.30.91.20.1
2000-20490.80.33.80.71.30.2
2050-20992.90.13.00.23.30.4
USA1950-20992.00.52.20.63.10.3
1950-19990.90.63.60.51.40.2
2000-20491.60.23.10.62.70.3
2050-20993.50.50.80.64.50.2
India1950-209919.50.320.80.620.70.3
1950-19998.10.318.31.111.10.3
2000-204923.80.414.60.521.20.3
2050-209929.40.325.70.424.40.4
  1. *

    Split Fertility Design

Data and code availability

The authors used publicly available data for this analysis. The 2019 Revision of the World Population Prospects are available as excel files to download from the archive section of the United Nations website: https://population.un.org/wpp/Download/Archive/Standard/. The data are available without registration and are not proprietary.

The model was built in AnyLogic version 8.7.6. AnyLogic can be downloaded and installed at: https://www.anylogic.com/downloads/. There is a free personal learning edition that can load the model and run the model at low agent counts. An AnyLogic Professional or Researcher license is required to run the model at the scale ran in the paper (100K initial agents).

Results from the AnyLogic models were compiled and made into graphs using R version 4.1.2. R is open-source and can be downloaded and installed at: https://www.r-project.org/.

The AnyLogic and R code for this paper are available as replication materials on a GitHub repository accessible at: https://github.com/centerformindandculture/UN-CCM-2-MICROSIM. Others can re-use the code under a Creative Commons license, and we request our paper be cited when borrowing portions of the code.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)