#delimit;
set more off;

estsimp ologit retNatEcon GWBtherm PartyID age yrseduc;

generate pworse = .;
generate psame = .;
generate pbetter = .;
generate gryred = .;

local a= 1;
while `a' < 18 {;
setx GWBtherm mean PartyID median age mean yrseduc `a';
quietly replace gryred = `a' if _n== `a';
simqi, prval(-1 0 1) genpr(plow pmed phigh);
quietly egen mplow = mean(plow);
quietly egen mpmed = mean(pmed);
quietly egen mphigh = mean(phigh);
quietly replace pworse = mplow if _n== `a';
quietly replace psame = mpmed if _n==`a';
quietly replace pbetter = mphigh if _n==`a';
quietly drop plow pmed phigh mplow mpmed mphigh;
local a = `a' + 1;
};


label var pworse "Worse";
label var psame "Same";
label var pbetter "Better";
label var gryred "Years of education";

graph twoway line pworse psame pbetter gryred, 
lwidth(medthick medthick medthick) lcolor(black black black) 
lpattern(dash solid dot) ytitle("Probability of Ret. Econ. Opinion");

drop b1 b2 b3 b4 b5 b6 pworse psame pbetter gryred;