gmacro
est2p
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=mean(c2)
let k83=n(c2)
let k84=k90*sqrt(k80*(1-k80)/k81 + k82*(1-k82)/k83)
let k85=k80-k82-k84
let k86=k80-k82+k84
copy k80 c89
copy k82 c90
copy k85 c91
copy k86 c92
print  c100 c89 c90 c91 c92;
format ('Confidence level =',I3,'%',//'p1-hat =',f9.5,//,'p2-hat =',f9.5,//,&
'Lower confidence limit =',f9.5,//,'Upper confidence limit =',f9.5).
erase k80-k90 c89-c100
endmacro