#delimit;
set more off;

estsimp mlogit vote limimp age yrseduc pid;

quietly generate pbush = .;
quietly generate pblo = .;
quietly generate pbhi = .;
quietly generate pclin = .;
quietly generate pclo = .;
quietly generate pchi = .;
quietly generate pperot = .;
quietly generate pplo = .;
quietly generate pphi = .;
quietly generate gryred = .;

local a= 18;
while `a' < 92 {;
setx limimp 1 age `a' yrseduc median pid median;
quietly replace gryred = `a' if _n== `a';
quietly simqi, prval(1 2 3) genpr(pb pc pp);
quietly egen mpb = mean(pb);
quietly egen mpc = mean(pc);
quietly egen mpp = mean(pp);
quietly replace pbush = mpb if _n== `a';
quietly replace pclin = mpc if _n==`a';
quietly replace pperot = mpp if _n==`a';
_pctile pb, p(2.5, 97.5);
quietly replace pblo = r(r1) if _n==`a';
quietly replace pbhi = r(r2) if _n==`a';
_pctile pc, p(2.5, 97.5);
quietly replace pclo = r(r1) if _n==`a';
quietly replace pchi = r(r2) if _n==`a';
_pctile pp, p(2.5, 97.5);
quietly replace pplo = r(r1) if _n==`a';
quietly replace pphi = r(r2) if _n==`a';

quietly drop pb pc pp mpb mpc mpp;
local a = `a' + 1;
};


label var pbush "Mean Bush";
label var pclin "Mean Clinton";
label var pperot "Mean Perot";
label var pblo "B 95% LL";
label var pbhi "B 95% UL";
label var pclo "C 95% LL";
label var pchi "C 95% UL";
label var pplo "P 95% LL";
label var pphi "P 95% UL";
label var gryred "Age";

graph twoway line pbush pclin pperot pblo pbhi pclo pchi pplo pphi gryred, 
lwidth(medthick medthick medthick thin thin thin thin thin thin) 
lcolor(black black black black black black black black black) 
lpattern(solid longdash_dot dash dot dot dot dot dot dot) 
ytitle("Probability of Vote");

drop b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 pbush pclin pperot pblo pbhi pclo pchi pplo pphi gryred;