Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Эээ Елена / CourseWork / Курсовая по ПСЭР1.doc
Скачиваний:
9
Добавлен:
19.04.2013
Размер:
333.82 Кб
Скачать

Двухшаговый метод наименьших квадратов.

Записанные выше уравнения, полученные по МНК, представляют собой структурную форму модели. Для получения оценок парметров модели по 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);

  1. Yt = 599,324*Irt + 3,51875*Mt - 2,20675*Mt-1 + 0,680562*Yt-1 - 7,64033*t;

  2. Ct=– 1832,52*Tt + 478,055*Irt + 2,7986*Mt - 1,8586*Mt-1 + 0,574976*Yt-1 - 4,48*t;

  3. Lft = -0,384668*Mt-1 + 2,53119*Pit-1 + 0,727815*Nt;

  4. It = 2,62736*Gt + 20,8567*Pit-1 - 0,352074*It-1 - 7,53186*Nt

  5. Wg = 0,0259616*G - 0,0506729*M + 0,0139026*Y1 + 0,186522*t;

  6. Ex = 232,531*Ir – 1,37859*M1 + 10,8018*Pi1 - 5,71101*t;

  7. Im = 106,071*Ir + 0,512589*M + 0,776596*M1 - 5,70969*Pit-1

  8. Un = 0,0107602*G - 0,0223452*M + 0,253997*Pit-1 - 0,00188577*Yt-1;

  9. Pit = 7,56143*Irt + 0,120031*Mt - 0,0853*Mt-1 + 0,852586*Pit-1 - 0,0032*Yt-1;

Шаг 2.

  1. Ln(Yt) = 0,825741*Ln(Kt) + 0,402234*Ln(Ett) - 0,0204713*t;

  2. Ct = 0,815648*Nit - 0,815648*Ir;

  3. Lft = 0,767357*Wgt + 1,30571*Pit;

  4. It = 1819,06-1819,06*Ir;

  5. Wgt = 27,7255 – 0,0514576*Lft - 0,339164*(Pit – Pit-1);

  6. Ext = 0,0896857*Yt

  7. Imt = 0,110128*Ct

  8. Unt = – 0,00247893*Yt-1 + 0,125305*Lft;

  9. 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МНК).

Соседние файлы в папке CourseWork