*
* Example from section 11.3.3 from pp 365-366
*
all 10000
compute alpha_1=2.0,se_a1=.25
compute alpha_2=0.5,se_a2=.08
*
* To get a Normal draw with mean alpha_1 and std dev se_a1, you can use
* alpha_1+%ran(se_a1), as the %ran(x) draws a Normal with mean 0 and standard
* deviation x.
*
set top    = alpha_1+%ran(se_a1)
set bottom = alpha_2+%ran(se_a2)
set ratio  = top/bottom
compute f=%fractiles(ratio,||.025,.975||)
disp "95% Interval" f(1) f(2)