STATA = ../../run.stata15.sh

# number of simulation reps
NREPS = 1000

# scenarios with no migration
NOMIGSCENARIOS = mortval_cohort_mortadj2016 mortval_cohort_minimal_mortadj2016 mortval_cohort_nomortadj mortval_cohort_minimal_nomortadj mortval_cohort_mortadj2004 mortval_cohort_minimal_mortadj2004

# scenarios with migration
MIGSCENARIOS = mortval_cohort_nomortadj_mig2004 mortval_cohort_minimal_nomortadj_mig2004 mortval_cohort_mortadj2004_mig2004 mortval_cohort_minimal_mortadj2004_mig2004

data/%_mortality_rates.dta : fem_mortality_rates.do $(wildcard ../../output/%/*_summary.dta)
	SCEN=$* $(STATA) fem_mortality_rates.do
	
mr : $(addprefix data/, $(addsuffix _mortality_rates.dta, $(NOMIGSCENARIOS)))

data/%_cdc_smr.dta : data/%_mortality_rates.dta fem_cdc_smr.do
	SCEN=$* $(STATA) fem_cdc_smr.do
	

data/%_hrs_smr.dta : data/%_mortality_rates.dta fem_hrs_smr.do
	SCEN=$* $(STATA) fem_hrs_smr.do
	
data/%_cdc_bellmiller_mortprobdiff.dta : bellmillermortprob.do data/%_cdc_smr.dta
	SCEN=$* $(STATA) bellmillermortprob.do

plots : plot_mortality.R $(addprefix data/, $(addsuffix _cdc_smr.dta, $(NOMIGSCENARIOS)) $(addsuffix _hrs_smr.dta, $(NOMIGSCENARIOS)) $(addsuffix _cdc_bellmiller_mortprobdiff.dta, $(NOMIGSCENARIOS)))
	Rscript plot_mortality.R

data/%_e65.dta : $(wildcard ../../output/%/*_by_rep.dta) fem_e65.do
	SCEN=$* $(STATA) fem_e65.do

SSA_e65.dta : SSA_e65.xlsx
	st SSA_e65.xlsx SSA_e65.dta -y

e65: $(addprefix data/, $(addsuffix _e65.dta, $(NOMIGSCENARIOS))) SSA_e65.dta fem_ssa_e65.do
	$(STATA) fem_ssa_e65.do

data/census_proj2004.dta : census_proj2004.do
	$(STATA) census_proj2004.do

data/census_popsizeest.dta : census_popsizeest.do
	$(STATA) census_popsizeest.do

data/%_census_popsize.dta : fem_census_popsize.do $(wildcard ../../output/%/*_by_rep.dta) data/census_proj2004.dta data/census_popsizeest.dta
	SCEN=$* NREPS=$(NREPS) $(STATA) fem_census_popsize.do

popsize : $(addprefix data/, $(addsuffix _census_popsize.dta, $(MIGSCENARIOS)))

data/cdc_survival.dta : cdc_survival.do 
	$(STATA) cdc_survival.do

data/bellmiller_survival.dta : bellmiller_survival.do 
	$(STATA) bellmiller_survival.do

data/%_cdc_bellmiller_survival.dta : fem_cdc_bellmiller_survival.do $(wildcard ../../output/%/*_summary.dta) data/cdc_survival.dta data/bellmiller_survival.dta
	SCEN=$* $(STATA) fem_cdc_bellmiller_survival.do

plot_survival : $(addprefix data/, $(addsuffix _cdc_bellmiller_survival.dta, $(NOMIGSCENARIOS))) plot_survival.R
	Rscript plot_survival.R