Scenario model of the future reformed population of the city of Zurich
This R code calculates population scenarios for the future reformed population of the city of Zurich. For data protection reasons the input and output data are not available here. The results of the calculations are population sizes by the ten church districts, years, and age classes (10 year classes).
Code to run all modules
- 00_run-all.r: run all modules
Calculate future rates
-
10_birth-fertility.r: future fertility rate
-
11_birth-rel.r: religion of mother and child
-
20_death.r: mortality rates and life expectancy by religion
-
30_immigration.r: future immigration rate
-
40_emigration.r: future emigration rate
-
50_conversion.r: future conversion rate
Core model
- 60_demography-popscen.r: combine demographic processes and population scenarios
Area correction factors
- 70_area.r: area factors
Functions
-
90_general.r: general functions
-
91_plot-functions.r: plot functions
-
92_constrained-regression.r: constrained regression
-
93_life-expectancy.r: functions to calculate life expectancy
We use an renv enabled RStudio project to manage the local environment with libraries in order to ensure reproducible results, correct working directories etc. You can use the included one or create your own.
When you pull this repository, make sure to restart your R session (this is done automatically if you open this as an RStudio project). R will read the project specific .Rprofile file and thereby activate renv.
You can also activate it manually by using renv::activate(). This might lead to the download of many (potentially older) packages.
Keep in mind that you will be working in this specific environment as long as you're not switching to another project or manually deactivate renv (renv::deactivate()).
If a new package is required, use renv::install("new package"). For updates, use renv::update("update_package").
Then add this to the lock file by running renv::snapshot() and commit the changes in renv.lock.
A good overview can be found here