gmacro 
test2p1
set 'terminal' c100;
nobs 1.
end
copy c100 k88
let k88=k88+2
let k75=mean(c1)
let k76=n(c1)
let k77=mean(c2)
let k78=n(c2)
stack c1-c2 c99
let k79=mean(c99)
let k87=(k75-k77)/sqrt(k79*(1-k79)*(1/k76 + 1/k78))
cdf k87 k86
let k85=1-k86
let k83=absolute(k87)
cdf k83 k82
let k84=2*(1-k82)
copy k88 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.
let k93=0
copy k93 c98
print c98 c93 c98 c92 c91 c97 c96;
format ('H0:p1-p2 = ',i1,//,'HA:p1-p2',a6,i1,//,'p1-hat =',f9.5,//,&
'p2-hat =',f9.5,//,'z =',f9.5,//,'p-value =',f9.5).
erase k80-k99 c70-c100
endmacro