sysuse auto, clear
qui regress mpg i.foreign
local df = e(df_r)
local t = round(((_b[1.foreign]/_se[1.foreign])), .01)
local p = round((2 * ttail(`df', abs(`t'))), .001)
margins, at(foreign = (0 1))
marginsplot, recast(bar) ytitle("Mileage") xtitle("Car type") ///
plotopts(fcolor(gs14) lcolor(black)) ///
title("{bf:A} Mileage", justification(left) bexpand span) ///
addplot(scatteri 28.95 0 28.95 1, recast(line) lwidth(medium) lpattern(solid) ///
text(29.5 .5 "{it:t}(`df') = `t', {it:p} = `p'") ///
xlabel(0 1) ///
|| scatteri 28.95 0 28.05 0, recast(line) lwidth(medium) lpattern(solid) ///
|| scatteri 28.95 1 28.05 1, recast(line) lwidth(medium) lpattern(solid)) ///
legend(off) ///
name(figureA, replace)
qui regress weight i.foreign
local df = e(df_r)
local t = round(((_b[1.foreign]/_se[1.foreign])), .01)
local double p = round((2 * ttail(`df', abs(`t'))), .0000001)
margins, at(foreign = (0 1))
marginsplot, recast(bar) ytitle("Weight") xtitle("Car type") ///
plotopts(fcolor(gs14) lcolor(black)) ///
title("{bf:B} Weight", justification(left) bexpand span) ///
addplot(scatteri 3700 0 3700 1, recast(line) lwidth(medium) lpattern(solid) ///
text(3760 .5 "{it:t}(`df') = `t', {it:p} = `p'") ///
xlabel(0 1) ///
|| scatteri 3700 0 3600 0, recast(line) lwidth(medium) lpattern(solid) ///
|| scatteri 3700 1 3600 1, recast(line) lwidth(medium) lpattern(solid)) ///
legend(off) ///
name(figureB, replace)
graph combine figureA figureB, col(2) name(figure, replace)
