Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

cgiirbis_64 (6)

.pdf
Скачиваний:
10
Добавлен:
28.01.2022
Размер:
4.3 Mб
Скачать

пояснюючих змінних (F-критерій); в) кожної пари пояснюючих змінних (t-критерій).

Усі ці критерії при порівнянні з їх критичними значеннями дають змогу зробити конкретні висновки щодо наявності чи відсутності мультиколінеарності пояснюючих змінних.

Опишемо алгоритм Фаррара—Глобера.

Крок 1. Нормалізація змінних.

Нехай Х1, Х2, Х3, ..., Хm — вектори пояснювальних змінних економетричної моделі. Елементи стандартизованих векторів обчислимо за формулою:

 

X ik

 

 

 

 

X ik

 

k

,

X ik*

X

k

, або

X ik*

X

 

 

 

 

 

 

 

n 2 X k

 

X k

де п — число спостережень i 1, n ; т — число пояснювальних змінних k 1, m ; X k — середня арифметична k-ї пояснювальної змінної; 2X k — дисперсія k-ї пояснювальної змінної.

Крок 2. Знаходження кореляційної матриці першого порядку

r X * X * ,

де X * матриця нормалізованих пояснювальних змінних; X * матриця, транспонована до матриці X * .

Крок 3. Визначення критерію «хі-квадрат»:

2 n 1 1 2m 5 ln r ,6

де |r| — визначник кореляційної матриці r.

Значення цього критерію порівнюється з табличним за таких умов: 12 m m 1 ступенів свободи і рівень значущості . Якщо

факт2 2табл , то в масиві незалежних змінних мультиколінеарність відсутня.

Крок 4. Визначення матриці, оберненої до r:

Cr 1 X * X * 1 .

29

1 rkj2

Крок 5. Розрахунок F-критеріїв:

Fk ckk 1 nm m1 ,

де ckk — діагональні елементи матриці С. Фактичні значення критеріїв Fk порівнюються з табличними, коли маємо т – 1 і п – т

ступенів свободи і рівень значущості α. Якщо Fk факт Fтабл , від-

повідна k-та пояснююча змінна мультиколінеарна з іншими. Коефіцієнт детермінації для кожної змінної обчислюється так:

R2

1

1

.

 

X k

 

ckk

 

 

Крок 6. Знаходження частинних коефіцієнтів кореляції:

rkj

ckj

,

ckk * c jj

 

 

де сkj елемент матриці С, що міститься в k-му рядку і j-му стовпці ( k 1, m , j l, m , ckk і cjj — діагональніелементи матриціС).

Крок 7. Розрахунок t-критеріїв:

tkj rkj n m .

Фактичні значення критеріїв tkj порівнюються з табличними, коли маємо nm ступенів свободи і рівень значущості α. Якщо

tkj факт tтабл , між пояснюючими змінними Xk і Xj існує мультиколінеарність.

Дослідження наявності мультиколінеарності згідно з алгоритмом Фаррара—Глобера

Розглянемо застосування алгоритму Фаррара—Глобера для розв’язування конкретної задачі.

Приклад. Оскільки між пояснюючими змінними можлива тісна залежність (мультиколінеарність), метод 1МНК може дати зміщені кількісні характеристики взаємозв’язку. Тому дослідимо наявність мультиколінеарності між пояснюючими змінними, скориставшись алгоритмом Фаррара—Глобера.

30

Якщо виявлено мультиколінеарність, потрібно перетворити певним чином вихідну інформацію, щоб уникнути мультиколінеарності, а далі знову оцінити параметри моделі методом 1МНК.

Для дослідження мультиколінеарності наведемо статистичну сукупність спостережень (табл. 1).

 

 

 

 

Таблиця 1

 

 

 

 

 

Питомі інвестиції Х3,

 

Місяць

Прибуток на

Фондовіддача

Продуктивність

 

місяць Y, грн.

Х1, грн.

праці Х2, грн.

грн.

 

 

 

1-й

40

12

5

15

 

 

 

 

 

 

 

2-й

45

17

7

18

 

 

 

 

 

 

 

3-й

40

13

6

16

 

 

 

 

 

 

 

4-й

43

14

7

17

 

 

 

 

 

 

 

5-й

48

16

6

20

 

 

 

 

 

 

 

6-й

39

15

5

15

 

 

 

 

 

 

 

7-й

42

14

6

16

 

 

 

 

 

 

 

8-й

45

17

9

18

 

 

 

 

 

 

 

9-й

38

12

5

19

 

 

 

 

 

 

 

10-й

48

18

10

20

 

 

 

 

 

 

 

11-й

50

20

11

22

 

 

 

 

 

 

 

12-й

48

17

10

21

 

 

 

 

 

 

 

13-й

49

18

12

21

 

 

 

 

 

 

 

14-й

45

19

8

20

 

 

 

 

 

 

 

15-й

49

20

9

22

 

 

 

 

 

 

 

16-й

52

22

14

23

 

 

 

 

 

 

 

17-й

54

24

15

24

 

 

 

 

 

 

 

18-й

51

21

13

20

 

 

 

 

 

 

 

19-й

55

25

16

24

 

 

 

 

 

 

 

20-й

56

27

18

25

 

 

 

 

 

 

 

Необхідно:

1. Обчислити середнє значення та стандартні відхилення.

31

2.Нормалізувати пояснюючі змінні.

3.Знайти кореляційнуматрицюr.

4.Визначити визначник матриці r (det r).

5.Обчислити критерій 2 .

6.Визначити F-критерії.

7.Визначити частинні коефіцієнти кореляції r.

8.Визначити t-критерії.

9.Зробити висновки щодо мультиколінеарності.

Розв’язання

Для розв’язання використаємо стандартні функції Excel. 1.Обчислимо середні значення та стандартні відхилення пояс-

нюючих змінних Х1, Х2, Х3.

Для цього в мастер функції знайдемо категорію «статистичні» і функції «СРЗНАЧ» та «СТАНДВІДХИЛ».

Х1

Х2

Х3

 

 

 

18,05

9,6

19,8

 

 

 

4,260899

3,965642

3,036619

 

 

 

2.Нормалізуємо пояснюючі змінні.

Серед статистичних функцій знайдемо функцію «НОРМАЛІЗАЦІЯ» та нормалізуємо Х1, Х2, Х3:

32

 

*

*

*

 

 

X1

X 2

X 3

 

 

1,41989

1,15996

1,58071

 

 

0,24643

0,65563

0,59276

 

 

1,1852

0,9078

1,25139

 

 

0,9505

0,65563

0,92208

 

 

 

 

0,48112

0,9078

0,065963

 

 

0,71581

1,15996

1,58071

 

 

0,9505

0,9078

1,25139

 

 

0,24643

0,1513

0,59276

 

 

1,41989

1,15996

0,26345

X *

0,01173

0,100866

0,065863

.

 

0,45765

0,353032

0,72449

 

 

 

 

0,24643

0,100866

0,395176

 

 

0,01173

0,605198

0,395176

 

 

0,222958

0,40347

0,065863

 

 

0,45765

0,1513

0,72449

 

 

0,927034

1,10953

1,053804

 

 

1,396419

1,361696

1,383117

 

 

0,692342

0,857364

0,065863

 

 

1,631111

1,613862

1,383117

 

 

 

 

2,100496

2,118194

1,712431

 

Транспонуємо матрицю Х* (нормалізовану) в матрицю X *

 

 

1,4199

0,2464

1,1852

0,9505

0,4811

0,7158

0,9505

 

 

 

1,16

0,6556

0,9078

0,6556

0,9078

1,16

0,9078

X *

 

 

 

 

1,5807

0,5928

1,2514

0,9221

0,06586

1,5807

1,2514

 

 

 

 

 

 

 

 

Продовження матриці X *

 

 

 

 

 

0,2464

1,4199

0,0117

0,45765

0,2464

0,0117

0,22296

0,45765

0,1513

1,16

0,10087

0,35303

0,10087

0,6052

0,4035

0,1513

0,5928

0,2635

0,06586

0,72449

0,39518

0,39518

0,06586

0,72449

 

 

 

 

 

 

 

 

 

 

Завкінчення матриці X *

 

 

 

 

 

 

 

 

0,92703

 

1,39642

0,69234

1,63111

2,1005

 

1,10953

 

1,3617

0,85736

1,61386

2,11819

.

1,0538

 

1,38312

0,06586

1,38312

1,71243

 

 

 

 

 

Перемножимо матриці X *

та Х*:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19

 

17,8964552

16,9413894

 

 

 

 

 

 

19

16,6415575

 

X * X *

17,8964552

 

.

 

16,9413894

16,6415575

 

19

 

 

3.Знайдемо кореляційну матрицю r.

33

Для знаходження кореляційної матриці r необхідно кожний

 

 

 

 

1

 

 

 

 

 

n 1 (у нашому випадку

елемент матриці X * X * помножити на

 

n – 1 = 19, оскільки п = 20):

 

 

 

 

 

 

1

0,941918693

0,891652074

 

r

0,941918693

1

0,875871449

.

 

0,891652074

0,875871449

1

 

 

 

4.Знайдемо визначник матриці r (det r).

Для знаходження det r необхідно серед математичних функцій 6Excel знайти функцію «МОПРЕД». Скориставшись нею, дістанемо:

det r

0,02182033

 

 

Оскільки det r наближається до 0, то в масиві пояснюючих змінних може існувати мультиколінеарність.

Прологарифмуємо визначник матриці r:

ln det r

–3,824913185

 

 

5.Обчислимо критерій 2

за формулою:

 

 

 

2

 

 

 

 

1

 

 

 

 

 

n 1

6

2m 5 ln det r .

Маємо:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

64,38603861

 

Знайдене

значення

2

порівнюємо з табличним значенням

kp

2 7,80 ,

коли маємо

1 m m 1 3 ступенів свободи та при

 

 

 

 

 

 

2

 

 

 

 

 

рівні значущості = 0,05.

34

Оскільки факт2 2крит , то в масиві пояснюючих змінних (про-

дуктивність праці, питомі інвестиції та фондовіддача) існує мультиколінеарність.

6.Обчислимо F-критерії.

Для визначення F-критеріїв необхідно знайти матрицю С, яка є оберненою до матриці r:

 

10,67120467

7,375970017

3,05460023

C

7,375970017

9,392918454

1,65019013

 

 

3,054600232

1,65019013

5,168995053

 

 

 

Безпосередньо F-критерії обчислюються за формулою:

FCkk 1 n m ,

m 1

де сkk — діагональний елемент матриці С.

F1

51,57975822

 

 

F2

44,76223175

 

 

F3

22,23464028

 

 

Обчислені F-критерії порівнюються з табличним значенням (F = 5,29), коли є n – m (16) ступенів свободи та при рівні значущос-

ті = 0,05.

У розглядуваному випадку F1 > Fкрит, F2 > Fкрит, F3 > Fкрит. Це означає, що кожна з пояснюючих змінних мультиколінеарна з іншими.

7.Визначимо частинні коефіцієнти кореляції r.

Частинні коефіцієнти кореляції показують тісноту зв’язку між двома пояснюючими змінними за умови, що інші змінні не впливають на цей зв’язок і обчислюються за формулою:

rkj.s

ckj

.

ckk c jj

 

 

r12.3

0,736736

 

 

r13.2

0,411287

 

 

r23.1

0,236827

 

 

 

35

Отже, спираючись на здобуті нами значення окремих (частинних) коефіцієнтів кореляції, можна сказати, що зв’язок між фондовіддачею та продуктивністю праці є тісним, якщо не враховувати вплив питомих інвестицій; зв’язок між фондовіддачею та питомими інвестиціями є слабким, якщо не брати до уваги вплив продуктивності праці. Зв’язок між продуктивністю праці та питомими інвестиціями також є слабким, якщо не враховувати фондовіддачу.

8.Визначимо t-критерії.

Ці критерії застосовуються для визначення мультиколінеарності двох пояснюючих змінних і обчислюються за формулою:

tkj

rkj.s

n m

 

1

r

2

 

 

 

 

 

 

 

 

 

 

kj.s

 

 

 

 

 

t12

 

 

6,445364

 

 

 

 

t13

 

 

1,980092

 

 

 

 

t23

 

 

1,003595

 

 

 

 

 

 

 

Обчислені t-критерії порівнюються з табличним значенням (t = 1,746), коли маємо n – m (16) ступенів свободи та при рівні значущості = 0,05.

Оскільки t12 > tкрит, то продуктивність праці та фондовіддача є відповідно мультиколінеарними між собою;

t13 > tкрит, тому відповідно продуктивність праці та питомі інвестиції є мультиколінеарними між собою;

t23 < tкрит, тому продуктивність праці та питомі інвестиції не є мультиколінеарними між собою.

Дослідження показали, що мультиколінеарність існує. Отже, для того, щоб можна було застосувати метод 1МНК для оцінювання параметрів моделі за цією інформацією, необхідно звільнитися від мультиколінеарності.

Покажемо, як побудувати економетричну модель методом, 1МНК на основі нормалізованих даних.

36

 

 

 

*

*

*

 

 

 

 

 

 

 

X1

X 2

X 3

 

 

 

 

 

 

 

1,41988804

1,15996353

1,58070558

 

 

1,30559111

 

 

 

 

0,24642685

0,65563156

0,59276459

 

 

0,3526049

 

 

 

 

 

 

 

 

 

 

1,1851958

0,90779755

1,25139192

 

 

1,30559111

 

 

 

 

0,95050356

0,65563156

0,92207826

 

 

0,73379939

 

 

 

 

0,48111909

0,90779755

0,06586273

 

 

0,21918683

 

 

 

 

0,71581133

1,15996353

1,58070558

 

 

1,49618836

 

 

 

 

0,95050356

0,90779755

1,25139192

 

0,92439663

 

 

 

0,24642685

0,15129959

0,59276459

 

 

0,3526059

 

 

 

 

1,41988804

1,15996353

0,26345093

 

1,6876856

 

X

*

 

0,01173461

0,100866394

0,06586273

 

 

0,21918683

 

 

 

0,457649864

0,353032379

0,72449006

.

Y *

0,60038132

.

 

 

 

 

 

 

 

 

 

0,24642685

0,100866394

0,3951764

 

 

0,21918683

 

 

 

 

0,01173461

0,605198364

0,3951764

 

 

0,40978407

 

 

 

 

0,222957626

0,40346558

0,06586273

 

 

0,3526049

 

 

 

 

0,457649864

0,151299559

0,72449006

 

 

0,40978407

 

 

 

 

0,927034339

1,109530335

1,05380372

 

 

 

 

 

 

1,396418814

1,36169632

1,38311738

 

 

0,9815758

 

 

 

 

 

 

1,36277029

 

 

 

 

0,692342101

0,85736435

0,06586273

 

 

0,79097856

 

 

 

 

1,631111052

1,613862305

1,38311738

 

 

1,55336753

 

 

 

 

2,100495528

2,118194275

1,71243105

 

 

1,74396477

 

 

 

 

1,631111052

1,613862305

1,38311738

 

 

 

 

 

 

 

2,100495528

2,118194275

1,71243105

 

 

 

 

 

 

 

 

 

 

 

Застосуємо функцію «лінійн» для знаходження всіх оцінок параметрів нової моделі, стандартні похибки цих оцінок, коефіцієнт детермінації, дисперсію залишків і т. ін.

b3

b2

b1

 

0,3183466

0,31026984

0,39186698

0

 

 

 

 

0,1529347

0,2061594

0,21974024

#Н/Д

 

 

 

 

0,9269234

0,2932106

#Н/Д

#Н/Д

 

 

 

 

71,877548

17

#Н/Д

#Н/Д

 

 

 

 

18,538468

1,46153177

#Н/Д

#Н/Д

 

 

 

 

Обчислюємо * — залежну змінну, підставивши пояснюючі

Y

змінні у відповідну формулу:

*

*

*

*

Y

b1 X1

b2 X 2

b3 X 3 .

37

1,41952111

0,48869381

1,1444776

0,86943409

0,44922971

1,14361669

1,05250946

0,33221481

1,00017763

0,04766457 . Y * 0,51951208

0,0605323

0,308979440,01684629

0,36303308

1,043002731,41001437

0,55828748

1,58022201

2,02547315

Щоб обчислити оцінки параметрів моделі в абсолютному виразі, скористаємося формулами:

 

 

Y

 

 

 

 

 

 

X 2 ;

 

 

 

 

X 3

;

a1

 

b1

 

X1 ; a2

b2 Y

 

a3

 

b3

Y

 

 

 

 

 

a0

Ycp a1 X cp1 a2 X cp2

a3 X cp3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

23,30887

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

 

0,482526

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

0,410496

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a3

 

 

0,550039

 

 

 

 

 

 

 

Отже, економетрична модель має вигляд:

 

 

0,410X 2

0,550X 3 ,

 

Y 23,31 0,483X1

де

— дохід; Х1 — фондовіддача; Х2 — продуктивність праці;

Y

Х3

— питомі інвестиції.

 

 

Завдання для самостійної роботи

Нехай на витрати обігу впливають: обсяг вантажообігу, грн., запаси з вантажообігу, грн., та трудомісткість одиниці вантажообігу, людино-год. Щоб побудувати економетричну модель цієї зале-

38

Соседние файлы в предмете Моделирование