new;
output file=c:\jae\ols\wamp2\wampler2.out  reset;

@Dataset Name:  Wampler-2 (wampler2.dat)
Data:          1 Response Variable (y)
               1 Predictor Variable (x)
               21 Observations
               Higher Level of Difficulty
Model:         Polynomial Class
               6 Parameters (B0,B1,...,B5)
     y = B0 + B1*x + B2*(x**2) + B3*(x**3)+ B4*(x**4) + B5*(x**5)
               Certified Regression Statistics

                                           Standard Deviation
     Parameter         Estimate               of Estimate
@
cbcse={        1.00000000000000         0.000000000000000,
               0.100000000000000        0.000000000000000,
               0.100000000000000E-01    0.000000000000000,
               0.100000000000000E-02    0.000000000000000,
               0.100000000000000E-03    0.000000000000000,
               0.100000000000000E-04    0.000000000000000};
cb=cbcse[.,1];
cstdb=cbcse[.,2];

cR2= 1.00000000000000;

@Data:          y       x@
yx={        1.00000    0,
            1.11111    1,
            1.24992    2,
            1.42753    3,
            1.65984    4,
            1.96875    5,
            2.38336    6,
            2.94117    7,
            3.68928    8,
            4.68559    9,
            6.00000   10,
            7.71561   11,
            9.92992   12,
           12.75603   13,
           16.32384   14,
           20.78125   15,
           26.29536   16,
           33.05367   17,
           41.26528   18,
           51.16209   19,
           63.00000   20};
y=yx[.,1];x=yx[.,2];
rhs=x~(x^2)~(x^3)~(x^4)~(x^5);
ans1=nistols(y,rhs,cb,cstdb,cr2);  
ans2=niolsqr2(y,rhs,cb,cstdb,cr2);  
format /ld  6,1;
"        " ans1;  
"   (          "  ans2   "   ) ";
"Wampler2";
end;
