Глава 9

Отклонения от исходных предположений

[17/100%]
Показать
LaTeX
Задача 9.1

Предположим, что постулируемая модель регрессии имеет вид

E[Y]=β0+β1x \mathbb {E}\left[Y\right] = \beta _{0}+\beta _{1} x

тогда как истинная модель на самом деле такова:

E[Y]=β0+β1x+β2x2+β3x3 \mathbb {E}\left[Y\right] = \beta _{0}+\beta _{1} x+\beta _{2} x^{2}+\beta _{3} x^{3}

Если мы используем наблюдения YY при x=−3,−2,−1,0,1,2,3x = -3,-2,-1,0,1,2,3 для оценивания β0\beta_{0} и β1\beta_{1} в постулируемой модели, какое смещение будет внесено в эти оценки?

?
Задача 9.2

Предположим, что модель

Yi=β0+β1xi1+εi Y_{i} = \beta _{0}+\beta _{1} x_{i 1}+\varepsilon _{i}

подгоняется, тогда как истинная модель на самом деле такова:

E[Yi]=β0+β1xi1+β2xi2 \mathbb {E}\left[Y_{i}\right] = \beta _{0}+\beta _{1} x_{i 1}+\beta _{2} x_{i 2}

Если eie_{i} — остаток подогнанной модели, докажите, что

E[ei]=β2(xi2+gxi1+h) \mathbb {E}\left[e_{i}\right] = \beta _{2}\left(x_{i 2}+g x_{i 1}+h\right)

где gg и hh — функции от xijx_{i j}.

?
Задача 9.3

Предположим, что мы хотим проверить гипотезу HH о том, что средние двух генеральных совокупностей равны, при наличии nin_{i} наблюдений из ii-й совокупности (i=1,2i = 1,2). Предполагая, что совокупности имеют одинаковую дисперсию и эксцесс (γ2\gamma_{2}), найдите приближённые выражения для E[Z]\mathbb {E}\left[Z\right] и Var⁡[Z]\operatorname {Var}\left[Z\right] в предположении, что HH верна (ср. теорему 9.2). Покажите, что в используемом порядке приближения эти выражения не зависят от γ2\gamma_{2}, если n1=n2n_{1} = n_{2}.

?
Задача 9.4

Докажите, что (Xi,Yi)(X_{i}, Y_{i}) в примере 9.3 раздела 9.6.5 имеет заявленное двумерное нормальное распределение.

?
§
Задача 9a.1

Предположим, что в примере 9.1 роли постулированной модели и истинной модели поменялись местами. Найдите смещения оценок наименьших квадратов.

?
Задача 9a.2

В разделе 9.2.2 наблюдаемый остаток равен Y−Xβ^\mathbf{Y}-\mathbf{X} \widehat{\mathbf{\beta }}. Найдите его истинное среднее и ковариационную матрицу.

?
Задача 9a.3

В разделе 9.2.2 предположим, что x0T=(x10T,x20T)\mathbf{x}_{0}^{T} = \left(\mathbf{x}_{10}^{T}, \mathbf{x}_{20}^{T}\right), где x0\mathbf{x}_{0} имеет размер p×1p \times 1, а x10\mathbf{x}_{10} — размер k×1k \times 1. Прогноз в точке x0\mathbf{x}_{0} равен Y^0=x0Tβ^\widehat{Y}_{0} = \mathbf{x}_{0}^{T} \widehat{\mathbf{\beta }}, тогда как правильный прогноз равен Y^10=x10Tβ^1\widehat{Y}_{10} = \mathbf{x}_{10}^{T} \widehat{\mathbf{\beta }}_{1}. Сравните Y^0\widehat{Y}_{0} и Y^10\widehat{Y}_{10} по их средним и дисперсиям.

?
§
Задача 9b.1

Если первый столбец X\mathbf{X} равен 1n\mathbf{1}_{n} и

V=(1−ρ)In+ρ1n1nT(0≤ρ<1) \mathbf{V} = (1-\rho ) \mathbf{I}_{n}+\rho \mathbf{1}_{n} \mathbf{1}_{n}^{T} \quad (0 \leq \rho < 1)

воспользуйтесь (9.8), чтобы показать, что

E[S2]=σ2(1−ρ) \mathbb {E}\left[S^{2}\right] = \sigma ^{2}(1-\rho )

Указание: рассмотрите P1n\mathbf{P 1}_{n}.

?
Задача 9b.2

Когда Var⁡[ε]=σ2V\operatorname {Var}\left[\varepsilon \right] = \sigma^{2} \mathbf{V}, подходящей оценкой β\mathbf{\beta } является обобщённая оценка наименьших квадратов β∗=(XTV−1X)−1XTV−1Y\beta * = \left(\mathbf{X}^{T} \mathbf{V}^{-1} \mathbf{X}\right)^{-1} \mathbf{X}^{T} \mathbf{V}^{-1} \mathbf{Y}. Если C(V−1X)=C(X)\mathcal{C}\left(\mathbf{V}^{-1} \mathbf{X}\right) = \mathcal{C}(\mathbf{X}), покажите, что β∗\mathbf{\beta }^{*} и β^\widehat{\mathbf{\beta }} совпадают. Указание: V−1X=XW\mathbf{V}^{-1} \mathbf{X} = \mathbf{X W}, где W\mathbf{W} невырождена.

?
§
Задача 9c.1

Заполните детали доказательства теоремы 9.2.

?
Задача 9c.2

Покажите, что теория раздела 9.5.2 может быть также применена к случаю H:Aβ=cH: \mathbf{A} \mathbf{\beta } = \mathbf{c}, где c≠0\mathrm{c} \neq \mathbf{0}.

?
Задача 9c.3

Рассмотрим регрессионную модель полного ранга

Yi=β0+β1xi1+⋯+βp−1xi,p−1+εi(i=1,2,…,n) Y_{i} = \beta _{0}+\beta _{1} x_{i 1}+\cdots +\beta _{p-1} x_{i, p-1}+\varepsilon _{i} \quad (i = 1,2, \ldots , n)

где εi\varepsilon_{i} независимо и одинаково распределены как N(0,σ2)N\left(0, \sigma^{2}\right), и предположим, что мы хотим проверить H:β1=β2=⋯=βp−1=0H: \beta_{1} = \beta_{2} = \cdots = \beta_{p-1} = 0.

Предполагая, что HH верна, найдите приближённое выражение для E[Z]\mathbb {E}\left[Z\right], где Z=12log⁡FZ = \frac{1}{2} \log F, через диагональные элементы X(XTX)−1XT\mathbf{X}\left(\mathbf{X}^{T} \mathbf{X}\right)^{-1} \mathbf{X}^{T}.

?
§
Задача 9d.1

Докажите, что в регрессии с двумя объясняющими переменными xx и zz

Var⁡[β^0]=σ2n[1+11−r2(1CVx2−2rCVxCVz+1CVz2)] \operatorname {Var}\left[\widehat{\beta }_{0}\right] = \frac{\sigma ^{2}}{n}\left[1+\frac{1}{1-r^{2}}\left(\frac{1}{\mathrm{CV}_{x}^{2}}-\frac{2 r}{\mathrm{CV}_{x} \mathrm{CV}_{z}}+\frac{1}{\mathrm{CV}_{z}^{2}}\right)\right]
?
Задача 9d.2

Докажите, что Rjy^=0R_{j}^{\widehat{y}} = 0 тогда и только тогда, когда x∗(j)\mathrm{x}^{*(j)} ортогонален столбцам X∗(j)\mathbf{X}^{*(j)}.

?
Задача 9d.3

Покажите, что Rxx\mathbf{R}_{x x} положительно определена при условии, что X\mathbf{X} имеет полный ранг и содержит столбец из единиц.

?
Задача 9d.4

Покажите, что собственные значения λj\lambda_{j} матрицы Rxx\mathbf{R}_{x x} удовлетворяют условию

1≤λj≤p−1 1 \leq \lambda _{j} \leq p-1

Отсюда докажите, что VIFj≤κ2\mathrm{VIF}_{j} \leq \kappa^{2}, где κ\kappa — число обусловленности матрицы Rxx\mathbf{R}_{x x}.

?
Задача 9d.5

Докажите, что

∂∂xi∑i=1n(xi−xˉ)(Yi−Yˉ)=Yi−Yˉ \frac{\partial }{\partial x_{i}} \sum _{i = 1}^{n}\left(x_{i}-\bar{x}\right)\left(Y_{i}-\bar{Y}\right) = Y_{i}-\bar{Y}

и

∂∂xi∑i=1n(xi−xˉ)2=2(xi−xˉ). \frac{\partial }{\partial x_{i}} \sum _{i = 1}^{n}\left(x_{i}-\bar{x}\right)^{2} = 2\left(x_{i}-\bar{x}\right) .
?