OPTIONS LS=64 PS=66 NONUMBER NODATE;

************************************;
*** Generate the pmf and cdf for ***;
*** a binomial random variable   ***;
************************************;
DATA binom;
     RETAIN tmp 0;
     N = 18;        *** N = number of trials ***;
     P = .4;        *** P = prob of success  ***;
     DO X = 0 TO N;
        cum_prob = probbnml(P,N,X);
        prob = cum_prob - TMP;
        tmp = cum_prob;
        OUTPUT;
     END;
PROC PRINT DATA=binom;
     ID N P; 
     VAR X prob cum_prob;
     FORMAT prob 11.10 cum_prob 11.10;
TITLE 'BINOMIAL DISTRIBUTION PROBABILITIES';
RUN;

**************************************;
*** Generate the pmf and cdf for a ***;
*** hypergeometric random variable ***;
**************************************;
DATA hyper;
     RETAIN tmp 0;
     nn= 18;                *** SAMPLE SIZE ;
     N = 49;                *** POPULATION TOTAL ;
     m = 22;                *** GROUP 1 TOTAL ;
     _max = MAX(0,m+nn-N);  *** LOWER BOUND FOR X;
     _min = MIN(m,nn);      *** UPPER BOUND FOR X;
     DO X = _max TO _min;
        cum_prob = PROBHYPR(N,m,nn,X);
        prob = cum_prob - tmp;
        tmp = cum_prob;
        OUTPUT;
     END;
PROC PRINT DATA=hyper;
     ID N m nn; 
     VAR X prob cum_prob;
     FORMAT prob 11.10 cum_prob 11.10;
TITLE 'HYPERGEOMETRIC DISTRIBUTION PROBABILITIES';

RUN;