log using logs/census_proj2004.log, replace

include ../../fem_env.do

/*** load 2004 census projections and collapse to age range of 2004 stock FEM cohort in each simulation year ***/
clear
insheet using ../../input_data/usproj2000-2050.csv

keep if racegroupcode=="RTT" & year >= 2004

forvalues i=6/106 {
	  local age = `i' - 6
	  rename v`i' age`age'
}

keep year sex age*
reshape long age, i(year sex) j(age_new)
rename age popsize
label var popsize "Census population size projection including migration"
rename age_new age

desc
tab year sex

drop if age < 51

gen stockagemin = 51 + year - 2004

replace popsize = 0.0 if age < stockagemin

collapse (sum) stockpopsize=popsize (first) stockagemin, by(year sex)

rename stockpopsize stockpopsize_census
keep year sex stockagemin stockpopsize_census

* drop odd-numbered years that are not simulated in FEM
drop if mod(year,2) == 1

* rescale census projection counts to 1000s in order to match FEM scale
replace stockpopsize_census = stockpopsize_census / 1000.0 

save data/census_proj2004.dta, replace