Prediction from a fitted regression

[7/100%]
Show
LaTeX

Show that the square of the total bias ∥μ−E[Xβ^]∥2\left\| \mathbf{\mu }-\mathbb {E}\left[\mathbf{X} \widehat{\mathbf{\beta }}\right]\right\|^{2} is equal to the first term in (12.8).

?

Consider the regression model

Yi=α+γ1xi+γ2zi+εi(i=1,…,n)(12.13) Y_{i} = \alpha +\gamma _{1} x_{i}+\gamma _{2} z_{i}+\varepsilon _{i} \quad (i = 1, \ldots , n) \tag {12.13}

where ∑ixi=∑izi=0\sum_{i} x_{i} = \sum_{i} z_{i} = 0 and ∑ixi2=∑izi2=1\sum_{i} x_{i}^{2} = \sum_{i} z_{i}^{2} = 1. Suppose that we want to predict the response YY corresponding to a vector (x,z)(x, z) using a predictor of the form

Y^=Yˉ+γ~1x \widehat{Y} = \bar{Y}+\tilde{\gamma }_{1} x

where γ~1\tilde{\gamma }_{1} is the least squares estimate obtained by fitting the model

Yi=α+γ1xi+εi(i=1,…,n) Y_{i} = \alpha +\gamma _{1} x_{i}+\varepsilon _{i} \quad (i = 1, \ldots , n)
?
(a)

Assuming that the model (12.13) is correct, show that the expected model error E[ME]=E[(Y^−α−γ1x−γ2z)2]\mathbb {E}\left[\mathrm{ME}\right] = \mathbb {E}\left[\left(\widehat{Y}-\alpha -\gamma_{1} x-\gamma_{2} z\right)^{2}\right] of the predictor Y^\widehat{Y} is given by

E[ME]=σ2(n−1+x2)+γ22(rx−z)2, \mathbb {E}\left[\mathrm{ME}\right] = \sigma ^{2}\left(n^{-1}+x^{2}\right)+\gamma _{2}^{2}(r x-z)^{2},

where r=∑ixizir = \sum_{i} x_{i} z_{i}.

(b)

Under what circumstances will Y^\widehat{Y} be a better predictor than the predictor based on both xx and zz?

Find the form of the posterior predictive density for straight-line regression. Assume that the explanatory variable has been centered and that the prior mean is zero. Also assume that you are predicting a single value (i.e., assume that m=1m = 1).

?

In fitting the straight line Yi=β0+β1xi+εi(i=1,2,…,n)Y_{i} = \beta_{0}+\beta_{1} x_{i}+\varepsilon_{i}(i = 1,2, \ldots , n), prove that Yˉ\bar{Y} and β^1\widehat{\beta }_{1} are uncorrelated. If Y^0=β^0+β^1x0\widehat{Y}_{0} = \widehat{\beta }_{0}+\widehat{\beta }_{1} x_{0}, deduce that

Var⁡[Y^0]=σ2{1n+(x0−xˉ)2∑(xi−xˉ)2} \operatorname {Var}\left[\widehat{Y}_{0}\right] = \sigma ^{2}\left\{ \frac{1}{n}+\frac{\left(x_{0}-\bar{x}\right)^{2}}{\sum \left(x_{i}-\bar{x}\right)^{2}}\right\}
?

In Section 9.2.2 suppose that x0T=(x10T,x20T)\mathbf{x}_{0}^{T} = \left(\mathbf{x}_{10}^{T}, \mathbf{x}_{20}^{T}\right), where x0\mathbf{x}_{0} is p×1p \times 1 and x10\mathbf{x}_{10} is k×1k \times 1. The prediction at x0\mathbf{x}_{0} is Y^0=x0Tβ^\widehat{Y}_{0} = \mathbf{x}_{0}^{T} \widehat{\mathbf{\beta }}, whereas the correct prediction is Y^10=x10Tβ^1\widehat{Y}_{10} = \mathbf{x}_{10}^{T} \widehat{\mathbf{\beta }}_{1}. Compare Y^0\widehat{Y}_{0} and Y^10\widehat{Y}_{10} with respect to their means and variances.

?

Given the predictor Y^=β^0+β^1x1+⋯+β^p−1xp−1\widehat{Y} = \widehat{\beta }_{0}+\widehat{\beta }_{1} x_{1}+\cdots +\widehat{\beta }_{p-1} x_{p-1}, show that Y^\widehat{Y} has a minimum variance of σ2/n\sigma^{2} / n at the xx point xj=xˉ.j(j=1,2,…,p−1)x_{j} = \bar{x}_{. j}(j = 1,2, \ldots , p-1). Hint: Consider the model

Yi=α0+β1(xi1−xˉ.1)+⋯+βp−1(xi,p−1−xˉ.p−1)+εi. Y_{i} = \alpha _{0}+\beta _{1}\left(x_{i 1}-\bar{x}_{.1}\right)+\cdots +\beta _{p-1}\left(x_{i, p-1}-\bar{x}_{. p-1}\right)+\varepsilon _{i} .
?

Generalize the argument given in Section 5.4; namely, show that the addition of several regressors to a regression model cannot decrease the variance of the prediction Y^\widehat{Y}.

?
Remark.
?

Such a proof is, of course, not necessary, as we can add in the regressors just one at a time and evoke (5.29).

References

[SeLe03] Seber G.A., Lee A.J. Linear Regression Analysis. 2nd ed., John Wiley & Sons, 2003.