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

## Years for which CDC life tables are available
DATAYEARS = 2001  2002  2003  2004  2005  2006  2007  2008  2009  2010  2011  2012  2013  2014  2015

# create HRS analytic files
data/hrs_analytic_%.dta : gen_hrs_analytic_file.do
	VERSION=$* $(STATA) gen_hrs_analytic_file.do 

hrs_analytic_all : data/hrs_analytic_P.dta data/hrs_analytic_2014v3.dta

# Replicate Weir mortality calculation

## compute actual and expected HRS death rates for each year and demographic group
data/weir_mortality_repl%.dta : data/hrs_analytic_P.dta ../../base_data/cdc_nvsr_lifetable%.dta weir_mortality_repl_yr.do
	YEAR=$* $(STATA) weir_mortality_repl_yr.do

WEIRDTAYRS = $(addprefix data/weir_mortality_repl, $(addsuffix .dta, $(DATAYEARS)))
weir_repl : $(WEIRDTAYRS)



# Mortality calculation using only HRS death dates (no NDI) dropping attritors

## compute actual and expected HRS death rates for each year and demographic group
data/hrsdeath_mortality%.dta : data/hrs_analytic_2014v3.dta ../../base_data/cdc_nvsr_lifetable%.dta hrsdeath_mortality_yr.do
	YEAR=$* $(STATA) hrsdeath_mortality_yr.do

HRSDEATHDTAYRS = $(addprefix data/hrsdeath_mortality, $(addsuffix .dta, $(DATAYEARS)))
hrsdeath : $(HRSDEATHDTAYRS)



# Mortality calculation using only HRS death dates (no NDI) and HRS weights

## years where nursing home weights are available in the analytic file
NHMYRS = 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012

## compute actual and expected HRS death rates for each year and demographic group
data/weighted_hrsdeath_mortality%.dta : data/hrs_analytic_2014v3.dta ../../base_data/cdc_nvsr_lifetable%.dta weighted_hrsdeath_mortality_yr.do
	YEAR=$* $(STATA) weighted_hrsdeath_mortality_yr.do

WTDHRSDEATHDTAYRS = $(addprefix data/weighted_hrsdeath_mortality, $(addsuffix .dta, $(NHMYRS)))
wtdhrsdeath : $(WTDHRSDEATHDTAYRS)




# aggregate death rates across years and demographic groups
FORMATTED_LIFE_TABLES = $(addprefix ../../base_data, $(addprefix cdc_nvsr_lifetable, $(CDCYRDTA)))
data/mortality_ratios_allyrs.dta : aggregate_hrsmort.do $(WEIRDTAYRS) $(HRSDEATHDTAYRS) $(WTDHRSDEATHDTAYRS) $(FORMATTED_LIFE_TABLES)
	$(STATA) aggregate_hrsmort.do



weighted_hrs_mortality%.dta : hrs_analytic.dta weighted_mortality_rates_yr.do
	YEAR=$* $(STATA) weighted_mortality_rates_yr.do




qx_compare.dta : weighted_hrs_mortality_compare.do $(addprefix weighted_hrs_mortality, $(addsuffix .dta, $(NHMYRS))) formatted_life_tables_dta
	$(STATA) weighted_hrs_mortality_compare.do
	
%.Rdata : %.dta
	st $*.dta $*.Rdata -y	
	
qx_compare : qx_compare.Rdata plot_qx_compare.R
	Rscript plot_qx_compare.R