Fixes #4
When doing a cholesky decomposition, the intercept is no more an intersept !
beta_0 is multiplied by L (if (LL'=Sigma).
Take this into account for the lm part of the gauss lasso, using lm.fit (instead of lm)