Глава 4

Проверка гипотез

[19/95%]
Показать
LaTeX
Задача 4.1

Проводятся аэронаблюдения Y1,Y2,Y3Y_{1}, Y_{2}, Y_{3} и Y4Y_{4} углов θ1,θ2,θ3\theta_{1}, \theta_{2}, \theta_{3} и θ4\theta_{4} соответственно у четырёхугольника на земле. Если наблюдения подвержены независимым нормальным ошибкам с нулевыми средними и общей дисперсией σ2\sigma^{2}, выведите тестовую статистику для гипотезы о том, что четырёхугольник является параллелограммом, то есть θ1=θ3\theta_{1} = \theta_{3} и θ2=θ4\theta_{2} = \theta_{4}. (Адаптировано из Silvey [1970].)

?
Задача 4.2

Даны две линии регрессии

Yki=βkxi+εki(k=1,2;i=1,2,…,n), Y_{k i} = \beta _{k} x_{i}+\varepsilon _{k i} \quad (k = 1,2 ; i = 1,2, \ldots , n),

покажите, что FF-статистику для проверки H:β1=β2H: \beta_{1} = \beta_{2} можно записать в виде

F=(β^1−β^2)22S2(∑ixi2)−1 F = \frac{\left(\widehat{\beta }_{1}-\widehat{\beta }_{2}\right)^{2}}{2 S^{2}\left(\sum _{i} x_{i}^{2}\right)^{-1}}

Найдите RSS и RSSH\mathrm{RSS}_{H} и проверьте, что

RSSH−RSS=∑ixi2(β^1−β^2)22 \mathrm{RSS}_{H}-\mathrm{RSS} = \frac{\sum _{i} x_{i}^{2}\left(\widehat{\beta }_{1}-\widehat{\beta }_{2}\right)^{2}}{2}
?
Задача 4.3

Покажите, что обычную регрессионную модель полного ранга и гипотезу H:Aβ=H: \mathbf{A} \mathbf{\beta } = 0 можно преобразовать к модели Z=μ+η\mathbf{Z} = \mathbf{\mu }+\mathbf{\eta }, где μp+1=μp+2=⋯=μn=0\mu_{p+1} = \mu_{p+2} = \cdots = \mu_{n} = 0 и η∼Nn(0,σ2In)\mathbf{\eta } \sim N_{n}\left(\mathbf{0}, \sigma^{2} I_{n}\right), и к гипотезе H:μ1=μ2=⋯=μq=0H: \mu_{1} = \mu_{2} = \cdots = \mu_{q} = 0. Указание: выберите ортонормированный базис из p−qp-q векторов {αq+1,αq+2,…,αp}\left\{ \mathbf{\alpha }_{q+1}, \mathbf{\alpha }_{q+2}, \ldots , \mathbf{\alpha }_{p}\right\} для C(XA)\mathcal{C}\left(\mathbf{X}_{A}\right), где XA\mathbf{X}_{A} определено в упражнениях 4a, № 4; дополните его до ортонормированного базиса {α1,α2,…,αp}\left\{ \alpha_{1}, \alpha_{2}, \ldots , \alpha_{p}\right\} для C(X)\mathcal{C}(\mathbf{X}); а затем дополните ещё раз до ортонормированного базиса {α1,α2,…,αn}\left\{ \alpha_{1}, \alpha_{2}, \ldots , \alpha_{n}\right\} для ℜn\Re_{n}. Рассмотрите преобразование Z=TTY\mathbf{Z} = \mathbf{T}^{T} \mathbf{Y}, где T=(α1,α2,…,αn)\mathbf{T} = \left(\mathbf{\alpha }_{1}, \mathbf{\alpha }_{2}, \ldots , \mathbf{\alpha }_{n}\right) ортогональна.

?
Задача 4.4

Серия из n+1n+1 наблюдений Yi(i=1,2,…,n+1)Y_{i}(i = 1,2, \ldots , n+1) берётся из нормального распределения с неизвестной дисперсией σ2\sigma^{2}. После первых nn наблюдений возникает подозрение, что произошло резкое изменение среднего распределения. Выведите тестовую статистику для проверки гипотезы о том, что (n+1)(n+1)-е наблюдение имеет то же среднее значение генеральной совокупности, что и предыдущие наблюдения.

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

Докажите, что RSSH−RSS≥0\mathrm{RSS}_{H}-\mathrm{RSS} \geq 0.

?
Задача 4a.2

Если H:Aβ=cH: \mathbf{A} \mathbf{\beta } = \mathbf{c} верна, покажите, что FF можно выразить в виде

n−pq⋅εT(P−PH)εεT(In−P)ε \frac{n-p}{q} \cdot \frac{\varepsilon ^{T}\left(\mathbf{P}-\mathbf{P}_{H}\right) \mathbf{\varepsilon }}{\mathbf{\varepsilon }^{T}\left(\mathbf{I}_{n}-\mathbf{P}\right) \varepsilon }
?
Задача 4a.3

Если λ^H\widehat{\lambda }_{H} — оценка наименьших квадратов множителя Лагранжа, связанного с ограничениями Aβ=c\mathbf{A} \mathbf{\beta } = \mathrm{c} (ср. раздел 3.8), покажите, что

RSSH−RSS=σ2λ^HT(Var⁡[λ^H])−1λ^H \mathrm{RSS}_{H}-\mathrm{RSS} = \sigma ^{2} \widehat{\lambda }_{H}^{T}\left(\operatorname {Var}\left[\widehat{\lambda }_{H}\right]\right)^{-1} \widehat{\lambda }_{H}

(Эта идея используется для построения тестов на основе множителей Лагранжа.)

?
Задача 4a.4

Предположим, что мы хотим проверить Aβ=0\mathbf{A} \mathbf{\beta } = \mathbf{0}, где A\mathbf{A} имеет размер q×pq \times p и ранг qq. Предположим, что последние qq столбцов A\mathbf{A}, скажем A2\mathbf{A}_{2}, линейно независимы, так что A=(A1,A2)\mathbf{A} = \left(\mathbf{A}_{1}, \mathbf{A}_{2}\right), где A2\mathbf{A}_{2} — невырожденная матрица. Выразив β2\mathbf{\beta }_{2} через β1\beta_{1}, найдите матрицу XA\mathbf{X}_{A} такую, что при HH линейная модель принимает вид E[Y]=XAγ\mathbb {E}\left[\mathbf{Y}\right] = \mathbf{X}_{A} \mathbf{\gamma }. Докажите, что XA\mathbf{X}_{A} имеет полный ранг.

?
Задача 4a.5

Рассмотрим модель полного ранга с Xβ=(X1,X2)(β1T,β2T)T\mathbf{X} \mathbf{\beta } = \left(\mathbf{X}_{1}, \mathbf{X}_{2}\right)\left(\mathbf{\beta }_{1}^{T}, \mathbf{\beta }_{2}^{T}\right)^{T}, где X2\mathbf{X}_{2} имеет размер n×qn \times q.

?
(a)

Получите статистику критерия для проверки H:β2=0H: \beta_{2} = 0 в виде правой части теоремы 4.1(i). Указание: используйте A.9.1.

(b)

Найдите E[RSSH−RSS]\mathbb {E}\left[\mathrm{RSS}_{H}-\mathrm{RSS}\right].

§
Задача 4b.1

Пусть Yi=β0+β1xi1+⋯+βp−1xi,p−1+εi,i=1,2,…,nY_{i} = \beta_{0}+\beta_{1} x_{i 1}+\cdots +\beta_{p-1} x_{i, p-1}+\varepsilon_{i}, i = 1,2, \ldots , n, где εi\varepsilon_{i} независимы и распределены как N(0,σ2)N\left(0, \sigma^{2}\right). Докажите, что FF-статистика для проверки гипотезы H:βq=βq+1=⋯=βp−1=0(0<q≤p−1)H: \beta_{q} = \beta_{q+1} = \cdots = \beta_{p-1} = 0(0 < q \leq p-1) не изменяется, если из каждого YiY_{i} вычесть константу cc.

?
Задача 4b.2

Пусть Yi=β0+β1xi+εi,(i=1,2,…,n)Y_{i} = \beta_{0}+\beta_{1} x_{i}+\varepsilon_{i},(i = 1,2, \ldots , n), где εi\varepsilon_{i} независимы и распределены как N(0,σ2)N\left(0, \sigma^{2}\right).

?
(a)

Покажите, что коэффициент корреляции β^0\widehat{\beta }_{0} и β^1\widehat{\beta }_{1} равен −nxˉ/(n∑xi2)-n \bar{x} /\left(n \sqrt{\sum x_{i}^{2}}\right).

(b)

Выведите FF-статистику для проверки H:β0=0H: \beta_{0} = 0.

Задача 4b.3

Считая, что xˉ=0\bar{x} = 0, выведите FF-статистику для проверки гипотезы HH : β0=β1\beta_{0} = \beta_{1} в упражнении № 2 выше. Покажите, что она эквивалентна некоторому tt-критерию.

?
Задача 4b.4

Пусть

Y1=θ1+θ2+ε1Y2=2θ2+ε2 \begin{aligned} & Y_{1} = \theta _{1}+\theta _{2}+\varepsilon _{1} \\ & Y_{2} = 2 \theta _{2}+\varepsilon _{2} \end{aligned}

и

Y3=−θ1+θ2+ε2 Y_{3} = -\theta _{1}+\theta _{2}+\varepsilon _{2}

где εi(i=1,2,3)\varepsilon_{i}(i = 1,2,3) независимы и распределены как N(0,σ2)N\left(0, \sigma^{2}\right). Выведите FF-статистику для проверки гипотезы H:θ1=2θ2H: \theta_{1} = 2 \theta_{2}.

?
Задача 4b.5

Дано Y=θ+ε\mathbf{Y} = \theta +\varepsilon, где ε∼N4(0,σ2I4)\varepsilon \sim N_{4}\left(\mathbf{0}, \sigma^{2} \mathbf{I}_{4}\right) и θ1+θ2+θ3+θ4=0\theta_{1}+\theta_{2}+\theta_{3}+\theta_{4} = 0, покажите, что FF-статистика для проверки H:θ1=θ3H: \theta_{1} = \theta_{3} равна

2(Y1−Y3)2(Y1+Y2+Y3+Y4)2 \frac{2\left(Y_{1}-Y_{3}\right)^{2}}{\left(Y_{1}+Y_{2}+Y_{3}+Y_{4}\right)^{2}}
?
§
Задача 4c.1

Предположим, что β1=β2=⋯=βp−1=0\beta_{1} = \beta_{2} = \cdots = \beta_{p-1} = 0. Найдите распределение R2R^{2} и тем самым докажите, что

E[R2]=p−1n−1 \mathbb {E}\left[R^{2}\right] = \frac{p-1}{n-1}
?
Задача 4c.2

Для общей линейной регрессионной модели полного ранга докажите, что R2R^{2} и FF-статистика для проверки H:βj=0(j≠0)H: \beta_{j} = 0(j \neq 0) не зависят от единиц измерения, в которых измеряются YiY_{i} и xijx_{i j}.

?
Задача 4c.3

Дана модель полного ранга; предположим, что мы хотим проверить H:βj=0H: \beta_{j} = 0, j≠0j \neq 0. Пусть RH2R_{H}^{2} — коэффициент детерминации для модели с βj=0\beta_{j} = 0.

?
(a)

Докажите, что FF-статистика для проверки HH имеет вид

F=R2−RH21−R2⋅n−p1 F = \frac{R^{2}-R_{H}^{2}}{1-R^{2}} \cdot \frac{n-p}{1}

(Этот результат показывает, что FF является тестом на значимое уменьшение R2R^{2}.)

(b)

Выведите отсюда, что R2R^{2} никогда не может увеличиться, если один из коэффициентов β\beta полагается равным нулю.

§
Задача 4d.1

Выразите гипотезы из примеров 4.4 и 4.5 в каноническом виде.

?
Задача 4d.2

Предположим, что у нас есть n1n_{1} наблюдений над w1,w2,…,wp−1w_{1}, w_{2}, \ldots , w_{p-1} и UU, что даёт модель

Ui=γ0(1)+γ1(1)wi1+⋯+γp−1(1)wi,p−1+ηi(i=1,2,…,n1) U_{i} = \gamma _{0}^{(1)}+\gamma _{1}^{(1)} w_{i 1}+\cdots +\gamma _{p-1}^{(1)} w_{i, p-1}+\eta _{i} \quad \left(i = 1,2, \ldots , n_{1}\right)

Теперь нам даны ещё n2(>pn_{2}( > p) дополнительных наблюдений, которые можно выразить таким же образом, а именно

Ui=γ0(2)+γ1(2)wi1+⋯+γp−1(2)wi,p−1+ηi(i=n1+1,n2+2,…,n1+n2) \begin{aligned} U_{i} = & \gamma _{0}^{(2)}+\gamma _{1}^{(2)} w_{i 1}+\cdots +\gamma _{p-1}^{(2)} w_{i, p-1}+\eta _{i} \\ & \left(i = n_{1}+1, n_{2}+2, \ldots , n_{1}+n_{2}\right) \end{aligned}

Выведите FF-статистику для проверки гипотезы HH о том, что дополнительные наблюдения получены из той же модели.

?