*** Predict the value of triglyceride in NHANES
*** And use the same prediction model to estimate the TG value in HRS
**  Based on the equation from NHANES  [LDL-cholesterol] = [total cholesterol] – [HDL-cholesterol] – [triglycerides/5]
set more off

quietly include "../../../fem_env.do"


use $outdata/nhanes_analytic.dta, clear 

keep if cohort>=2005 & age >=51


mkspline age_l60 60 age_6070 70 age_7080 80 age80p = age

*** interaction between tx and diseases

gen rxchol_stroke = rxchol*stroke
gen rxchol_hearte = rxchol*hearte
gen rxchol_diabe  = rxchol*diabe
gen rxchol_cancre  = rxchol*cancre
gen tchol_hdl      = tchol*hdl

*** interaction between hdl tchol and health conditions

gen tchol_hearte = tchol*hearte
gen tchol_storke = tchol*stroke
gen tchol_diabe = tchol*diabe
gen tchol_hibpe = tchol*hibpe
gen tchol_cancre = tchol*cancre

gen hdl_hearte = hdl*hearte
gen hdl_storke = hdl*stroke
gen hdl_diabe = hdl*diabe
gen hdl_hibpe = hdl*hibpe
gen hdl_cancre = hdl*cancre

gen tcholmhdl = tchol - hdl 

corr tcholmhdl ldl



regress ldl tchol hdl c.age/*age_l60 age_6070 age_7080 age80p*/ male black hispan male_black male_hispan smoken logbmi hearte stroke diabe hibpe cancre rxchol rxchol_stroke rxchol_hearte diabe  ///
        [aw=ldl_weight]

regress ldl tchol hdl tchol_hdl age_l60 age_6070 age_7080 age80p male black hispan male_black male_hispan smoken logbmi hearte stroke diabe hibpe cancre rxchol  ///
        rxchol_stroke rxchol_hearte rxchol_diabe rxchol_cancre tchol_hearte tchol_storke tchol_diabe tchol_hibpe tchol_cancre hdl_hearte hdl_storke hdl_diabe   ///
        hdl_hibpe hdl_cancre   [aw=ldl_weight]



ereturn list

est save $local_path/Estimates/nhanes_biomarkers/nhanes_ldl.ster, replace


regress ldl tchol hdl tchol_hdl age_l60 age_6070 age_7080 age80p male black hispan male_black male_hispan [aw=ldl_weight]



ereturn list

est save $local_path/Estimates/nhanes_biomarkers/nhanes_ldl_minimal.ster, replace