title 'PMF: Y~bin(n=10,p=0.26)';
data args;
p=.26;
n=10;
do y=0 to n;  
output; end;
run; 
data bin1;
set args;
py=PDF('BINOMIAL',y,p,n);
run;
proc print data=bin1;
var y py;
run;
proc sgplot data=bin1;
needle x=y y=py;
xaxis values=(0 to 10);
run; quit;
title;

title 'CDF: Y~bin(n=10,p=0.26)';
data bin2;
set args;
cdfy=probbnml(p,n,y);
run;
proc print data=bin2 label;
var y cdfy;
label cdfy='py';
run;
proc sgplot data=bin2;
scatter x=y y=cdfy;
series x=y y=cdfy;
xaxis values=(0 to 10);
yaxis label='py';
run; quit;
title;

title 'Y~bin(n=20,p=0.26)';
data rbin;
do x=0 to 20;
 px=ranbin(7777,30,0.26);
 output;
end;
run;
proc sgplot data=rbin;
needle x=x y=px;
xaxis values=(0 to 20);
run; quit;
title;