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

Ольков_С_Г_Аналитическая юриспруденция

.pdf
Скачиваний:
280
Добавлен:
13.05.2015
Размер:
9 Мб
Скачать
☆

 

 

 

 

 

rx1x1rx1x2rx1x3

 

111

 

 

 

 

 

 

 

 

D

 

R

 

=

r

rer

 

=111t

=0.

 

 

 

 

 

 

 

x2x1

x2x2x2x3

 

 

 

 

 

 

 

 

rx3x1rx2x3rx3x3

 

111

 

Все остальные комбинации находятся в пределах между нулем и единицей. Чем ближе определитель матрицы парных коэффициентов корреляции к единице, тем выше мультиколлинеарность факторов.

Для выявления мультиколлинеарности факторов следует проводить проверку гипотезы о независимости переменных: H0:

Det|R|=1. Величина

é

-

1

ù

имеет приближенное

ên -1

6

×(2×m +5) ×lg DetR ú

 

 

ë

 

 

û

 

æ

1

× n ×(n -1)

ö

степенями свободы. Если фактическое

распределение χ2 с ç

2

÷

è

 

 

ø

 

 

 

значение χ2 больше табличного, то принимается альтернативная гипотеза, а нулевая отклоняется. Это значит, что Det|R|≠1. В противном случае, принимается нулевая гипотеза.

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

гомоскедастичности имеет место неравенство: σ 2ε x1 ¹ σ 2ε x2 ¹ σ 2ε x k . При гомоскедастичности, напротив, имеет место равенство дисперсий остатков.

276

При построении регрессионных моделей с помощью метода наименьших квадратов (МНК) важное значение имеет соблюдение определенных требований или условий относительно случайной составляющей εi = yi − yi , которая по определению всегда присутствует в регрессионных уравнениях. По существу, это разница между эмпирическими и теоретическими значениями результативного признака (зависимой переменной регрессионной модели). Она может быть установлена только после получения регрессионного уравнения и подстановки в него соответствующих значений факторной (факторных) переменной (переменных). При этом нельзя забывать, что остатки, как и параметры регрессионного уравнения – это лишь оценки того, что есть в реальности. То есть их следует рассматривать как определенные выборочные значения остатков, характеризующих генеральную совокупность остатков.

Требования, предъявляемые к случайному члену регрессионных уравнений:

1)Остатки должны носить случайный характер.

2)Средняя величина остатков должна равняться нулю.

3)Гомоскедастичность остатков (равный разброс для всех значений факторной переменной) или отсутствие гетероскедастичности (неравный разброс остатков).

4)Отсутствие автокорреляции остатков.

5)Подчинение остатков закону нормального распределения.

Соблюдение данных условий обеспечивает несмещенность, эффективность и состоятельность полученных с помощью МНК оценок.

Несмещенность говорит о том, что математическое ожидание остатков равно нулю. Эффективность о малой дисперсии остатков. Оценка является состоятельной, если с увеличением объема выборки (с ростом числа наблюдений) увеличивается точность оценки.

277

Именно поэтому недостаточно просто получить оценку параметров уравнения и проверить их с помощью соответствующих тестов. Кроме того, необходимо исследовать остатки.

Выявление гетероскедастичности при малых объемах выборки.

Для выявления гетероскедастичности при сравнительно небольшом числе наблюдений удобно использовать тест Гольдфельда-Квандта, названный так в честь его создателей Гольдфельда и Квандта. Тест выполняется по нижеследующему алгоритму:

1)По первичным объединенным данным получаем регрессионное уравнение y = a +bx .

2)Наблюдения (n) по независимой переменной (х) ранжируются от минимума к максимуму.

3)В соответствие этим наблюдениям приводятся значения зависимой переменной у и её оценки y , полученные по регрессионному уравнению: y = a +bx , а также остатки εi = yi − yi . Все это сводится в единую таблицу.

4)Вариационный ряд факторной переменной (x) разбивается на две равные части (первая - с малыми значениями фактора, а вторая - с большими). Далее из вариационного ряда исключается несколько значений, находящихся в центре данного ряда. Например, для ряда, включающего 22 наблюдения, можно удалить 4-5 наблюдений. При этом важно, чтобы соблюдалось неравенство:

n − Ц

> k , где k – число оцениваемых параметров, Ц -

2

 

количество удаленных центральных наблюдений, n – первоначальное число наблюдений.

5) Для

обеих выборок строим новые регрессионные

уравнения.

 

6)Ищем остаточную сумму квадратов отклонений

S= å( yi − yi )2 , для первой (S1) и второй выборки (S1).

7)Находим G, как частное от деления остаточных сумм

квадратов: G = SS1 .

2

278

8) Проверяем гипотезу H1 о наличии гетероскедастичности (отсутствии гомоскедастичности), сравнивая полученное эмпирическое значение G c табличным значением F-критерия с

числом степеней свободы df= n − Ц − 2k .

2

9) Интерпретируем полученный результат. Если G>F, то имеет место гетероскедастичность, и она тем больше, чем больше разница между G и F. Соответственно и наоборот, если G<F, то имеет место гомоскедастичность, и она тем больше, чем больше разница между G и F.

Напомним, что гомоскедастичность говорит о равенстве дисперсий остаточных величин и обеспечивает выполнение одной из предпосылок МНК (метода наименьших квадратов), что тоже самое - выполнение условия Гаусса-Маркова о том, что дисперсия случайного члена в каждом наблюдении должна быть постоянной. Имеется в виду возможное поведение случайного члена до того как сделана выборка, поскольку мы имеем дело не с самим случайным членом, а его оценками. Другими словами, вероятность того, что величина ε примет какое-то данное положительное или отрицательное значение будет одинаковой для всех наблюдений.

Рассмотрим конкретный пример. Пусть начальника следственного управления волнует вопрос о влиянии стажа следственной работы на число, совершаемых следователями уголовно-процессуальных правонарушений. Он собирает первичные статистические данные, и составляет исходную таблицу.

УПП, Стаж, штук месяцев

25 12

29 14

23 17

31 13

22 12

10 49

8 53

279

18

25

11

30

16

24

41

5

9

67

3

95

1

129

7

36

4

93

17

24

32

17

2

145

6

115

21

18

15

14

 

Далее.

Начальник

 

 

 

управления

получает

следственного

оценочное регрессионное уравнение.

 

 

 

 

 

Коэффициент

 

 

t-

P-

 

 

 

ы

Стандартна

 

статистик

Значени

 

 

 

 

 

я ошибка

 

а

е

 

 

 

 

 

 

 

 

 

 

 

Y-

 

 

 

 

 

 

 

 

пересечение

25,44014

2,166866

 

11,74052

2E-10

 

 

 

 

 

 

 

 

 

 

 

Перем

 

 

 

 

 

 

 

 

енная X

-0,20723

0,035065

 

-5,91

8,83E-06

 

 

 

 

 

 

 

 

 

 

Рис. Зависимость случайных остатков от факторной переменной.

280

Кроме того, можно построить график остатков от теоретических значений зависимой переменной y .

Рис. Зависимость случайных остатков от теоретических значений зависимой переменной.

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

Зная, что гетероскедастичность может привести к уменьшению эффективности оценки параметра b (возникают проблемы с использованием формулы стандартной ошибки коэффициента регрессии mb, предполагающей равную дисперсию остатков при каждом значении фактора). Принимается решение проверить модель на гетероскедастичность, для чего используется метод Гольфельда-Квандта.

1)По первичным объединенным данным получаем регрессионное уравнение y = a +bx . В нашем случае: y =25,44 −0,207 x .

2)Наблюдения (n) по независимой переменной (х) ранжируются от минимума к максимуму.

3)В соответствие этим наблюдениям приводятся значения зависимой переменной у и её оценки y , полученные по

регрессионному

уравнению:

y = a +bx (в

нашем

случае:

y =25,44 −0,207 x ), а

также

остатки

εi = yi − yi . Все это

сводится в

единую таблицу.

 

 

 

 

 

 

 

 

 

 

 

Х,

 

Y,

 

 

εi = yi − yi

 

 

 

 

 

 

 

 

 

 

 

 

 

281

 

 

 

 

 

Ф.И.О.

месяцы

 

Y ,

 

 

 

шт.

шт.

 

 

 

 

 

 

1

5

41

24,404

16,59601

 

 

 

 

 

2

12

25

22,953

2,04662

 

 

 

 

 

3

12

22

22,953

-0,95338

 

 

 

 

 

4

13

31

22,746

8,25385

 

 

 

 

 

5

14

29

22,539

6,46108

 

 

 

 

 

6

14

15

22,539

-7,53892

 

 

 

 

 

7

17

23

21,917

1,08277

 

 

 

 

 

8

17

32

21,917

10,08277

 

 

 

 

 

9

18

25

21,71

3,29

 

 

 

 

 

10

24

16

20,467

-4,46662

 

 

 

 

 

11

24

17

20,467

-3,46662

 

 

 

 

 

12

25

18

20,259

-2,25939

 

 

 

 

 

13

30

11

19,223

-8,22324

 

 

 

 

 

14

36

7

17,98

-10,9799

 

 

 

 

 

15

49

10

15,286

-5,28587

 

 

 

 

 

16

53

8

14,457

-6,45695

 

 

 

 

 

17

67

9

11,556

-2,55573

 

 

 

 

 

18

93

4

6,1678

-2,16775

 

 

 

 

 

19

95

3

5,7533

-2,75329

 

 

 

 

 

20

115

6

1,6087

4,39131

 

 

 

 

 

21

129

1

-1,293

2,29253

 

 

 

 

 

22

145

2

-4,608

6,60821

 

 

 

 

 

ИТОГО

1007

355

351

3,99753

 

 

 

 

 

4) Вариационный ряд факторной переменной (x) разбивается на две равные части (первая - с малыми значениями фактора, а вторая - с большими). Далее из вариационного ряда исключается несколько значений, находящихся в центре данного ряда. Например, для ряда, включающего 22 наблюдения, можно удалить 4-5 наблюдений. При этом важно, чтобы соблюдалось неравенство:

n − Ц

> k , где k – число оцениваемых параметров, Ц -

2

 

количество удаленных центральных наблюдений, n – первоначальное число наблюдений.

Для нашего примера: 22 −4 > 2 .

2

282

 

5) Для

обеих выборок строим

новые регрессионные

уравнения.

 

 

 

 

 

 

Для малых значений факторной переменной.

 

 

 

 

 

Станда

t-

P-

 

 

 

Коэффиц

ртная

статистик

Значени

 

 

 

иенты

ошибка

а

е

 

 

 

 

 

 

 

 

 

Y-

 

8,51306

4,7479

0,0

 

 

пересечение

40,42

8

94

02088

 

 

 

 

 

 

 

 

 

Перем

 

0,60601

-

0,1

 

 

енная X 1

-0,99

9

1,63361

46357

 

 

 

 

 

 

 

 

Для больших значений факторной переменной.

 

 

 

 

 

 

 

Станда

t-

P-

 

Коэффиц

ртная

статистик

Значени

 

иенты

ошибка

а

е

 

 

 

 

 

Y-

 

1,67507

6,9321

0,0

пересечение

11,61191

7

64

00225

 

 

 

 

 

Перем

 

0,01781

-

0,0

енная X 1

-0,0697

6

3,91235

05806

 

 

 

 

 

6)

Ищем

остаточную

сумму

квадратов

отклонений

S = å( yi

− yi )

 

, для первой (S1) и второй выборки (S1).

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Таблица. Для первой выборки с малыми значениями

факторной переменной.

 

 

 

 

 

 

 

 

Х,

 

 

Y,

Y ,

 

 

 

 

 

 

месяц

 

 

εi = yi − yi

 

εi

2

 

 

ы

 

 

 

 

шт.

шт.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

41

35,47

5,53

 

30,5809

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

25

28,54

-3,54

 

12,5316

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

22

28,54

-6,54

 

42,7716

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

 

31

27,55

3,45

 

11,9025

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

 

29

26,56

2,44

 

5,9536

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

 

15

26,56

-11,56

 

133,6336

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

17

 

23

23,59

-0,59

 

0,3481

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

17

 

32

23,59

8,41

 

70,7281

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

 

25

22,6

2,4

 

5,76

 

 

 

 

 

 

 

 

 

 

 

 

Сумма

 

 

 

 

 

 

314,21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

283

 

Таблица. Для

второй

выборки с

большими

значениями

факторной переменной.

 

 

 

 

 

 

 

Х,

 

Y,

 

Y ,

 

εi = yi − yi

 

εi 2

 

 

месяцы

 

шт.

 

шт.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

36

 

7

 

9,10271

 

-2,10271

 

4,421389

 

 

 

 

 

 

 

 

 

 

 

 

 

49

 

10

 

8,19661

 

1,80339

 

3,252215

 

 

 

 

 

 

 

 

 

 

 

 

 

53

 

8

 

7,91781

 

0,08219

 

0,006755

 

 

 

 

 

 

 

 

 

 

 

 

 

67

 

9

 

6,94201

 

2,05799

 

4,235323

 

 

 

 

 

 

 

 

 

 

 

 

 

93

 

4

 

5,12981

 

-1,12981

 

1,276471

 

 

 

 

 

 

 

 

 

 

 

 

 

95

 

3

 

4,99041

 

-1,99041

 

3,961732

 

 

 

 

 

 

 

 

 

 

 

 

 

115

 

6

 

3,59641

 

2,40359

 

5,777245

 

 

 

 

 

 

 

 

 

 

 

 

 

129

 

1

 

2,62061

 

-1,62061

 

2,626377

 

 

 

 

 

 

 

 

 

 

 

 

 

145

 

2

 

1,50541

 

0,49459

 

0,244619

 

 

 

 

 

 

 

 

 

 

 

 

 

Сумма

 

 

 

 

 

 

 

25,80213

 

 

 

 

 

 

 

 

 

 

 

 

7) Находим G, как частное от деления остаточных сумм

квадратов: G =

S1

. Для нашего примера: G =

S1

=

314,21

=12,177.

S2

S2

25,80213

 

 

 

 

8) Проверяем гипотезу H1 о наличии гетероскедастичности (отсутствии гомоскедастичности), сравнивая полученное эмпирическое значение G c табличным значением F-критерия с

числом степеней свободы df= n − Ц − 2k . Для нашего примера число

2

степеней свободы: df= 22 - 4 - 2×2 = 7 .

2

По таблице критических значений находим Fтаб. для 7 степеней свободы для каждой остаточной суммы квадратов. При уровне значимости 5% Fтаб.=3,79 (в данном случае берем 7 по числителю и по знаменателю βтаблицы).

9) Интерпретируем полученный результат. В нашем примере G>F . То есть имеет место гетероскедастичность, и она достаточно велика: 3,79<12,177.

Ответ: в исходной регрессионной модели, полученной начальником следственного управления, имеет место гетероскедастичность остатков, а, следовательно, модель нужно

284

совершенствовать, дабы получить более точную оценку параметра β (параметр b в регрессионном уравнении).

При построении регрессионных моделей с использованием МНК весьма важным условием является отсутствие автокорреляции остатков. То есть должна отсутствовать корреляция между остатками текущих (εi) и предыдущих значений (εi-1):

r = K (εi ,εi −1) .

Низкая

автокорреляция

обеспечивает

σεi

×σεi 1

 

 

 

 

−

 

 

 

состоятельность и эффективность, полученных оценок параметров регрессионного уравнения.

Для примера со следственным управлением имеем:

εi-1

εi

16,59601

2,04662

2,04662

-0,9534

-0,95338

8,25385

8,25385

6,46108

6,46108

-7,5389

-7,53892

1,08277

1,08277

10,0828

10,08277

3,29

3,29

-4,4666

-4,46662

-3,4666

-3,46662

-2,2594

-2,25939

-8,2232

-8,22324

-10,98

-10,9799

-5,2859

-5,28587

-6,457

-6,45695

-2,5557

-2,55573

-2,1678

-2,16775

-2,7533

-2,75329

4,39131

4,39131

2,29253

2,29253

6,60821

r= K (εi ,εi−1 ) = 0,4 .

σεi ×σεi−1

285