/* this do file conduct the ROC analysis for all potential obs, taking in info on proxy type 

*/
log using roc_all_2012.log, replace
quietly include ../../../fem_env.do

local maxrep 500
local minyr 2004
local maxyr 2014

local scenarios imp_totcog

foreach scn of local scenarios {

    * append all of the simulations
    forvalues rep = 1/`maxrep' {
        append using ../../../output/`scn'/detailed_output/y2012_rep`rep'.dta
    }
    gen totcog_si = totcog - 1

    * one observation per person
    bys hhidpn mcrep (year): keep if _n == _N
    
    gen demented = (totcog_si <= 6 ) if !missing(totcog_si)
    gen mci = (totcog_si <= 11) if !missing(totcog_si)
    
    foreach var in demented mci {
        rename `var' `var'_fem
    }
    
    keep hhidpn demented_fem mci_fem mcrep
    sum demented_fem mci_fem, de
    collapse demented_fem mci_fem, by(hhidpn)
    sum demented_fem mci_fem, de
    
    tempfile fem
    save `fem'
    
    /*all case*/
    use ../../../input_data/hrs_selected.dta, clear
    keep if wave==11
    keep if iwstat==1
    keep if age>=59
    merge m:1 hhidpn using `fem'
    keep if _m==3
    gen demented = (totcog <= 6)&!missing(totcog)|missing(totcog)&proxy_limitcog==1
    gen mci = (totcog <= 11)&!missing(totcog)|missing(totcog)&(proxy_limitcog==1|proxy_limitcog==2)

    di "ALL"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_mci, replace) title("`label1'") scheme(s1mono)

    preserve
    keep if male==1
    di "MALE"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_male_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_male_mci, replace) title("`label1'") scheme(s1mono)

    restore
    preserve
    keep if male==0
    di "FEMALE"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_female_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_female_mci, replace) title("`label1'") scheme(s1mono)

    restore
    preserve
    keep if hsless==1
    di "HSLESS"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_hsless_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_hsless_mci, replace) title("`label1'") scheme(s1mono)
    restore

    preserve
    keep if hsless==0
    di "NO HSLESS"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_nonhsless_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_nonhsless_mci, replace) title("`label1'") scheme(s1mono)
    restore

    preserve
    keep if college==1
    di "COLLEGE"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_college_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_college_mci, replace) title("`label1'") scheme(s1mono)
    restore

    preserve
    keep if college==0
    di "NO COLLEGE"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_noncollege_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_noncollege_mci, replace) title("`label1'") scheme(s1mono)
    restore

    preserve
    keep if age >=83
    di "STARTING AGE >= 75"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_age75_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_age75_mci, replace) title("`label1'") scheme(s1mono)
    restore

    preserve
    keep if age >= 78
    di "STARTING AGE >= 70"
    roctab demented demented_fem
    roctab demented demented_fem, graph saving(roc_plots_2012/all/`scn'_age70_demented, replace) title("`label1'") scheme(s1mono)

    roctab mci mci_fem
    roctab mci mci_fem, graph saving(roc_plots_2012/all/`scn'_age70_mci, replace) title("`label1'") scheme(s1mono)
    restore

}
    

log close