options nocenter; data hyper; N = 41; *<--TOTAL NUMBER OF OBSERVATIONS; M = 36; *<--TOTAL NUMBER OF SUCCESSES; k = 23; *<--TOTAL NUMBER OF GROUP 1 OBSERVATIONS; * Based on the fixed margins of Table 2.17 the possible values of n11 are 18, 19, ..., 23; do n11 = 18 to 23; prob=PDF('hyper',n11,41,36,23); output; end; run; proc print data=hyper; var n11 prob; run;