gmacro
estp
set 'terminal' c100;
nobs=1.
end
let k88 = (c100(1)/100)+(1-c100(1)/100)/2
invcdf k88 k90
let k80=mean(c1)
let k81=n(c1)
let k82=k90*sqrt((k80*(1-k80)/k81))
let k83=k80-k82
let k84=k80+k82
copy k80 c90
copy k83 c91
copy k84 c92
print c100 c90 c91 c92;
format ('Confidence level =',i3,'%',//,'p-hat =',f9.5,//,&
'Lower confidence limit =',f9.5,//,'Upper confidence limit =',f9.5).
erase k80-k90 c90-c100
endmacro