/* Example 3.4.1-----Kolmogorov's Test */
/* p=2, n=32, alpha=5% */
proc iml;
xx={1.4 1.4 1.3 1.5 1.4 1.7 1.4 1.5 1.4 1.5
    1.5 1.6 1.4 1.1 1.2 1.5 1.3 1.4 1.7 1.5
    1.7 1.5 1.0 1.7 1.9 1.6 1.6 1.5 1.4 1.6
    1.6 1.5 1.5 1.4 1.5 1.2 1.3 1.4 1.3 1.5
    1.3 1.3 1.3 1.6 1.9 1.4 1.6 1.4 1.5 1.4};
xb=J(1,50,1/50)*xx`;
s=sqrt((xx-xb)*(xx-xb)`/49);
n=50;
d1=1-probnorm((2-xb)/s);
d3=1-probnorm((xb-1)/s);
d2=0;
do i=1 to 200;
b= (xx-(1+i/200))<>0;
c=^(^b);
Fn=1-sum(c)/n;
di=abs(Fn-probnorm((1+i/200-xb)/s));
if di>d2 then d2=di;
end;
d=max(d1, d2, d3);
print d xb s;
quit;
             D        XB             S
     0.1533977     1.462      0.173664