data kaz;
N_yes=1162;
N_no=381;
mean_yes=0.129088;
mean_no=0.170604;
DEG_FD=N_Yes + N_No -2;
/*QC’ed
tValue=2.228;
DEG_FD=10;
*/
tValue=(Mean_NO-Mean_Yes)/(SQRT((Mean_NO*(1-Mean_NO) /
N_No )+Mean_Yes*(1-Mean_Yes)/ N_Yes));
/*2 tail test*/
P_value=(1-probt(abs(tValue),DEG_FD))*2;
run;