One more question:
Can i include a histogram in a coefplot? How?
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
hist Total_CO2 if(SEX==0)
coefplot crude full, 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') ") pos(6)) recast(line) noci
Anna
Can i include a histogram in a coefplot? How?
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
hist Total_CO2 if(SEX==0)
coefplot crude full, 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') ") pos(6)) recast(line) noci
Anna