Skip to main content
3 of 5
added 11 characters in body
Heisenberg
  • 4.6k
  • 4
  • 31
  • 65

Why does ridge estimate become better than OLS by adding a constant to the diagonal?

I understand that the ridge regression estimate is the $\beta$ that minimizes residual sum of square and a penalty on the size of $\beta$

$$\beta_{ridge} = (\lambda I_D + X'X)^{-1}X'y = argmin RSS + \lambda ||\beta||^2_2$$

However, I don't fully understand the significance of the fact that $\beta_{ridge}$ differs from $\beta_{OLS}$ by only adding a small constant to the diagonal of $X'X$. Indeed,

$$\beta_{OLS} = (X'X)^{-1}X'y$$

  1. My book mentions that this makes the estimate more stable numerically -- why?

  2. Is numerical stability related to the shrinkage towards 0 of the ridge estimate, or it's just a coincidence?

Heisenberg
  • 4.6k
  • 4
  • 31
  • 65