gmacro 
test2p2
set 'terminal' c100;
nobs 2.
end
copy c100 k89 k88
let k88=k88+2
let k75=mean(c1)
let k76=n(c1)
let k77=mean(c2)
let k78=n(c2)
let k87=(k75-k77-k89)/sqrt(k75*(1-k75)/k76 + k77*(1-k77)/k78)
cdf k87 k86
let k85=1-k86
let k83=absolute(k87)
cdf k83 k82
let k84=2*(1-k82)
copy k89 c98
copy k87 c97
copy k75 c92
copy k77 c91
copy k86 k84 k85 c96
set c93;
format (a6).
   <
 not=
   >
end
copy c93 c96 c93 c96;
use k88.
print c98 c93 c98 c92 c91 c97 c96;
format ('H0: p1 - p2 =',f7.4,//'HA: p1 - p2',a6,f7.4,//,&
'p1-hat =',f9.5,//,'p2-hat =',f9.5,//,'z =',f9.5,//,'p-value =',f9.5).
erase k80-k99 c70-c100
endmacro