
360 Statistical Programming in SAS
Display 10.6 PROC PRINT output with the numerical integral of P(-1.645<Z<1.645)
Trapezoidal Rule Area Estimate for P(-1.645<Z<1.645)
(based on 50 equally spaced points)
Obs low high incr area_est
1
-1.645 1.645 0.067143 0.89990
If you repeatedly run this program and change &npts to different values, then you get
the output in Display 10.7.
Display 10.7 Output written to the file est.out
est. P(-1.645 < Z < 1.645) =0.8999027424 (based on 50 points)
est. P(-1.645 < Z < 1.645) =0.8994989035 (based on 25 points)
est. P(-1.645 < Z < 1.645) =0.8984685898 (based on 15 points)
est. P(-1.645 < Z < 1.645) ...