new;
output file=c:\jae\ols\wamp1\wampler1.out  reset;
@
Dataset Name:  Wampler-1 (wampler1.dat)
Data:          1 Response Variable (y)
               1 Predictor Variable (x)
               21 Observations
               Higher Level of Difficulty
               Generated Data

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,
                1.00000000000000        0.000000000000000,
                1.00000000000000        0.000000000000000,
                1.00000000000000        0.000000000000000,
                1.00000000000000        0.000000000000000,
                1.00000000000000        0.000000000000000};

cb=cbcse[.,1];
cstdb=cbcse[.,2];
cR2= 1.00000000000000;

@Data:            y     x@
yx={             1     0,
                 6     1,
                63     2,
               364     3,
              1365     4,
              3906     5,
              9331     6,
             19608     7,
             37449     8,
             66430     9,
            111111    10,
            177156    11,
            271453    12,
            402234    13,
            579195    14,
            813616    15,
           1118481    16,
           1508598    17,
           2000719    18,
           2613660    19,
           3368421    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   "   ) ";
"Wampler1";
end;

