#include "libscl.h" 

using namespace scl;

using std::cout;

int main(int argc, char** argp, char** envp)
{

  realmat data;
  vecread("regr.dat",data);

  realmat y = data("",1);
  realmat X = data("",seq(2,data.ncol()));
  
  realmat b = inv(T(X)*X)*T(X)*y;

  realmat sse = T(y - X*b)*(y - X*b);

  realmat V = sse[1]*inv(T(X)*X)/(y.nrow()-X.ncol());

  cout << starbox("/Estimate of b and its variance//");

  cout << b << V;

  return 0;
}