- •Условные обозначения, принятые в работе.
- •Разделение переменных на эндогенные и экзогенные.
- •Описание показателей и формулирование рабочих гипотез.
- •Структурная форма модели в общем виде.
- •Анализ развития экономики страны в ретроспективном периоде.
- •Корреляционный анализ переменных модели.
- •Двухшаговый метод наименьших квадратов.
- •Прогноз по 1мнк.
- •Прогноз по 2мнк.
Двухшаговый метод наименьших квадратов.
Записанные выше уравнения, полученные по МНК, представляют собой структурную форму модели. Для получения оценок парметров модели по 2МНК на первом шаге необходимо записать, и восстановить параметры, приведенной формы модели.
Шаг 1. Приведенная форма модели.
Yt = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Ct = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Lft = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
It = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Wgt = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Ext = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Imt = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Unt = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Pit = F(Gt, Tt, Irt, Mt, Mt-1, Pit-1, Yt-1, Kt-1, It-1, Nt, t);
Yt = 599,324*Irt + 3,51875*Mt - 2,20675*Mt-1 + 0,680562*Yt-1 - 7,64033*t;
Ct=– 1832,52*Tt + 478,055*Irt + 2,7986*Mt - 1,8586*Mt-1 + 0,574976*Yt-1 - 4,48*t;
Lft = -0,384668*Mt-1 + 2,53119*Pit-1 + 0,727815*Nt;
It = 2,62736*Gt + 20,8567*Pit-1 - 0,352074*It-1 - 7,53186*Nt
Wg = 0,0259616*G - 0,0506729*M + 0,0139026*Y1 + 0,186522*t;
Ex = 232,531*Ir – 1,37859*M1 + 10,8018*Pi1 - 5,71101*t;
Im = 106,071*Ir + 0,512589*M + 0,776596*M1 - 5,70969*Pit-1
Un = 0,0107602*G - 0,0223452*M + 0,253997*Pit-1 - 0,00188577*Yt-1;
Pit = 7,56143*Irt + 0,120031*Mt - 0,0853*Mt-1 + 0,852586*Pit-1 - 0,0032*Yt-1;
Шаг 2.
Ln(Yt) = 0,825741*Ln(Kt) + 0,402234*Ln(Ett) - 0,0204713*t;
Ct = 0,815648*Nit - 0,815648*Ir;
Lft = 0,767357*Wgt + 1,30571*Pit;
It = 1819,06-1819,06*Ir;
Wgt = 27,7255 – 0,0514576*Lft - 0,339164*(Pit – Pit-1);
Ext = 0,0896857*Yt
Imt = 0,110128*Ct
Unt = – 0,00247893*Yt-1 + 0,125305*Lft;
Pit = – 0,0193702*Yt-1 + 0,815745*Wgt + 0,166507*M + 0,065405*Mt-1;
-
Модель:
Характеристика
S2
R2
F-stat.
Y
0,0140943
99,9997
2340645,55
C
10,2215
99,9956
524862,16
Lf
2,94556
99,9039
11952,08
I
133,831
98,0123
1183,41
Wg
0,39507
71,1915
29,42
Ex
4,1922
99,95
46005,80
Im
4,31986
99,9355
35640,82
Un
0,286114
99,8263
6321,95
Pi
0,237653
99,9984
349235,07
Итого модель мо 2МНК:
Статистические зависимости:
Ln(Yt) = 0,825741*Ln(Kt) + 0,402234*Ln(Ett) - 0,0204713*t
Yt = Kt0,825741Ett0,402234 e- 0,0204713*t;
Ct = 0,815648*Nit - 0,815648*Ir;
Lft = 0,767357*Wgt + 1,30571*Pit;
It = 1819,06-1819,06*Ir;
Wgt = 27,7255 – 0,0514576*Lft - 0,339164*(Pit – Pit-1);
Ext = 0,0896857*Yt
Imt = 0,110128*Ct
Unt = – 0,00247893*Yt-1 + 0,125305*Lft;
Pit = – 0,0193702*Yt-1 + 0,815745*Wgt + 0,166507*M + 0,065405*Mt-1;
Gt = 340,54;
Mt = 1/(0,00355437 - 0,0000436785*t);
Tt = 0,094212;
Nt = 1/(0,0099676 - 0,0000550187*t);
Irt = 0,494;
Балансовы уравнения:
Nit = Yt – Taxt;
Taxt = Tt*Yt;
Ut = (Unt/Lft)*100%;
Emt = Lft – Unt;
Et = Ext – Imt;
Xt = Ct + Gt + It + Et;
![]()
Сравнение моделей:
-
Модель:
Характеристика
S2
R2
F-stat.
2МНК: Y
0,0140943
99,9997
2340645,55
МНК: Y
0,00735974
99,9999
8942829,58
2МНК: C
10,2215
99,9956
524862,16
МНК: C
6,40144
99,9983
1393957,68
2МНК: Lf
2,94556
99,9039
11952,08
МНК: Lf
3,01168
99,898
11751,92
2МНК: I
133,831
98,0123
1183,41
МНК: I
133,831
98,0123
1183,41
2МНК: Wg
0,39507
71,1915
29,42
МНК: Wg
0,325805
80,4076
48,20
2МНК: Ex
4,1922
99,95
46005,80
МНК: Ex
3,9698
99,9552
53538,94
2МНК: Im
4,31986
99,9355
35640,82
МНК: Im
4,57495
99,9278
33225,56
2МНК: Un
0,286114
99,8263
6321,95
МНК: Un
0,426174
99,6052
2901,64
2МНК: Pi
0,237653
99,9984
349235,07
МНК: Pi
0,282685
99,9977
253156,86
Ни один метод не имеет явного приимущества перед другим. Из 9-и моделей 4-е 2МНК-модели являются лучшими и столько же 1МНК-моделей. Модель инвестиций является функцией экзогенной переменной – ставки процента, и поэтому одинакова в обоих моделях. Для прогноза будем всё же использовать 1МНК, так как модель Wg построенная по 1МНК гораздо лучше чем та же модель по 2МНК, при близости R2 в других моделях.
Проверка точности прогнозирования выбранных классов моделей на основе коэффициента тейла.
Модель 1МНК построеннтая по усеченной выборке в 22 значения периода 1976-1997гг.:
Yt = K0,801118Et0,444812 e- 0,0217497*t;
Ct = 0,814671*Nit + 0,814671*Irt;
Lft = 0,846041*Wgt + 1,2709*Pit;
It = 1775,2 – 1775,2*Irt;
Wgt = 28,3437 - 0,0592236*Lft + 0,274963*(Pit – Pit-1);
Ext = 0,0897539*Yt;
Imt = 0,110168*Ct
Unt = -0,00217527*Yt-1 + 0,118465*Lft
Pit = -0,0130677*Yt-1 + 0,292003*Wgt + 0,162779*Mt + 0,0660956*Mt-1;
Gt = 340,54;
Tt = 0,094212;
Irt = 0,494;
Mt = 1/(0,00354946 - 0,0000430572*t)
Nt = 1/(0,00993241 - 0,000050896*t)
Представим систему в виде, удобном последовательного расчета показателей:
(Здесь мы решим систему структурных уравнений относительно экзогенных и лаговых переменных. Записав первые 3-и уравнения мы можем в указанном ниже порядке получить такое выражение.)
Pit = 8.3449- 0.08095*Pit-1 - 0.013836*Yt-1+ 0.1723495*Mt+ 0.069982*Mt-1;
Wgt = 28.57822 - 0.2772381*Pit-1 - 0.0026312*Yt-1 + 0.03278*M+ 0.01331*Mt-1;
Lft = 34.78391 - 0.33744*Pit-1 - 0.01981*Y1+ 0.24677*Mt + 0.1002*Mt-1;
Unt = -0,00217527*Yt-1 + 0,118465*Lft
Ett = Lft – Unt;
Yt = K0,801118Et0,444812 e- 0,0217497*t;
It = 1775,2 – 1775,2*Irt;
Taxt = Tt*Yt;
Nit = Yt – Taxt;
Ct = 0,814671*Nit + 0,814671*Irt;
Ext = 0,0897539*Yt;
Imt = 0,110168*Ct
Et = Ext – Imt;
Ut = (Unt/Lft)*100%;
Xt = Ct + Gt + It + Et;
Значения коэффициентров Тэйла
|
Показатель |
Значение Кт. |
|
Y |
0,06072 |
|
C |
0,06655 |
|
Lf |
0,01754 |
|
Wg |
0,05115 |
|
Ex |
0,05880 |
|
Im |
0,07824 |
|
Un |
0,04435 |
|
Pi |
0,01793 |
Модель 2МНК построеннтая по усеченной выборке в 22 значения периода 1976-1997гг.:
Шаг 1. Прогнозирование эндогенных переменных.
Yt* = 593,515*Irt + 3,66584*Mt - 2,00348*Mt-1 + 0,635715*Yt-1 - 9,34186*t;
Ct*=-1795,37*Tt+472,93*Irt+2,90444*Mt-1,6959*Mt-1+0,539041*Yt-1-5,778*t;
Lft* = 1,00676*Pit-1 + 0,325944*Nt;
It* = 2,74006*Mt;
Wgt* = 0,0310719*Gt - 0,0642466*Mt + 0,0149036*Yt-1 + 0,246186*t;
Ext* = 217,679*Irt - 1,1833*Mt-1 + 9,70697*Pit-1 - 5,46336*t
Imt* = 0,547591*Gt + 2,59127*Mt-1 - 17,2266*Pit-1 + 5,85674*t
Unt* =0,019099*Gt – 0,0295377*Mt+0,053839*Mt-1-0,0064*Yt-1+0,0466258*Nt
Pit* = 8,06697*Irt + 0,116295*Mt - 0,0894732*Mt-1 + 0,885228*Pit-1 - 0,00292363*Yt-1
Шаг 2. Восстановление коэффициентов структурной формы.
Yt = K1,28307Et-0,379686 e- 0,015362*t; - отрицательный коэффициент при Et не может быть удовлетворительно интерпритирован.
Ct = 0,815042*Nit - 0,815042*Irt
Lft = 0,894735*Wgt + 1,25076*Pit
Wgt = 27,7484 - 0,0524772*Lft - 0,303455*(Pit – Pit-1); (низкий R2 и плохая интерпритируемость коэффициента при dPit);
Pit = -0,0182869*Yt-1 + 0,708228*Wgt + 0,163126*Mt + 0,0695614*Mt-1;
Unt = -0,00284301*Yt-1 + 0,133628*Lft
Imt = 0,109922*Ct
Ext = 0,0897042*Yt
It = 1775,2 – 1775,2*Irt;
Gt = 340,54;
Tt = 0,094212;
Irt = 0,494;
Mt = 1/(0,00354946 - 0,0000430572*t)
Nt = 1/(0,00993241 - 0,000050896*t)
Запишем систему, решенную относительно эндогенных пременных:
Lft =37.76321+0.4129764*Pit-1- 0.013687*Yt-1+0.1221*Mt+0.052064*Mt-1;
Pit = 15.021+ 0.164264*Pit-1- 0.0146333*Yt-1+0.1305345*Mt+0.055663*Mt-1;
Wgt = 21.2086+ 0.23194*Pit-1+0.0051588*Yt-1-0.04602*Mt- 0.01962*Mt-1;
Unt = -0,00284301*Yt-1 + 0,133628*Lft-1
Imt = 0,109922*Ct
Ext = 0,0897042*Yt
It = 1775,2 – 1775,2*Irt;
Gt = 340,54;
Tt = 0,094212;
Irt = 0,494;
Mt = 1/(0,00354946 - 0,0000430572*t)
Nt = 1/(0,00993241 - 0,000050896*t)
Ct = 0,815042*Nit - 0,815042*Irt
Yt = K1,28307Et-0,379686 e- 0,015362*t;
Сравнительная таблица по коэффициенту Тэйла:
|
Модель |
Коэффициент Тэйла. | |
|
1МНК. |
2МНК. | |
|
Y |
0,06072 |
0,08618 |
|
C |
0,06655 |
0,09177 |
|
Lf |
0,01754 |
0,02758 |
|
Wg |
0,05115 |
0,02817 |
|
Ex |
0,05880 |
0,08390 |
|
Im |
0,07824 |
0,10450 |
|
Un |
0,04435 |
0,05703 |
|
Pi |
0,01793 |
0,02190 |
Как видно из таблици – в большинстве случаев 1МНК дает лучший результат, чем 2МНК. По данной причине далее для целей прогнозирования мы будем использовать 1МНК. (хотя для сравнения построим прогноз и по 2МНК).
