Create census.R