/* module 4: problems 1, 4ab, 5 */

/* 4.1 */
filename well 'Well#1.dat';
data one;
  infile well;
  input monthday $ 1-5;
  month = substr(monthday,1,3);
  day = substr(monthday,4,2);
run;
proc print data=one;
run;

/* 4.4ab */
data two;
  input x;
  cdf = probbnml(.4,5,x);
  pmf = cdf - lag(cdf);
  if x=0 then pmf = cdf;
  cards;
    0
    1
    2
    3
    4
    5
  ;
run;
proc print data=two;
run;

/* 4.5 */
* ammendment use <= instead of < in requested distributions;
data three;
  mean = 12.6;
  sd = 2.3;
  x1 = (10 - mean)/sd;
  d1 = probnorm(x1);
  x2 = (15 - mean)/sd;
  x3 = (7.6 - mean)/sd;
  d2 = probnorm(x2) - probnorm(x3);
run;
proc print data=three;
*  var d1 d2;
run;