Hi!
One more question. How do I make one line thicker in a graph? I want the line for the fully adjusted model (number 2, red) to be thicker than the rest. Can I also manually change color of the lines?
regress HGB Total_CO2 if(SEX==0), robust
local betahat1= `:di %5.3f _b[Total_CO2]'
quietly eststo crude: margins, at(Total_CO2=(0(1)10)) post
regress HGB Total_CO2 AGE DATUM i.SEASON if(SEX==0), robust
local betahat2= `:di %5.3f _b[Total_CO2]'
quietly eststo full: margins, at(Total_CO2=(0(1)10)) post
regress HGB Total_CO2 AGE DATUM i.SEASON if(SEX==0 & HT115==0), robust
local betahat3= `:di %5.3f _b[Total_CO2]'
quietly eststo dietchange: margins, at(Total_CO2=(0(1)10)) post
regress HGB Total_CO2 AGE DATUM i.SEASON if(SEX==0 & ENER_RAP_ADEQUATE==1), robust
local betahat4= `:di %5.3f _b[Total_CO2]'
quietly eststo misreporters: margins, at(Total_CO2=(0(1)10)) post
coefplot crude full dietchange misreporters, at title(Haemoglobin (women)) ytitle(Haemoglobin) xtitle(CO{sub:2} (kg/day)) leg(order(1 "Crude model (B{sub:1}=`betahat1')"2 "Fully adjusted model (B{sub:2}=`betahat2') "3 "Excluding recent diet changers (B{sub:3}=`betahat3')"4 "Excluding potential mis-reporters (B{sub:4}=`betahat4')") pos(6)) recast(line) noci
One more question. How do I make one line thicker in a graph? I want the line for the fully adjusted model (number 2, red) to be thicker than the rest. Can I also manually change color of the lines?
regress HGB Total_CO2 if(SEX==0), robust
local betahat1= `:di %5.3f _b[Total_CO2]'
quietly eststo crude: margins, at(Total_CO2=(0(1)10)) post
regress HGB Total_CO2 AGE DATUM i.SEASON if(SEX==0), robust
local betahat2= `:di %5.3f _b[Total_CO2]'
quietly eststo full: margins, at(Total_CO2=(0(1)10)) post
regress HGB Total_CO2 AGE DATUM i.SEASON if(SEX==0 & HT115==0), robust
local betahat3= `:di %5.3f _b[Total_CO2]'
quietly eststo dietchange: margins, at(Total_CO2=(0(1)10)) post
regress HGB Total_CO2 AGE DATUM i.SEASON if(SEX==0 & ENER_RAP_ADEQUATE==1), robust
local betahat4= `:di %5.3f _b[Total_CO2]'
quietly eststo misreporters: margins, at(Total_CO2=(0(1)10)) post
coefplot crude full dietchange misreporters, at title(Haemoglobin (women)) ytitle(Haemoglobin) xtitle(CO{sub:2} (kg/day)) leg(order(1 "Crude model (B{sub:1}=`betahat1')"2 "Fully adjusted model (B{sub:2}=`betahat2') "3 "Excluding recent diet changers (B{sub:3}=`betahat3')"4 "Excluding potential mis-reporters (B{sub:4}=`betahat4')") pos(6)) recast(line) noci
Comment