11

Обобщенная линейная модель регрессии

[9/0%]
LaTeX
Задача 11.1

Найдите ОМНК-оценку и дисперсию ее компонент в случае, когда ковариационная матрица ошибок VnV_{n} является диагональной.

?
Задача 11.2

Получите матрицу вида (11.12).

?
Задача 11.3

Пусть Xk=θ1+θ2hk+εk,k=1,2,3X_{k} = \theta_{1}+\theta_{2} h_{k}+\varepsilon_{k}, k = 1,2,3, где вектор ошибок E={ε1,ε2,ε3}TE = \left\{ \varepsilon_{1}, \varepsilon_{2}, \varepsilon_{3}\right\}^{T} имеет математическое ожидание E[E]=0\mathbb {E}\left[E\right] = 0, а Cov⁡[E,E]=σ2V\operatorname {Cov}\left[ E, E \right] = \sigma^{2} V,

V=[1ρaρρaa2ρaρρa1], V =\left[\begin{smallmatrix} 1 & \rho a & \rho \\ \rho a & a^{2} & \rho a \\ \rho & \rho a & 1 \end{smallmatrix}\right],

числа aa и 0<ρ<10 < \rho < 1 неизвестны, h1=−1,h2=0,h3=1h_{1} = -1, h_{2} = 0, h_{3} = 1. Покажите, что ОМНК-оценка вектора θ\theta имеет вид

θ^=[1r((a2−aρ)X1+(1−2aρ+ρ)X2+(a2−aρ)X3)−0,5X1+0,5X3], \widehat{\theta } =\left[\begin{smallmatrix} \frac{1}{r}\left(\left(a^{2}-a \rho \right) X_{1}+(1-2 a \rho +\rho ) X_{2}+\left(a^{2}-a \rho \right) X_{3}\right) \\ -0,5 X_{1}+0,5 X_{3} \end{smallmatrix}\right],

где r=1+ρ+2a2−4aρr = 1+\rho +2 a^{2}-4 a \rho. Найдите такие значения ρ,a,Xk,k=1,2,3\rho , a, X_{k}, k = 1,2,3, при которых линия регрессии лежит целиком ниже и целиком выше всех наблюдавшихся значений.

?
Задача 11.4

Пусть Xi=θhi+εi,i=1,2X_{i} = \theta h_{i}+\varepsilon_{i}, i = 1,2, где ε1∼N(0,σ2),ε2∼N(0,4σ2)\varepsilon_{1} \sim \mathcal{N}\left(0, \sigma^{2}\right), \varepsilon_{2} \sim \mathcal{N}\left(0,4 \sigma^{2}\right), причем ε1\varepsilon_{1} и ε2\varepsilon_{2} независимы. Для двух векторов h:h(1)=[−1,1]⊤,h(2)=[1,−1]⊤h: h^{(1)} = [-1,1]^{\top }, h^{(2)} = [1,-1]^{\top } найдите ОМНК-оценку и ее дисперсию для двух векторов hh.

?
Задача 11.5

Пусть модель регрессии имеет вид

Xk=θ1+θ2hk+εk,k=1,…5, X_{k} = \theta _{1}+\theta _{2} h_{k}+\varepsilon _{k}, \quad k = 1, \ldots 5,

где εk\varepsilon_{k} — независимые, центрированные СВ с дисперсиями σk2=hk2\sigma_{k}^{2} = h_{k}^{2}. По наблюдениям, приведенным в табл. 11.5 найдите реализации МНК-оценки и ОМНК-оценки неизвестных параметров [θ1,θ2]⊤\left[\theta_{1}, \theta_{2}\right]^{\top }. Для каждой оценки вычислите ковариационную матрицу ее ошибки.

Таблица 11.5

hkh_{k}12345
XkX_{k}12234
?
Задача 11.6

Покажите, что использование ОМНК для модели простой линейной регрессии с дисперсиями ошибок, равными Var⁡[ek]=σ2hk2\operatorname {Var}\left[e_{k}\right] = \sigma^{2} h_{k}^{2}, равносильно использованию обычного МНК для модели наблюдений вида

Xkhk=θ1hk+θ2+δk,k=1,…,n \frac{X_{k}}{h_{k}} = \frac{\theta _{1}}{h_{k}}+\theta _{2}+\delta _{k}, \quad k = 1, \ldots , n

где δk\delta_{k} — независимые СВ с дисперсиями Var⁡[δk]=σ2\operatorname {Var}\left[\delta_{k}\right] = \sigma^{2} для k=1,…,nk = 1, \ldots , n.

?
Задача 11.7

Пусть задана модель регрессии

Xk=θ1+θ2hk+εk,k=1,…10 X_{k} = \theta _{1}+\theta _{2} h_{k}+\varepsilon _{k}, \quad k = 1, \ldots 10

где ошибки εk\varepsilon_{k} удовлетворяют уравнению авторегрессии первого порядка

εk=−0,4εk−1+vk,k=2,…,10 \varepsilon _{k} = -0,4 \varepsilon _{k-1}+v_{k}, \quad k = 2, \ldots , 10

где vkv_{k} — последовательность независимых нормально распределенных случайных величин с нулевым средним и дисперсией σ02=1\sigma_{0}^{2} = 1. Найдите реализацию ОМНК-оценки вектора неизвестных параметров θ=[θ1,θ2]⊤\theta = \left[\theta_{1}, \theta_{2}\right]^{\top } по наблюдениям, приведенным в табл. 11.6.

Таблица 11.6

hkh_{k}5,02,51,86,89,03,86,59,01,03,5
XkX_{k}5,04,83,18,28,65,56,511,12,14,5

Сравните полученный вектор с реализацией МНК-оценки, построенной в предположении некоррелированности ошибок εk\varepsilon_{k}.

?
Задача 11.8

Рассмотрите модель, связывающую количество вакансий ww и уровень безработицы uu

Xk=θ1+θ2hk+εk,k=1,…,24, X_{k} = \theta _{1}+\theta _{2} h_{k}+\varepsilon _{k}, \quad k = 1, \ldots , 24,

где Xk=ln⁡wk,hk=ln⁡uk,εk∼N(0;σ2)X_{k} = \ln w_{k}, h_{k} = \ln u_{k}, \varepsilon_{k} \sim \mathcal{N}\left(0 ; \sigma^{2}\right). Используя данные из табл. 11.7, найдите МНК-оценки параметров θ1\theta_{1} и θ2\theta_{2}. Проверьте с помощью критерия Дарбина-Уотсона гипотезу H0:ρ=0H_{0}: \rho = 0 о некоррелированности ошибок наблюдений на уровне значимости α=0,05\alpha = 0,05. Если гипотеза H0H_{0} отвергается, то предположите, что последовательность случайных ошибок {εk}\left\{ \varepsilon_{k}\right\} удовлетворяет уравнению авторегрессии первого порядка

εk=ρεk−1+ek,k=2,…,n \varepsilon _{k} = \rho \varepsilon _{k-1}+e_{k}, \quad k = 2, \ldots , n

где eke_{k} — центрированные гауссовские СВ\mathrm{СВ} с дисперсией Var⁡[ek]=σ02\operatorname {Var}\left[e_{k}\right] = \sigma_{0}^{2}. Найдите оценку неизвестного коэффициента ρ\rho в уравнении авторегрессии.

?
Задача 11.9

Исследуйте потребительскую функцию Японии X(h)X(h) за период 1951-1962 гг. на наличие коррелированности ошибок наблюдений. Если ошибки в наблюдениях являются коррелированными, то найдите реализацию ОМНК-оценки вектора неизвестных параметров. Найдите оценку ковариационной матрицы вектора ошибок EE. Данные для расчетов приведены в табл. 11.8 (данные взяты из: U. N. Yearbook of National Accounts Statistics, 1963). Здесь XX — совокупное индивидуальное потребление (трлн иен), hh совокупные личные доходы (трлн иен).

kkXkX_{k}hkh_{k}kkXkX_{k}hkh_{k}
12,863,5275,897,00
23,514,2686,207,32
34,224,7796,708,24
44,675,25107,519,46
54,975,73118,5811,10
65,436,35129,9612,69
?