Within this very nice piece, Rob drops this bomb of mathematical knowledge:
It is not necessary to actually fit 
separate models when computing the CV statistic for linear models.
separate models when computing the CV statistic for linear models.Say what?
Here is a broader excerpt and the method itself (after the jump).
While cross-validation can be computationally expensive in general, it is very easy and fast to compute LOOCV for linear models. A linear model can be written as
![\[<br />
\mathbf{Y} = \mathbf{X}\mbox{\boldmath$\beta$} + \mathbf{e}.<br />
\]](https://lh3.googleusercontent.com/blogger_img_proxy/AEn0k_ufX7zcpnmg8ZOgXs0OA8m0V8MXGh29WrXyU46ff3Ci3Ji3trQcCkBcJJZHEguV8j_neibTwrZ2LjBH1e65xwbL7_poAICk4cIbHWN3Z1Cb11co5i1YfXQRsMe9DDw-y9folYwMiLYCVSmV4xXXWtBU3oKQoqjO-bDw6GB9-zlqARRDOmmWa012PxqdP_xb1VM=s0-d "Rendered by QuickLaTeX.com")
Then
![\[<br />
\hat{\mbox{\boldmath$\beta$}} = (\mathbf{X}'\mathbf{X})^{-1}\mathbf{X}'\mathbf{Y}<br />
\]](https://lh3.googleusercontent.com/blogger_img_proxy/AEn0k_sR4YpQ4Ao2GTw7P3n35fBuvvdz_Sazvu5xJX6rSGXNNnf9wjcYSb8MWQ4EUJMqXcXxCw3QW0Fe8CoVVjfPX8F-lv_FLyj4jTEkNzaPE3HV3jex3JoZo8xZDqKR1-ZR9gQBUKofpTUhWxEWCM92323odLkEhzCfKqsAo5kuGUbmJDQZTQPhQEM28kYhga9A9j8I=s0-d "Rendered by QuickLaTeX.com")
and the fitted values can be calculated using
![\[<br />
\mathbf{\hat{Y}} = \mathbf{X}\hat{\mbox{\boldmath$\beta$}} = \mathbf{X}(\mathbf{X}'\mathbf{X})^{-1}\mathbf{X}'\mathbf{Y} = \mathbf{H}\mathbf{Y},<br />
\]](https://lh3.googleusercontent.com/blogger_img_proxy/AEn0k_slyWy-U4-ds8lVUoWrUjIsuE5aqNyB7HPCaB73Fc8OqjKGK_YLA6I55Vgm8T5dmER9M0cIAwh8H_4_bj_PvKMDnR9Hf9rXtNh_ZX5Z8B1eQjB2OKS4o-JAPT5MXYsnHMPtm6OjEF8M65WCcjhtpHCm_khIwzva_vwKYRf7YNZr0sE3ae_3VFjUEF3XgQn0dtw=s0-d "Rendered by QuickLaTeX.com")
where 
is known as the “hat-matrix” because it is used to compute 
(“Y-hat”).
is known as the “hat-matrix” because it is used to compute (“Y-hat”).If the diagonal values of 
are denoted by 
, then the cross-validation statistic can be computed using
are denoted by , then the cross-validation statistic can be computed using![\[<br />
\text{CV} = \frac{1}{n}\sum_{i=1}^n [e_{i}/(1-h_{i})]^2,<br />
\]](https://lh3.googleusercontent.com/blogger_img_proxy/AEn0k_sWyxB4z_EpR1UzhmtsGDXVNXZbPcUxCV4ZfbXYWFdVjLD_Sol2rj7WK7LQffScDrSXe6QU8o8QHwT1XD6lsjgZM3Ykwql-fdYt5Zh7bGeslRPoJS8x3EZcAYs26g5QOR5TQ0ClyLIX08yyKWWsGl9G5gu-4OspCGR5SZRMo3jHkT3fJaTpv5OMFSQXfG2mAVq2=s0-d "Rendered by QuickLaTeX.com")
where 
is the residual obtained from fitting the model to all observations. See Christensen’s book Plane Answers to Complex Questions for a proof. Thus, it is not necessary to actually fit separate models when computing the CV statistic for linear models. This remarkable result allows cross-validation to be used while only fitting the model once to all available observations.
is the residual obtained from fitting the model to all observations. See Christensen’s book Plane Answers to Complex Questions for a proof. Thus, it is not necessary to actually fit separate models when computing the CV statistic for linear models. This remarkable result allows cross-validation to be used while only fitting the model once to all available observations.
No comments:
Post a Comment