%let p=0.5; %let n=2; data Bino_PMF; do k=0 to &n; PMF=pdf('Binomial', k, &p, &n); output; end; run; /* Plot PMF */ title "Binomial PMF with p=&p and n=&n"; proc sgplot data=Bino_PMF noautolegend; vbar k / response = PMF ; xaxis values=(0 to 2) label='NUMBER OF SUCCESSES' labelattrs=(size=12 weight=Bold); yaxis display=(nolabel); ***keylegend / position=NE location=inside across=1 noborder valueattrs=(Size=12 Weight=Bold); run; title; %let p=0.5; %let n=10; data Bino2_PMF; do k=0 to &n; PMF=pdf('Binomial', k, &p, &n); output; end; run; /* Plot PMF */ title "Binomial PMF with p=&p and n=&n"; proc sgplot data=Bino2_PMF noautolegend; vbar k / response = PMF ; xaxis values=(0 to 10) label='NUMBER OF SUCCESSES' labelattrs=(size=12 weight=Bold); yaxis display=(nolabel); ***keylegend / position=NE location=inside across=1 noborder valueattrs=(Size=12 Weight=Bold); run; title; %let p=0.5; %let n=200; data Bino3_PMF; do k=0 to &n; PMF=pdf('Binomial', k, &p, &n); output; end; run; /* Plot PMF */ ods pdf file='/home/mstear0/Distributions/Binomial3.pdf' dpi=600; ***odstitle=NONE; title "Binomial PMF with p=&p and n=&n"; proc sgplot data=Bino3_PMF noautolegend; vbar k / response = PMF ; xaxis values=(0 to 200) label='NUMBER OF SUCCESSES' labelattrs=(size=12 weight=Bold) Fitpolicy=THIN; yaxis display=(nolabel); ***keylegend / position=NE location=inside across=1 noborder valueattrs=(Size=12 Weight=Bold); run; title;