Update population.R