### Use the last file as the target for now. Should split this file up at some point
$(ESTIMHRS)_extval/estimatesHRS.xls: $(COMMON) fem_env.do FEM_Stata/Estimation/init_transition.do $(TRANSDATA) $(UTILITIES)/takestring.ado FEM_Stata/Estimation/hrs_covariate_definitionsHRS.do $(DATADIR)/crossvalidation.dta FEM_Stata/Estimation/define_modelsHRS.do 
	cd FEM_Stata/Estimation && USERVARS=$(USERVARS) DATAIN=$(TRANSDATA) suffix=HRS EXTVAL=1 $(STATA) init_transition.do
	
$(ESTIMHRS)_extval/ghregHRS_extval.xls: $(COMMON) fem_env.do FEM_Stata/Estimation/ghreg_estimations.do $(TRANSDATA) $(UTILITIES)/takestring.ado FEM_Stata/Estimation/hrs_covariate_definitionsHRS.do FEM_Stata/Estimation/define_modelsHRS.do $(UTILITIES)/_gh.ado $(GHREG) $(UTILITIES)/_ggh.ado $(UTILITIES)/ghreg_p.ado $(DATADIR)/crossvalidation.dta $(UTILITIES)/opti_unc.mo 
	cd FEM_Stata/Estimation && USERVARS=$(USERVARS) DATAIN=$(TRANSDATA) suffix=HRS EXTVAL=1 $(STATA) ghreg_estimations.do

$(ESTIMHRS)_extval/qaly.xls: fem_env.do FEM_Stata/Estimation/qaly_estimations.do $(DATADIR)/hrs_analytic_recoded.dta $(DATADIR)/MEPS_EQ5D.dta $(DATADIR)/MEPS_cost_est.dta
	cd FEM_Stata/Estimation && EXTVAL=1 $(STATA) qaly_estimations.do

$(ESTIMATES)/minimal_extval/estimatesminimal.xls: $(COMMON) fem_env.do FEM_Stata/Estimation/init_transition.do $(TRANSDATA) $(UTILITIES)/takestring.ado FEM_Stata/Estimation/hrs_covariate_definitionsminimal.do $(DATADIR)/crossvalidation.dta FEM_Stata/Estimation/define_modelsminimal.do 
	cd FEM_Stata/Estimation && USERVARS=$(USERVARS) DATAIN=$(TRANSDATA) suffix=minimal EXTVAL=1 $(STATA) init_transition.do
	
extval_estimation: 	$(ESTIMHRS)_extval/estimatesHRS.xls $(ESTIMHRS)_extval/ghregHRS_extval.xls $(ESTIMHRS)_extval/qaly.xls $(ESTIMATES)/minimal_extval/estimatesminimal.xls
 
extval_estimation_cpp: $(COMMON)
	cd FEM_Stata/Estimation && datain=$(ESTIMHRS)_extval dataout=$(ROOT)/FEM_CPP_settings/hrs_extval/models $(STATA) save_est_cpp.do
	cd FEM_Stata/Estimation && datain=$(ESTIMATES)/minimal_extval dataout=$(ROOT)/FEM_CPP_settings/minimal_extval/models $(STATA) save_est_cpp.do

mort_assumptions_ref2001: $(COMMON)
	cd FEM_Stata/Makedata/POPULATION && SSATABLE=SSA2004VA1 REFYEAR=2001 STARTYEAR=1990 $(STATA) mortality_assumptions_any_start_year.do
	cd FEM_Stata/Makedata/POPULATION && SSATABLE=SSA2016VA1 REFYEAR=2001 STARTYEAR=1990 $(STATA) mortality_assumptions_any_start_year.do

$(DATADIR)/usproj2000-2050.csv : $(CENSUS)/2004_Pop_Projection/usproj2000-2050.xls
	st $(CENSUS)/2004_Pop_Projection/usproj2000-2050.xls $(DATADIR)/usproj2000-2050_temp.csv -y && cd $(DATADIR) && awk ' NR > 2 ' usproj2000-2050_temp.csv > usproj2000-2050.csv && rm $(DATADIR)/usproj2000-2050_temp.csv

$(DATADIR)/population_projection_2004.dta: FEM_Stata/Makedata/POPULATION/census_pop_projection_2004.do  $(DATADIR)/usproj2000-2050.csv
	cd FEM_Stata/Makedata/POPULATION && $(STATA) census_pop_projection_2004.do

$(DATADIR)/all2004_pop_adjusted_extval.dta: $(COMMON) FEM_Stata/Makedata/HRS/all2004_weightadjust.do $(DATADIR)/population_projection_2004.dta $(DATADIR)/all2004.dta $(DATADIR)/death_counts.dta
	cd FEM_Stata/Makedata/HRS && DATAIN=$(DATADIR)/all2004.dta DATAOUT=$(DATADIR)/all2004_pop_adjusted_extval.dta EXTVAL=1 $(STATA) all2004_weightadjust.do
	
$(DATADIR)/stock_hrs_2004_extval.dta: $(COMMON) $(RESTIMATES)/fraime2004.ster $(RESTIMATES)/fraime_nonzero2004.ster $(RESTIMATES)/flogq2004.ster $(RESTIMATES)/rpia2004.ster FEM_Stata/Makedata/HRS/save_stock_hrs_2004_faked.do $(DATADIR)/all2004_pop_adjusted_extval.dta $(UTILITIES)/multiply_persons.ado expansion.makefile
	cd FEM_Stata/Makedata/HRS && EXPAND=$(EXPANSION_FACTOR) EXTVAL=1 $(STATA) save_stock_hrs_2004_faked.do