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

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

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

Когда исследуется корреляционная связь между переменными, полученными на основе пространственных рядов, то есть при фиксированном времени, у нас отсутствует риск получения так называемой ложной корреляции (spurious correlation) второго рода, которая может возникать при исследовании зависимостей по временным рядам. Назовем этот вид ложной корреляции ложной корреляцией второго рода, чтобы отграничить её от другого вида ложной корреляции, который назовем ложной корреляцией первого рода. Ложная корреляция первого рода возникает тогда, когда между переменными х и y существует высокий коэффициент корреляции, но реально независимая переменная не может детерминировать (обуславливать) поведение зависимой переменной, а поведение обоих переменных обусловлено действием «третьей» силы или связь случайна. Например, между умышленными убийствами и суицидами в РФ существует сильная корреляционная связь, хотя убийства не детерминируют суицидальное поведение, и, наоборот. Дело в том, что у убийств и суицидов, как разновидностей социально-патологического поведения, существует схожий причинный комплекс. Эндрю Ф. Сигл правильно замечает: «Одним из возможных оснований для существования «корреляции без причинной обусловленности» является наличие некоторого скрытого, ненаблюдаемого, третьего фактора, создающего впечатление, будто одна переменная является причиной другой переменной, тогда как на самом деле причиной для каждой из этих двух переменных является эта неизвестная третья переменная. Термином ложная корреляция обозначают высокую корреляцию, которая на самом деле обеспечивается действием некоего третьего фактора»26.

Ложная корреляция второго рода27 – это часть корреляции, которая обусловлена не чистым влиянием исследуемых в модели переменных, а эффектами, содержащимися в структуре изучаемых временных рядов. Очевидно, что временные ряды могут включать в

26 Сигел Э.Ф. Практическая бизнес-статистика: Пер. с англ. – М.: Издательский дом «Вильямс», 2002. С. 545.

27

216

себя тенденции, сезонные и циклические колебания, содержат случайную компоненту, что при определенных условиях негативно сказывается на коэффициенте корреляции – может быть завышен или занижен. Так, наличие тенденции в изучаемых временных рядах означает, что на переменные модели оказывает влияние фактор времени, оставшийся за пределами модели. Его влияние приведет к тому, что между остатками (ui) за текущий и предшествующий периоды времени возникнет корреляция (автокорреляция). То есть остатки не будут случайно распределены, а это ведет к нарушению важного условия для применения метода наименьших квадратов. Совершенно неслучайно существует потребность избавления, как от тенденции, так и других структурных компонентов временных рядов при проведении регрессионно-корреляционного анализа нацеленного на установление связи между переменными, поиск закономерностей. Если во временном ряду выявлена сезонная или иная периодическая компонента, то от неё можно избавиться с помощью процедуры сглаживания, использования сезонных индексов и т.д. От тенденции обычно избавляются включением в модель фактора времени, который «поглощает» часть корреляции или применяют иные методы, в частности метод последовательных разностей или метод отклонений от трендов, когда из временных рядов устраняется трендовая составляющая.

Покажем применение упомянутых методов на конкретных примерах.

I. Метод включения в модель регрессии фактора времени

является одним из наиболее эффективных, поскольку учитывает всю информацию, содержащуюся в исходных данных. Суть метода проста – строится новая модель, включающая дополнительную переменную - время (t).

Задача.

Дано: временные ряды переменных – коэффициент Джини (независимая переменная), умышленные убийства (шт.) (зависимая переменная).

Таблица №.

217

 

 

Коэффициент

Умышленные

Годы

t

Джини

убийства

 

по Российской

в Российской

 

 

 

 

Федерации

Федерации, ед.

1990

1

0,218

15600

1991

2

0,26

16200

1992

3

0,318

23000

1993

4

0,398

29200

1994

5

0,409

32300

1995

6

0,381

31700

1996

7

0,375

29400

1997

8

0,381

29300

1998

9

0,398

29600

1999

10

0,394

31100

Требуется: 1) проверить временные ряды на наличие в них тенденций; 2) найти коэффициент корреляции между переменными х (коэффициент Джини) и y (умышленные убийства); 3) найти параметры регрессионного уравнения до включения в модель фактора времени и дать их интерпретацию; 4) найти коэффициент корреляции между переменными после включения в модель фактора времени; 5) найти параметры регрессионного уравнения после включения в модель фактора времени и дать их интерпретацию.

1). Ниже представлен временной ряд коэффициентов Джини за период с 1990 по 1999 годы, из которого наглядно видно, что в ряду присутствует положительная тенденция. В среднем ежегодно коэффициент увеличивался от начала отсчета на 0,016.

218

Аналогичная ситуация наблюдается и во временном ряду умышленных убийств. Здесь присутствует положительная тенденция: в среднем ежегодно число умышленных убийств от начала отсчета увеличивалось на 1604 штук.

2. Коэффициент корреляции между переменными составил 0,975, что говорит о сильной положительной связи между переменными.

219

3.Очевидно, что свободный член в данном случае криминологического смысла не имеет и не интерпретируется. Коэффициент регрессии показывает, что при изменении коэффициента Джини на 0,1 число умышленных убийств в среднем изменяется на 9261,5 штук.

4.После включения в модель фактора времени коэффициент корреляции составил 0,976, то есть, практически не изменился.

5.Свободный член также не имеет криминологического смысла (отрицательная величина), а вот коэффициент регрессии снизился. При изменении коэффициента Джини на 0,1 число умышленных убийств стало в среднем изменяться на 8813,3 штук. То есть уменьшилось на 448,2 штуки. Это количество умышленных убийств определяется влиянием иных факторов, не учтенных в модели.

Таблица №. Регрессионная статистика для модели:

УУ=а+b∙кДж.

ВЫВОД

ИТОГОВ

Регрессионная

статистика

Множеств

 

енный R

0,976463

 

 

R-квадрат

0,95348

 

 

220

Нормиров

 

 

 

 

 

 

анный R-

 

 

 

 

 

 

квадрат

0,940188

 

 

 

 

 

 

 

 

 

 

 

 

Стандартн

 

 

 

 

 

 

ая ошибка

1528,533

 

 

 

 

 

 

 

 

 

 

 

 

Наблюден

 

 

 

 

 

 

ия

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дисперсио

 

 

 

 

 

 

нный

 

 

 

 

 

 

анализ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

Значим

 

 

f

SS

MS

F

ость F

 

 

 

 

 

 

 

 

Регрессия

2

3,35E+08

1,68E+08

71,73582

2,17E-05

 

 

 

 

 

 

 

 

Остаток

7

16354896

2336414

 

 

 

 

 

 

 

 

 

 

Итого

9

3,52E+08

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

Ста

t-

 

 

 

 

оэффиц

ндартная

статисти

P-

Нижни

Верхн

 

иенты

ошибка

ка

Значение

е 95%

ие 95%

 

 

 

 

 

 

 

Y-

 

 

 

 

 

 

пересечен

 

 

 

 

 

 

ие

-5083,71

3340,442

-1,52187

0,171854

-12982,6

2815,18

 

 

 

 

 

 

 

Переменна

1

 

 

0,6

-

751,671

я X 1

26,3373

264,4541

0,477729

47401

498,997

9

 

 

 

 

 

 

 

Переменна

 

 

 

 

 

116888,

я X 2

88133,8

12160,38

7,247621

0,00017

59379,08

5

 

 

 

 

 

 

 

II. Метод последовательных разностей – позволяет исключить тенденции во временных рядах путем перехода от исходных (первичных, эмпирических, фактических) данных к их абсолютным приростам: Δy= yi − yi−1 и Δх= xi −xi−1 . То есть в качестве переменных для регрессионно-корреляционного анализа берутся цепные абсолютные приросты. Если тренд выражен линейной функцией, то достаточно взять первые цепные абсолютные приросты. Если тренд выражен параболой второго порядка, то нужно брать цепные абсолютные приросты, но от первых цепных абсолютных приростов (вторые цепные абсолютные приросты). Когда тенденция выражена степенной или

221

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

Рассмотрим решение этой задачи на конкретном примере.

 

 

Коэффициент

Умышленные

 

Годы

 

Джини

по убийства

в

t

Российской

Российской

 

 

 

 

 

Федерации

Федерации, ед.

 

1990

1

0,218

15600

 

1991

2

0,26

16200

 

1992

3

0,318

23000

 

1993

4

0,398

29200

 

1994

5

0,409

32300

 

1995

6

0,381

31700

 

1996

7

0,375

29400

 

1997

8

0,381

29300

 

1998

9

0,398

29600

 

1999

10

0,394

31100

 

Коэффициент Δх=

−

Умышленные

Δy=

−

 

 

xi −xi 1

 

yi − yi 1

Джини

по

 

убийства

в

 

 

 

Российской

 

 

 

Российской

 

 

 

 

 

Федерации,

 

 

 

Федерации

 

 

 

 

 

ед.

 

 

 

 

 

 

 

 

 

 

0,218

0,042

15600

 

 

600

 

0,26

0,058

16200

 

 

6800

 

0,318

0,08

23000

 

 

6200

 

0,398

0,011

29200

 

 

3100

 

0,409

-0,028

32300

 

 

-600

222

Коэффициент

Δх=

xi −xi 1

Умышленные

Δy=

yi − yi 1

 

 

 

−

убийства в

 

−

Джини

по

 

 

 

 

 

 

Российской

 

 

Российской

 

 

 

 

 

 

Федерации,

 

 

Федерации

 

 

 

 

 

 

ед.

 

 

 

 

 

 

 

 

 

0,381

 

-0,006

31700

 

-2300

 

0,375

 

0,006

29400

 

-100

 

0,381

 

0,017

29300

 

300

 

0,398

 

-0,004

29600

 

1500

 

0,394

 

 

31100

 

 

Рис. №. Зависимость абсолютного прироста числа умышленных убийств от абсолютного прироста коэффициента Джини.

Таким образом, если абсолютный прирост по коэффициенту Джини меняется на 0,01, абсолютный прирост по числу умышленных убийств меняется на 721,26 (72126∙0,01=721,26) умышленных убийств. Видно, что за счет частичной потери информации заметно ухудшился коэффициент корреляции.

Таблица №. Итоговая регрессионная статистика к задаче.

ВЫВОД ИТОГОВ

Регрессионная статистика

Множественный R

0,801011

223

R-квадрат

0,641618

 

 

 

 

 

 

 

 

Нормированный

 

 

 

 

R-квадрат

0,590421

 

 

 

 

 

 

 

 

Стандартная

 

 

 

 

ошибка

1971,985

 

 

 

 

 

 

 

 

Наблюдения

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дисперсионный анализ

 

 

 

 

 

 

 

 

 

df

SS

MS

F

 

 

 

 

 

Регрессия

1

48734484

48734484

12,53225

 

 

 

 

 

Остаток

7

27221072

3888725

 

 

 

 

 

 

Итого

8

75955556

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t-

P-

 

Коэффициент

Стандартная

статистик

Значени

 

ы

ошибка

а

е

 

 

 

 

 

Y-пересечение

311,7655

768,6494

0,405602

0,697146

 

 

 

 

 

Переменная X 1

72125,63

20373,94

3,540092

0,009469

 

 

 

 

 

III. Метод отклонений от тренда связан с заменой исходных временных рядов (х и y) временными рядами их отклонений от полученных для них трендовых уравнений ( yi − yi и xi − xi ). То есть первоначально находим два уравнения: y = а +bt и x = a + bt ). Затем последовательно получаем оценки для переменных икс и игрек, подставляя в оба уравнения значения времени (1, 2 и т.д.). После этого из фактических данных по икс и игрек вычитаем их оценки – находим остатки для икс и игрек, а далее строим зависимость остатков по игрек от остатков по икс.

Пример.

 

Тренд для

временного ряда коэффициента Джини

описывается уравнением: кДж=0,016t+0,26.

224

Рис. №. Тренд во временном ряду коэффициентов Джини за период с 1990 по 1999 годы.

Тренд для временного ряда умышленных убийств описывается уравнением: УУ=1604t+17913.

Рис. №. Тренд во временном ряду умышленных убийств за период с 1990 по 1999 годы.

 

 

 

 

 

KДж

УУ

 

 

 

(X)

(Y)

t

xi − xi

yi − yi

 

 

 

 

 

0,218

15600

1

-0,058

-3917

 

 

 

 

 

0,26

16200

2

-0,032

-4921

 

 

 

 

 

0,318

23000

3

0,01

275

 

 

 

 

 

225