@---begin file named filip.txt  @
new; output file=c:\jae\ols\filip\filip.out reset;

@Make a two-column matrix of true certified b and SE(b)@
cbcse={
              -1467.48961422980         298.084530995537,
              -2772.17959193342         559.779865474950,
              -2316.37108160893         466.477572127796,
              -1127.97394098372         227.204274477751,
               -354.478233703349         71.6478660875927,
                -75.1242017393757         15.2897178747400,
                -10.8753180355343         2.23691159816033,
                -1.06221498588947         0.221624321934227,
                -0.670191154593408E-01    0.142363763154724E-01,
                -0.246781078275479E-02    0.535617408889821E-03,
                -0.402962525080404E-04    0.896632837373868E-05};


cb=cbcse[.,1]; @first column of cbcse has certified coefficients b@
cstdb=cbcse[.,2]; 
cr2= 0.996727416185620; @this is certified R-square@
format /le 23,15;

load xx[]=c:\jae\ols\filip\filip.dat;
x0=reshape(xx,82,2);
y=x0[.,1];
x1=x0[.,2];
ane=ones(82,1);
x12=(x1)^2;
x13=(x1)^3;
x14=(x1)^4;
x15=(x1)^5;
x16=(x1)^6;
x17=(x1)^7;
x18=(x1)^8;
x19=(x1)^9;
x110=(x1)^(10);

@rhs=ane~x1~x12~x13~x14~x15~x16~x17~x18~x19~x110;@
rhs=x1~x12~x13~x14~x15~x16~x17~x18~x19~x110;

ans1=nistols(y,rhs,cb,cstdb,cr2);  
ans2=niolsqr2(y,rhs,cb,cstdb,cr2);  
format /ld  6,1;
"        " ans1;  
"   (          "  ans2   "   ) ";
" Filippe  or Filip dataset from NIST-StRD";
"  ";
end;
@---end file named filip.txt  @

