%macro bin_probs(nobs); data bin_comp; do k=0 to &nobs; n_minus_k=&nobs-k; alpha0=2*probbnml(0.5,&nobs,n_minus_k); output; end; run; title 'Values of alpha0=2P(X<='&nobs'-k)'; title2'Where X has a b('&nobs',0.5) Distribution'; proc print noobs;run; %mend bin_probs; %bin_probs(11);