new;
output file=c:\jae\ols\wamp5\wampler5.out  reset;

@Dataset Name:  Wampler-5 (wampler5.dat)
               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                               21523262.4678170,
               1.00000000000000         23635517.3469681,
               1.00000000000000         7793435.24331583,
               1.00000000000000         1014755.07550350,
               1.00000000000000         56456.6512170752,
               1.00000000000000         1123.24854679312};
cb=cbcse[.,1];
cstdb=cbcse[.,2];
cR2= 0.224668921574940E-02;

@Data:            y     x@
yx={                      7590001     0,
           -20479994     1 ,
            20480063     2,
           -20479636     3,
            25231365     4,
           -20476094     5,
            20489331     6,
           -20460392     7,
            18417449     8,
           -20413570     9,
            20591111    10,
           -20302844    11,
            18651453    12,
           -20077766    13,
            21059195    14,
           -19666384    15,
            26348481    16,
           -18971402    17,
            22480719    18,
           -17866340    19,
            10958421    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   "   ) ";
"Wampler5";
end;

