Эконометрика. 6

                  Содержание: 
 
 
 
 
 
 
 
 

  1. Регрессионный анализ: понятие, задачи, основные цели……………3
 
  1. Прогнозирование, основанное на использовании  моделей временных  рядов…………………………………………………………………………5
 
 

     3.Задача№1…………………………………………………………………….8 

     4.Задача№2……………………………………………………………………25 

     5.Список используемой литературы………………………………………42 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

1.Регрессионный  анализ: понятие,  задачи, основные  цели. 

Регрессио́нный (линейный) анализ — статистический метод исследования зависимости  между зависимой переменной Y и одной или несколькими независимыми переменными X1,X2,...,Xp. Независимые переменные иначе называют регрессорами или предикторами, а зависимые переменные — критериальными. Терминология зависимых и независимых переменных отражает лишь математическую зависимость переменных (см. Ложная корреляция), а не причинно-следственные отношения. Цели регрессионного анализа

Определение степени  детерминированности вариации критериальной (зависимой) переменной предикторами (независимыми переменными)

Предсказание  значения зависимой переменной с  помощью независимой(-ых)

Определение вклада отдельных независимых переменных в вариацию зависимой

Регрессионный анализ нельзя использовать для определения  наличия связи между переменными, поскольку наличие такой связи  и есть предпосылка для применения анализа.

Регрессионным анализом называется определение аналитического выражения связи между исследуемыми переменными, в котором изменение  результативной переменной происходит под влиянием факторной переменной.

Модель регрессии  или уравнение регрессии позволяет  количественно оценить взаимосвязь  между исследуемыми переменными.

Предположим, что  имеется набор значений двух переменных: yi (результативная переменная) и xi (факторная переменная). Между этими переменными существует зависимость вида: y = f (x).

Задача регрессионного анализа состоит в том, чтобы  по данным наблюдений определить такую  функцию ỹ = f (x), которая наилучшим образом описывала исследуемую зависимость между переменными.

Для определения  аналитической формы зависимости  между исследуемыми переменными  применяются следующие методы:

1) графический  метод или визуальная оценка  характера связи. В этом случае  на линейном графике по оси  абсцисс откладываются значения  факторной переменной х, а по  оси ординат – значения результативной  переменной у. Затем на пересечении  соответствующих значений отмечаются  точки. Полученный точечный график  в системе координат (х, у)  называется корреляционным полем.  Линия, которая соединяет точки  на графике, называется эмпирической  линией. По её виду можно судить  не только о наличии, но и  о форме зависимости между  изучаемыми переменными;

2) на основе  теоретического и логического  анализа природы изучаемых явлений,  их социально-экономической сущности;

3) определение  аналитической формы зависимости  между переменными экспериментальным  путём.

При исследовании зависимости между двумя переменными  чаще всего используется линейная форма  связи. Это связано с двумя  обстоятельствами:

1) чёткая экономическая  интерпретация параметров линейной  модели регрессии;

2) в большинстве  случаев нелинейные модели регрессии  преобразуются к линейному виду.

Общий вид модели парной регрессии зависимости переменной у от переменной х: 

yi=β0+β1xi+εi,

где yi– результативные переменные,

xi– факторные переменные,

β0, β1 – параметры модели регрессии, подлежащие оцениванию;

εi – случайная ошибка модели регрессии. Данная величина является случайной, она характеризует отклонения реальных значений результативных переменных от теоретических, рассчитанных по уравнению регрессии.

Присутствие случайной  ошибки в модели регрессии порождено  следующими источниками:

1) нерепрезентативность  выборки. Модель парной регрессии  в большинстве случаев является  большим упрощением истинной  зависимости между переменными,  потому что в модель входит  только одна факторная переменная, не способная полностью объяснить  вариацию результативной переменной. При этом результативная переменная  может быть подвержена влиянию  множества других факторных переменных  в гораздо большей степени;

2) ошибки, возникающие  при измерении данных;

3) неправильная  функциональная спецификация модели.

Коэффициент β1, входящий в модельпарной регрессии, называется коэффициентом регрессии. Он характеризует, на сколько в среднем изменится результативная переменная у при условии изменения факторной переменной х на единицу своего измерения. Знак коэффициента регрессии указывает на направление связи между переменными:

1) если β1›0, то связь между изучаемыми переменными (с уменьшением факторной переменной х уменьшается и результативная переменная у, и наоборот);

2) если β1‹0, то связь между изучаемыми переменными (с увеличением факторной переменной х результативная переменная у уменьшается, и наоборот).

Коэффициент β0, входящий в модель парной регрессии, трактуется как среднее значение результативной переменной у при условии, что факторная переменная х равна нулю. Но если факторная переменная не имеет и не может иметь нулевого значения, то подобная трактовка коэффициента β0 не имеет смысла.

Общий вид модели парной регрессии в матричном  виде: 

Y= X* β+ ε, Где

– случайный  вектор-столбец значений результативной переменной размерности n x 1;

– матрица значений факторной переменной размерности  n x 2. Первый столбец является единичным, потому что в модели регрессии коэффициент β0 умножается на единицу;

– вектор-столбец  неизвестных коэффициентов модели регрессии размерности 2 x 1

– случайный  вектор-столбец ошибок модели регрессии  размерности n x 1.

2. Прогнозирование,  основанное на  использовании моделей  временных рядов.

    Характерной чертой адаптивных методов прогнозирования  является их способность непрерывно учитывать эволюцию динамических характеристик  изучаемых процессов, «адаптироваться» к этой эволюции, придавая тем больший вес, тем более высокую информационную ценность имеющимся наблюдениям, чем ближе они к текущему моменту прогнозирования.

    В основе процедуры адаптации лежит  метод проб и ошибок. По модели делается прогноз на один интервал по времени. Через один шаг моделирования  анализируется результат: насколько  он далек от фактического значения. Затем в соответствии с моделью  происходит корректировка. После этого  процесс повторяется. Таким образом, адаптация осуществляется рекуррентно  с получением каждой новой фактической  точки ряда.

    Методы  экспоненциального  сглаживания. Модель Брауна.

    Пусть анализируемый временной ряд  x(t)  представлен в виде:

x(t) = a0 + ε(t),

где  a0 – неизвестный параметр, не зависящий от времени, ε(t) - случайный остаток со средним значением, равным нулю, и конечной дисперсией.

    В соответствии с методом Брауна прогноз  x*(t+τ) для неизвестного значения  x(t+τ) по известной до момента времени t  траектории ряда  x(t)  строится по формуле:

x*(t; τ) = S(t),

где значение экспоненциально взвешенной скользящей средней  S(t)  определяется по рекуррентной формуле:

S(t)= αx(t) + (1-α) S(t-1).

    Коэффициент сглаживания α  можно интерпретировать как коэффициент дисконтирования, характеризующий меру обесценивания информации за единицу времени. Из формулы следует, что экспоненциально взвешенная скользящая средняя является взвешенной суммой всех уровней ряда  x(t), причем веса уменьшаются экспоненциально по мере удаления в прошлое.

    В качестве  S(0) берется, как правило, среднее значение ряда динамики или среднее значение нескольких начальных уровней ряда.

    Случай  линейного тренда:  x(t) = a0 + a1t + ε(t).

    В этом случае прогноз x*(t; τ) будущего значения определяется соотношением:

x*(t; τ) =

,

а пересчет коэффициентов  осуществляется по формулам:

    

    Начальные значения коэффициентов берутся  из оценки тренда линейной функцией.

    Модель  Хольта.

    В модели Хольта введено два параметра  сглаживания  α1  и α 2 (0< α 1, α 2 <1). Прогноз x*(t;l)  на  l  шагов по времени определяется формулой:

x*(t; τ) =

,

а пересчет коэффициентов  осуществляется по формулам:

    Модель  Хольта-Уинтерса.

    Эта модель помимо линейного тренда учитывает  и сезонную составляющую. Прогноз  x*(t;τ)  на  τ  шагов по времени определяется формулой:

x*(t;τ) =

,

где f(t) – коэффициент сезонности, а T – число временных тактов (фаз), содержащихся в полном сезонном цикле.

    Видно, что в данной модели сезонность представлена мультипликативно. Формулы обновления коэффициентов имеют вид:

    Модель  Тейла-Вейджа.

    Если  исследуемый временной ряд имеет  экспоненциальную тенденцию с мульти-пликативной  сезонностью, то после логарифмирования обеих частей уравнения получается модель с линейной тенденцией и аддитивной сезонностью или модель Тейла-Вейджа.

    Имеется модель:

x(t) = a0(t) + g(t) + δ(t),

a0(t) = a0(t-1) + a1(t).

Здесь a0(t) – уровень процесса после устранения сезонных колебаний, a1(t) – аддитивный коэффициент роста, ω(t) – аддитивный коэффициент сезонности и δ(t) – белый шум.

    Прогноз x*(t;τ)  на  τ  шагов по времени определяется формулой:

x*(t;τ) =

.

    Коэффициенты  вычисляются рекуррентным способом по формулам:

    Для определения оптимальных значений параметров адаптации перебирают различные  наборы их значений и сравнивают получающиеся при этом среднеквадратические ошибки прогнозов. 
 
 

 

Задача 1.

В табл. 1 приведены данные о личных потребительских расходах и располагаемом личном доходе населения США (млрд. дол., в ценах 1972 г.) за период с 1959 г. по 1983 г.

а) Постройте уравнение для функции спроса (Y) на данный товар или вид услуг в зависимости от располагаемого личного дохода (X). Постройте таблицу дисперсионного анализа. Вычислите коэффициент детерминации R2.

Проверьте регрессию на значимость с помощью  F-теста (α = 0,05 - критерий значимости). Вычислите 95% доверительные интервалы для истинных значений коэффициентов в этом уравнении. Изобразите диаграмму рассеяния и прямую регрессии.

б) Постройте линейный временной тренд для функции спроса.

в) Рассмотрите функцию спроса (Y) как функцию двух переменных: располагаемого дохода (X) и реальной цены на товар или вид услуг (P). (Реальная цена вычисляется по формуле  
P = (Z / L) ∙ 100%). (Z, L см. табл. 1).

Постройте уравнение множественной регрессии  Y на X и P. С помощью МНК оцените коэффициенты в этом уравнении (замените Y, X и P на их логарифмы).

Дайте экономическую интерпретацию коэффициентов  в уравнении. Вычислите коэффициенты эластичности функции спроса. Интерпретируйте  эти коэффициенты.

Таблица 1

Исходные  данные о личных потребительских  расходах и располагаемом личном доходе

Год Личный

располагаемый

доход,

млрд. $

Текущие

расходы

на бензин

Дефляторы цен 

для личных

потребительских

расходов 

Цены
T X Y Z L
1959 440,4 13,7 82,2 70,6
1960 452,0 14,2 84,5 71,9
1961 461,4 14,3 83,9 72,6
1962 482,0 14,9 84,5 73,7
1963 500,5 15,3 84,5 74,8
1964 528,0 16,0 84,4 75,9
1965 557,5 16,8 87,5 77,2
1966 646,8 17,8 89,5 79,4
1967 673,5 18,4 92,4 81,4
1968 701,3 19,9 93,8 84,6
1969 722,5 21,4 97,0 88,4
1970 751,6 22,9 97,9 92,5
1971 779,2 24,2 98,7 96,3
1972 810,3 25,4 109,0 100,0
1973 865,3 26,2 109,4 105,7
1974 858,4 24,8 147,7 116,3
1975 875,8 25,6 157,7 125,2
1976 906,8 26,8 164,3 131,7
1977 942,9 27,7 173,7 139,3
1978 988,8 28,3 181,3 149,1
1979 1015,5 27,4 243,2 162,5
1980 1021,6 25,1 337,9 179,0
1981 1049,3 25,1 376,4 194,3
1982 1058,3 25,3 356,6 206,0
1983 1095,4 26,1 344,9 213,6
 

Решение

а) Постройте уравнение для функции спроса (Y) на данный товар или вид услуг в зависимости от располагаемого личного дохода (X). Постройте таблицу дисперсионного анализа. Вычислите коэффициент детерминации R2.

Проверьте регрессию на значимость с помощью  F-теста (α = 0,05 - критерий значимости). Вычислите 95% доверительные интервалы для истинных значений коэффициентов в этом уравнении. Изобразите диаграмму рассеяния и прямую регрессии. 

Идентифицируем  переменные: x - независимая переменная (фактор) y - зависимая переменная (показатель).

Пусть эконометрическая модель специфицирована в линейной форме:

y = ax + b + u

где a, b - параметры модели u - стохастическая составляющая (остатки).

Используем метод  наименьших квадратов [ЛЕЩ, с.29]. 
Запишем систему нормальных уравнений, используя в качестве неизвестной переменной - переменную x:

где n - количество наблюдений.

Построим вспомогательную  таблицу 1.

Табл. 1.

x y x2 xy
1 440,4 13,7 193952,16 6033,48
2 452,0 14,2 204304,00 6418,40
3 461,4 14,3 212889,96 6598,02
4 482,0 14,9 232324,00 7181,80
5 500,5 15,3 250500,25 7657,65
6 528,0 16,0 278784,00 8448,00
7 557,5 16,8 310806,25 9366,00
8 646,8 17,8 418350,24 11513,04
9 673,5 18,4 453602,25 12392,40
10 701,3 19,9 491821,69 13955,87
11 722,5 21,4 522006,25 15461,50
12 751,6 22,9 564902,56 17211,64
13 779,2 24,2 607152,64 18856,64
14 810,3 25,4 656586,09 20581,62
15 865,3 26,2 748744,09 22670,86
16 858,4 24,8 736850,56 21288,32
17 875,8 25,6 767025,64 22420,48
18 906,8 26,8 822286,24 24302,24
19 942,9 27,7 889060,41 26118,33
20 988,8 28,3 977725,44 27983,04
21 1015,5 27,4 1031240,25 27824,70
22 1021,6 25,1 1043666,56 25642,16
23 1049,3 25,1 1101030,49 26337,43
24 1058,3 25,3 1119998,89 26774,99
25 1095,4 26,1 1199901,16 28589,94
Σ 19185,1 543,6 15835512,07 441628,55

 
Получим систему уравнений:

Решение системы  найдем по формулам Крамера [ГЕТ, с.30]:

где Δ - главный определитель системы.

Следовательно, уравнение линейной модели имеет  вид:

y = 4,8705 + 0,0220x.

Это значит, что  при увеличении или уменьшении значения фактора на 1 у.е. показатель увеличивается или уменьшается на 0,022 у.е., то есть между эконометрическими параметрами существует прямая пропорциональная или положительная зависимость.

Свободный член регрессии b = 4,8705 указывает значение показателя при нулевом значении фактора. Он имеет лишь расчетное значение, поскольку такой случай невозможен в реальной экономической ситуации.  

Для проведения исследования модели построим вспомогательную  таблицу 2.

Табл. 2.

x y yx
u2 = (y - yx)2
440,4 13,7 14,5539 64,7059 106931,6160 2630,4202 51,6973 0,7292
452,0 14,2 14,8090 56,9119 99479,6832 2379,4078 48,0945 0,3709
461,4 14,3 15,0157 55,4131 93638,4480 2277,8938 45,2705 0,5122
482,0 14,9 15,4686 46,8403 81455,4432 1953,3050 39,3805 0,3233
500,5 15,3 15,8754 41,5251 71237,7452 1719,9294 34,4407 0,3311
528,0 16,0 16,4800 32,9935 57314,2752 1375,1366 27,7092 0,2304
557,5 16,8 17,1287 24,4431 44059,6892 1037,7654 21,3011 0,1080
646,8 17,8 19,0922 15,5551 14545,3248 475,6622 7,0321 1,6698
673,5 18,4 19,6793 11,1823 8817,9612 314,0150 4,2631 1,6365
10  701,3 19,9 20,2905 3,4003 4369,7388 121,8958 2,1126 0,1525
11  722,5 21,4 20,7567 0,1183 2016,3692 15,4470 0,9748 0,4139
12  751,6 22,9 21,3965 1,3363 249,7664 -18,2694 0,1208 2,2605
13  779,2 24,2 22,0034 6,0319 139,1456 28,9710 0,0673 4,8252
14  810,3 25,4 22,6872 13,3663 1840,0668 156,8278 0,8896 7,3594
15  865,3 26,2 23,8965 19,8559 9583,6268 436,2246 4,6333 5,3061
16  858,4 24,8 23,7448 9,3391 8280,2720 278,0838 4,0032 1,1135
17  875,8 25,6 24,1274 14,8687 11749,6928 417,9750 5,6805 2,1686
18  906,8 26,8 24,8090 25,5631 19431,2448 704,7862 9,3942 3,9641
19  942,9 27,7 25,6028 35,4739 30798,8460 1045,2542 14,8900 4,3984
20  988,8 28,3 26,6120 42,9811 49016,1888 1451,4722 23,6974 2,8493
21  1015,5 27,4 27,1991 31,9903 61551,6252 1403,2310 29,7578 0,0404
22  1021,6 25,1 27,3332 11,2627 64615,6064 853,0818 31,2391 4,9872
23  1049,3 25,1 27,9423 11,2627 79465,3548 946,0430 38,4184 8,0784
24  1058,3 25,3 28,1401 12,6451 84620,4828 1034,4262 40,9107 8,0664
25  1095,4 26,1 28,9559 18,9747 107581,3760 1428,7506 52,0114 8,1561
Σ 19185,1 543,6 - 608,0416 1112789,5896 24467,7356 537,9904 70,0512
Σ/n 767,404 21,744 - 24,3217 44511,5836 978,7094 21,5196 2,8020

 
Оценим параметры модели альтернативным способом:

Линейное уравнение  регрессии аналогично: yx = 0,022x + 4,8705.  

Вычислим для  зависимой переменной y общую дисперсию, дисперсию, что объясняет регрессию, дисперсию ошибок [ЛУК, с. 56]:

 

Суммы квадратов  связанные с определенным источником вариации, а также со степенями  свободы и средними квадратами. Сведем их всех в таблице, которая называется базовой таблицей дисперсионного анализа - ANOVA-таблицей [ЛУК, с. 61].

Построим ANOVA-таблицу о зависимости между показателем и фактором:

Источник 
вариации
Количество 
степеней 
свободы
Сумма квадратов Средние квадраты
Предопределено  регрессией (модель) 1     
Необъяснимо с помощью 
регрессии (ошибки)
n - 2 = 23      
Общее n - 1 = 24 -
 

Определим коэффициенты детерминации R2 и корреляции r [ЛУК, c. 57]:

Этот результат  значит, что 88,48% вариации результативного  признака зависит от вариации уровня факторного признака, а 11,52% приходится на другие факторы.

Поскольку 0,7 < r < 1, то между факторным и результативным признаком корреляционная связь сильная.  

Проверим адекватность модели по критерию Фишера или F-критерию, 
который вычисляется по формуле:

Поскольку табличное  значение F(0,05; 1; 23) = 4,28 и |F| > Fтаб, то делаем вывод об адекватности эконометрической модели.  

Оценим статистическую значимость параметров регрессии.

Табличное значение t-критерия Стьюдента при заданном уровне значимости 0,95 и n - 2 = 23 степенях свободы равно 2,07.

Найдем матрицу  погрешностей C-1, обратную к матрице системы уравнений:

Δ = |C| = 27819739,74,

Определим стандартные  погрешности оценок параметров модели, учитывая дисперсию остатков:

где

Рассчитаем  t-критерий Стьюдента для каждого из коэффициентов

Поскольку tm = 2,07 и это значение больше t-критериев для каждого из коэффициентов, то делаем вывод об их статистической значимости.  

Определим интервалы  доверия для параметров регрессии [ЕЛИ, с. 57].

Для расчета  доверительного интервала определяем предельную погрешность для каждого  коэффициента:

Δa = tтаб.ma = 2,0690 · 0,0017 = 0,0034,

Δb = tтаб.mb = 2,0690 · 1,3167 = 2,7242,

Следовательно, экстремальные значения для каждого  коэффициента следующие:

min a = a - Δa = 0,0220 - 0,0034 = 0,0186, max a = a + Δa = 0,0220 + 0,0034 = 0,0254,  
min b = b - Δb = 4,8705 - 2,7242 = 2,1463, max b = b + Δb = 4,8705 + 2,7242 = 7,5947.

Таким образом, доверительные интервалы для  коэффициентов регрессии следующие:

a Î (0,0186; 0,0254),

b Î (2,1463; 7,5947).   
Изобразим диаграмму рассеяния и прямую регрессии.

 
 

б) Постройте линейный временной тренд для функции спроса.

 

Идентифицируем  переменные: t - независимая временная переменная (фактор), y - зависимая переменная (показатель). Пусть модель специфицирована в линейной форме:

y = at + b + u,

где a, b - параметры модели, u - стохастическая составляющая (остатки).

Используем метод  наименьших квадратов [ЛЕЩ, с.29]. 
Запишем систему нормальных уравнений, используя в качестве неизвестную переменную - переменную t:

Эконометрика. 6