log using "roc_mortality_5fcv.log", replace

***
* product: none
* generate ROC curves for prediction of cognitive status
* indata: mortality_confirmed_2006_5fcv.dta; mortality_hrs_confirmed.dta
* outdata: none
* exhibits: Table5; Graph2
* dependent files: "../fem/mortality_confirmed_2006_5fcv.dta"; "../hrs/mortality_hrs_confirmed.dta"
***

use "../fem/mortality_confirmed_2006_5fcv.dta", clear
gen fem16_de_or_died_w_de = status2016==1|status2016==5
gen fem16_mci_or_died_w_de = (status2016==1)|(status2016==2)|(status2016==5)
gen fem16_mci_or_died_w_mci_or_worse = (status2016==1)|(status2016==2)|(dead_w_mci_or_worse_16==1)
gen fem_dementia = (status2016==1)
gen fem_mci_or_worse = (status2016==1)|(status2016==2)
collapse fem16_de_or_died_w_de fem16_mci_or_died_w_de fem_dementia fem_mci_or_worse fem16_mci_or_died_w_mci_or_worse, by(hhidpn)
keep fem16_de_or_died_w_de fem16_mci_or_died_w_de fem16_mci_or_died_w_mci_or_worse fem_dementia fem_mci_or_worse hhidpn

tempfile fem
save `fem'

use "../hrs/mortality_hrs_confirmed.dta", clear
gen hrs16_de_or_died_w_de = (status13_ab==1)|(status13_ab==5)
gen hrs16_mci_or_died_w_de = (status13_ab==1)|(status13_ab==2)|(status13_ab==5)
gen hrs16_mci_or_died_w_mci_or_worse = (status13_ab==1)|(status13_ab==2)|(dead_w_mci_or_worse_2016==1)
gen hrs_dementia = (status13_ab==1)
gen hrs_mci_or_worse = (status13_ab==1)|(status13_ab==2)
merge 1:1 hhidpn using `fem'
keep if _m==3
keep fem16_de_or_died_w_de fem16_mci_or_died_w_de fem16_mci_or_died_w_mci_or_worse hrs16_mci_or_died_w_mci_or_worse hrs16_de_or_died_w_de hrs16_mci_or_died_w_de hhidpn confirmed_cogcat_hrs2006_nomi hrs_mci_or_worse hrs_dementia fem_mci_or_worse fem_dementia college8 hsless8 age8 diabe8 stroke8 hibpe8 hearte8 

di "ALL - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/all_mortality, replace) title("ALL - dementia or dead with dementia") scheme(s1mono)

di "ALL - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/all_de, replace) title("ALL - dementia (conditional on alive)") scheme(s1mono)

di "ALL - MCI or worse (conditional on alive)"
roctab hrs_mci_or_worse fem_mci_or_worse
roctab hrs_mci_or_worse fem_mci_or_worse, graph saving(roc_plots_2016/5fcv_mortality/all_mci_or_worse, replace) title("ALL - MCI or worse (conditional on alive)") scheme(s1mono)

di "ALL - MCI or worse or dead with MCI or worse"
roctab hrs16_mci_or_died_w_mci_or_worse fem16_mci_or_died_w_mci_or_worse
roctab hrs16_mci_or_died_w_mci_or_worse fem16_mci_or_died_w_mci_or_worse, graph saving(roc_plots_2016/5fcv_mortality/all_mortality_mci_or_worse, replace) title("ALL - MCI or worse or dead with MCI or worse") scheme(s1mono)

preserve
keep if hsless8==1
di "HSLESS - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/hsless_mortality, replace) title("HSLESS - dementia or dead with dementia") scheme(s1mono)
di "HSLESS - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/hsless_de, replace) title("hsless - dementia (conditional on alive)") scheme(s1mono)

restore
preserve
keep if college8==1
di "COLLEGE - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/college_mortality, replace) title("COLLEGE - dementia or dead with dementia") scheme(s1mono)
di "COLLEGE - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/college_de, replace) title("college - dementia (conditional on alive)") scheme(s1mono)

restore
preserve
keep if age8>=75
di "STARTING AGE >= 75 - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/age75p_mortality, replace) title("2006 age 75p - dementia or dead with dementia") scheme(s1mono)
di "starting age >= 75 - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/age75p_de, replace) title("age75p - dementia (conditional on alive)") scheme(s1mono)

restore
preserve
keep if diabe8==1
di "diabetes - dementia or dead with dementia "
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/diabe_mortality, replace) title("diabetes - dementia or dead with dementia") scheme(s1mono)
di "diabetes - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/diabe_de, replace) title("diabetes - dementia (conditional on alive)") scheme(s1mono)

restore
preserve
keep if stroke8==1
di "stroke - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/stroke_mortality, replace) title("stroke - dementia or dead with dementia") scheme(s1mono)
di "stroke - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/stroke_de, replace) title("stroke - dementia (conditional on alive)") scheme(s1mono)

restore
preserve
keep if hibpe8==1
di "hypertension - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/hibpe_mortality, replace) title("hypertension - dementia or dead with dementia") scheme(s1mono)
di "hypertension - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/hibpe_de, replace) title("hypertension - dementia (conditional on alive)") scheme(s1mono)

restore
preserve
keep if hearte8==1
di "heart disease - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/hearte_mortality, replace) title("heart disease - dementia or dead with dementia") scheme(s1mono)
di "heart disease - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/hearte_de, replace) title("heart disease - dementia (conditional on alive)") scheme(s1mono)

keep if confirmed_cogcat_hrs2006_nomi==2
di "MCI2006 - dementia or dead with dementia"
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de
roctab hrs16_de_or_died_w_de fem16_de_or_died_w_de, graph saving(roc_plots_2016/5fcv_mortality/mci2006_mortality, replace) title("MCI2006 - dementia or dead with dementia") scheme(s1mono)


di "MCI2006 - dementia (conditional on alive)"
roctab hrs_dementia fem_dementia
roctab hrs_dementia fem_dementia, graph saving(roc_plots_2016/5fcv_mortality/mci2006_de, replace) title("MCI2006 - dementia (conditional on alive)") scheme(s1mono)

log close