Added example for population forecast: test_population.R