new;
output file=c:\jae\ols\wamp4\wampler4.out  reset;

@Dataset Name:  Wampler-4 (wampler4.dat)
Data:          1 Response Variable (y)
               1 Predictor Variable (x)
               21 Observations
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         215232.624678170,
               1.00000000000000         236355.173469681,
               1.00000000000000         77934.3524331583,
               1.00000000000000         10147.5507550350,
               1.00000000000000         564.566512170752,
               1.00000000000000         11.2324854679312};
cb=cbcse[.,1];
cstdb=cbcse[.,2];
cR2= 0.957478440825662;

@Data:            y     x@
yx={          75901    0,
            -204794    1,
             204863    2,
            -204436    3,
             253665    4,
            -200894    5,
             214131    6,
            -185192    7,
             221249    8,
            -138370    9,
             315911   10,
             -27644   11,
             455253   12,
             197434   13,
             783995   14,
             608816   15,
            1370781   16,
            1303798   17,
            2205519   18,
            2408860   19,
            3444321   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   "   ) ";
"Wampler4 ";
end;
