Func void main()

{
   Integer n,i,k;
   Real    ESS,TSS,R2,s2,modR2;
   Matrix  t,x,y,b,e,s2b,sb,tval;


   n=10;
   t=Matrix([1:n]');
   x=[ ONE(n,1) t];
   k=Cols(x);
   y=[4 5 6 7 2 3 3 1 8 5]';


   b=inv(x'*x)*x'*y;
      print b;


   e=y-x*b;
   ESS=[e'*e](1,1);
   TSS=[(y.-mean(y))'*(y.-mean(y))](1,1);

   R2=1-ESS/TSS;
      print R2;

   modR2=1-(ESS/(n-k))/(TSS/(n-1));
      print modR2;

   s2=ESS/(n-k);
   s2b=diag2vec(inv(x'*x))*s2;
   sb=s2b.^0.5;
      print sb;

   tval=b./s2b;
      print tval;

}



main()