Ольков_С_Г_Аналитическая юриспруденция
.pdf
Москва |
212050 |
10509000 |
2018 |
|
|
|
|
Московская |
|
|
|
область |
120995 |
6713000 |
1802 |
|
|
|
|
Санкт-Петербург |
77146 |
4582000 |
1684 |
|
|
|
|
Ленинградская |
|
|
|
область |
32254 |
1632000 |
1976 |
|
|
|
|
Республика |
|
|
|
Татарстан |
70623 |
3769000 |
1874 |
|
|
|
|
Республика |
|
|
|
Башкортостан |
75471 |
4057000 |
1860 |
|
|
|
|
ХМАО |
37372 |
1520000 |
2459 |
|
|
|
|
ЯНАО |
10052 |
544000 |
1848 |
|
|
|
|
Курганская |
|
|
|
область |
24922 |
953000 |
2615 |
|
|
|
|
Кировская |
|
|
|
область |
23581 |
1401000 |
1683 |
|
|
|
|
Чувашская |
|
|
|
республика |
21653 |
1279000 |
1693 |
|
|
|
|
Республика |
|
|
|
Марий Эл |
17787 |
700000 |
2541 |
|
|
|
|
Республика |
|
|
|
Бурятия |
30273 |
961000 |
3150 |
|
|
|
|
Калининградская |
|
|
|
область |
19756 |
937000 |
2108 |
|
|
|
|
Смоленская |
|
|
|
область |
22840 |
974000 |
2345 |
|
|
|
|
Ярославская |
|
|
|
область |
26650 |
1310000 |
2034 |
|
|
|
|
Требуется: найти бисериальный коэффициент корреляции и ответить на вопрос, о коррелированности уровней преступности по субъектам РФ в зависимости от уровня их благосостояния.
Решение:
1. От таблицы первичных статистических данных перейдем к рабочей таблице.
Рабочая таблица.
Субъект РФ |
X |
Y |
|
|
|
|
|
|
Москва |
1 |
2018 |
|
|
|
Московская |
|
|
область |
1 |
1802 |
|
|
|
Санкт-Петербург |
1 |
1684 |
|
|
|
Ленинградская |
1 |
1976 |
|
|
|
266
область |
|
|
|
|
|
Республика |
|
|
Татарстан |
1 |
1874 |
|
|
|
Республика |
|
|
Башкортостан |
1 |
1860 |
|
|
|
ХМАО |
1 |
2459 |
|
|
|
ЯНАО |
1 |
1848 |
|
|
|
Курганская |
|
|
область |
0 |
2615 |
|
|
|
Кировская |
|
|
область |
0 |
1683 |
|
|
|
Чувашская |
|
|
республика |
0 |
1693 |
|
|
|
Республика |
|
|
Марий Эл |
0 |
2541 |
|
|
|
Республика |
|
|
Бурятия |
0 |
3150 |
|
|
|
Калининградская |
|
|
обл. |
0 |
2108 |
|
|
|
Смоленская |
|
|
область |
0 |
2345 |
|
|
|
Ярославская |
|
|
область |
0 |
2034 |
|
|
|
2. Найдем значения переменных, входящих в формулу для расчета бисериального коэффициента корреляции.
Y1 =(2018+1802+1684+1976+1874+1860+2459+1848)/8=1940 Y0 =(2615+1683+1693+2541+3150+2108+2345+2034)/8=
2271
y =(1940+2271)/2=2106.
|
|
1 |
N |
|
|
|
sy = |
× å(yi |
- y)2 = 399 |
||||
N |
||||||
|
|
i=1 |
|
|
||
Р.S. Промежуточное вычисление:
((2018 |
−2106 )2 +(1802 |
−2106 )2 +... |
+(2034 |
−2106 )2 ) /16 |
. |
N=16; n0 =8; n1 =8.
267
2. Рассчитаем |
|
|
бисериальный коэффициент |
корреляции по |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
n1 ×n0 |
|
|
|
|||
формуле: R = Y1 |
-Y0 × |
|
|
|
|
|
||||||||||
|
б |
|
|
sy |
|
|
|
|
|
N ×(N -1) |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
R =1940 -2271 |
|
|
|
|
|
|
|
|
|
|||||||
× |
8 ×8 |
|
|
|
|
|
|
|||||||||
б |
|
399 |
|
|
|
|
16 ×(16 -1) |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
R = -0,82957 × |
|
|
|
|
|
|
R = -0,82957 ×0,5163977795 |
= -0,428 . |
||||||||
0,2666666667 |
|
|
|
|
||||||||||||
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б |
|
Поскольку, как мы знаем, знак при бисериальном коэффициенте корреляции смысла не имеет, запишем: Rб = 0,428 .
Проверим, значим ли статистически полученный бисериальный коэффициент корреляции. Для этого рассчитаем tэмп.
по формуле: |
tэмп . = |
|
-0,428 |
|
× |
16 -2 |
|
=1,772 |
. |
|
|
|
|||||||||
|
|
|
1-0,428 ×0,428 |
|||||||
ì |
2,1 4 д |
л αя= |
0 |
, |
tт а=бí |
1 , 3 4д |
5л αя= |
0 |
, |
î |
Критическое значение t-критерия Стьюдента в соответствии с таблицей при уровне значимости α=0,05 при 14 степенях свободы составляет 2,14, что говорит о статистической не значимости бисериального коэффициента корреляции на этом уровне значимости. Между тем если уменьшить требовательность и согласиться на вероятность большей ошибки, например, принять уровень значимости α=0,2, то табличное критическое значение составит 1,345, и тогда бисериальный коэффициент корреляции становится значимым.
Ответ: Rб = 0,428 . Бисериальный коэффициент корреляции показывает, что между изучаемыми переменными имеется слабая корреляционная связь. То есть с вероятностью ошибки 20% можно утверждать, что уровень благосостояния субъекта РФ все же коррелирован с уровнем преступности.
Рангово-бисериальный коэффициент корреляции
используется в тех случаях, когда переменная игрек измерена по ранговой шкале. В данном случае расчет производится по формуле:
268
|
|
|
|
|
|
|
|
|
|
|
|
|
(Y1 |
-Y0 ) ×2 |
|
|
|
|
|
||||
Rбр = |
, где Y1 |
- |
средний ранг по элементам |
||||||||
|
|
N |
|||||||||
|
|
|
|
|
|
|
|
||||
переменной Y, которым соответствуют значения единицы (признак единица) по переменной X, Y0 - средний ранг по элементам переменной Y, которым соответствуют значения нуль (признак нуль) по переменной X, N – число наблюдений по переменной X.
Статистическая значимость рангово-бисериального коэффициента корреляции проверяется с помощью t-критерия Стьюдента с числом степеней свободы: df=N-2. Для этого
первоначально вычисляется значение |
|
|
|
|
N - 2 |
|
|
|
|
|
|||||
tэмп . = |
Rбр |
× |
1- R 2 . |
||||
|
|
бр |
|
||||
|
|
|
|
|
|||
Корреляционное отношение Пирсона.
Вообще говоря, все коэффициенты корреляции, как то: параметрический - линейный коэффициент корреляции или непараметрические - коэффициенты Спирмена и Кендалла, бисериальный коэффициент корреляции, коэффициенты ассоциации и контингенции, коэффициент конкордации, коэффициент Чупрова-Пирсона выявляют линейные зависимости между переменными, в то время как корреляционное отношение Пирсона нацелено на выявление нелинейной зависимости между переменными. Корреляционное отношение Пирсона исследует связь между двумя переменными в перекрестном порядке относительно друг друга, и всегда представляется двумя
коэффициентами: |
η |
и |
η |
. Если переменные Х и Y связаны |
|
хy |
|
yx |
|
строгой линейной зависимостью, то имеет место равенство: ηхy =ηyx
и в таком случае оба упомянутых коэффициента равны по величине линейному коэффициенту корреляции Пирсона, то есть
ηхy =ηyx = r = |
Kx, y |
|
|
|
|
|
σ x ×σ y . |
|
|
||||
|
|
|
||||
Формулы |
для вычисления корреляционного отношения |
|||||
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
Пирсона: ηхy = |
å f y (xy − x)2 |
|
|
|||
N |
|
|
||||
|
|
|
|
i=1 |
(1.1) |
|
|
|
|
|
å fx (xi − x)2 |
||
|
|
|
|
i=1 |
|
|
269
|
|
N |
|
|
|
|
|
|
η yx = |
å fx (yx |
− y)2 |
|
|||||
N |
|
|
(1.2), где х и y - средние арифметические по |
|||||
|
|
i=1 |
|
|
|
|
|
|
|
|
å f y (yi |
− y)2 |
|
||||
|
|
i=1 |
|
|
|
|
|
|
переменным х и y соответственно; х y и |
yx - групповые средние |
|||||||
арифметические; |
f y и f x - частоты в рядах Х и Y. |
|||||||
Решим конкретную задачу. Пусть |
нам нужно установить, |
|||||||
существует ли связь между сроком лишения свободы (Х), измеряемом в годах, и уровнем злокачественной агрессивности личности (Y), измеряемой в баллах от нуля до двадцати. Криминолог обследует 12 ранее судимых к различным срокам лишения свободы граждан на предмет их злокачественной
агрессивности. |
|
|
|
|
|
|
|
|
|
|
||
|
Дано: таблица исходных данных. |
|
|
|
|
|
||||||
|
Таблица №1. |
|
|
|
|
|
|
|
|
|
||
№ |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
Х |
4 |
6 |
6 |
10 |
5 |
3 |
10 |
2 |
2 |
7 |
6 |
7 |
Y |
10 |
12 |
17 |
14 |
11 |
10 |
19 |
9 |
8 |
19 |
12 |
14 |
|
Решение: |
|
|
|
|
|
|
|
|
|
|
|
Х |
Таблица №2. Для ранжированной переменной X. |
|
|
|||||||||
2 |
2 |
3 |
4 |
5 |
6 |
6 |
6 |
7 |
7 |
10 |
10 |
|
Y |
9 |
8 |
10 |
10 |
11 |
12 |
17 |
12 |
19 |
14 |
14 |
19 |
Таблица №3. Для значений, используемых в формулах.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
1 |
1 |
|
1 |
3 |
2 |
2 |
||||
f x |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
7 |
10 |
|
хi |
|
|
|
|
|
|
|
|
|
|
|
|
|
8,5 |
|
10 |
|
10 |
|
11 |
|
13,6 |
16,5 |
16,5 |
|
yx |
|
|
|
|
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
å fx = 2 +1+1+1+ 3 + 2 + 2 =12 . |
|
|
|
|
|
|
||||||
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
yx2 |
= |
9 +8 |
= 8,5 ; yx3 |
=10 ; yx 4 =10 |
; yx5 =11 ; |
|
|
|||||
|
|
|
||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
yx6 |
= 12 +17 +12 =13,6 ; yx7 = 19 +14 |
=16,5 |
; |
|
|
|
||||||
|
|
|
3 |
|
|
|
2 |
|
|
|
|
|
270
|
yx10 |
= |
14 +19 =16,5 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица №4. Для ранжированной переменной Y. |
|
|
|
||||||||||||||||||||
Y |
8 |
|
|
9 |
|
10 |
10 |
|
11 |
|
|
12 |
|
12 |
14 |
|
14 |
17 |
19 |
|
19 |
|||
X |
2 |
|
|
2 |
|
3 |
4 |
|
|
5 |
|
|
6 |
|
6 |
7 |
|
10 |
6 |
7 |
|
10 |
||
|
Таблица №5. Для значений, используемых в формулах. |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
1 |
|
1 |
2 |
|
1 |
2 |
|
2 |
1 |
|
2 |
|
|||||||||
f y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
8 |
|
|
9 |
|
10 |
|
|
11 |
|
12 |
|
14 |
|
17 |
|
19 |
|||||
yi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
3,5 |
|
5 |
|
6 |
|
8,5 |
|
6 |
|
8,5 |
||||||
x y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
å f y =1+1+ 2 +1+ 2 + 2 +1+ 2 =12 . |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xy8 =2 ; xy9 =2 ; xy10 |
= |
3 +4 |
= 3,5 ; xy11 |
=5 ; |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xy12 |
= |
6 +6 = 6 ; xy14 |
= 7 +10 |
= 8,5 ; |
xy17 |
=6 ; |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xy19 |
= |
7 +10 |
= 8,5 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X =(4+6+6+10+5+3+10+2+2+7+6+7)/12=5,66.
Y =(10+12+17+14+11+10+19+9+8+19+12+14)/12=12,9. Далее решаем по формулам.
|
|
|
|
|
N |
|
|
|
|
|
|
η yx |
= |
|
å fx (yx - y)2 |
||||
|
|
|
i=1 |
|
|
|
|||
|
|
|
N |
|
|
|
|||
|
|
|
|
|
å f y (yi - y)2 |
||||
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|||||
|
|
2× (8,5-12,9)2 + (10-12,9)2 + (10-12,9)2 + |
|||||||
|
|
(11-12,9)2 + 3× (13,6 -12,9)2 + 2× (16,5-12,9)2 |
|||||||
η yx = |
+ 2× (16,5 -12,9)2 |
|
|
||||||
(8 -12,9)2 + (9 -12,9)2 + 2× (10-12,9)2 + |
|||||||||
|
|
||||||||
|
|
(11-12,9)2 + 2× (12-12,9)2 + 2× (14-12,9)2 + |
|||||||
|
|
+ (17-12,9)2 + 2× (19-12,9)2 |
|||||||
|
|
ηyx |
=0,852 |
. |
|
|
|||
|
|
|
|
|
|||||
|
|
|
|
|
N |
|
|
|
|
|
|
ηхy |
= |
|
å f y (xy - x)2 |
||||
|
|
|
i=1 |
|
|
|
|||
|
|
|
N |
|
|
|
|||
|
|
|
|
|
å fx (xi - x)2 |
||||
|
|
|
|
|
i=1 |
|
|
|
|
271
(2- 5,6 6)2 + (2- 5,6 6)2 + 2× (3,5- 5,6 6)2 + (5- 5,6 6)2 + 2× (6- 5,6 6)2 + 2× (8,5- 5,6 6)2
+ (6- 5,6 6)2 + 2× (8,5- 5,6 6)2
η x y = 2× (2- 5,6 6)2 + (3- 5,6 6)2 + (4- 5,6 6)2 + (5- 5,6 6)2 + 3× (6- 5,6 6)2 + 2× (7- 5,6 6)2 + 2× (1 0- 5,6 6)2
ηxy =0,938 .
Проверим статистическую значимость коэффициентов
ηyx =0,852 |
|
|
|
|
ηxy |
=0,938 для df=12-2=10. |
|||||||||
ì |
3,1 6 9 д л яα = 0 ,0 |
|
|
|
|||||||||||
ï |
2 ,2 2 д8 л я |
α = 0 ,0 |
|
|
|||||||||||
tт а =б í |
|
|
|||||||||||||
ï |
1 ,8 |
1 д2 л я |
|
α = 0,1 |
|
|
|
||||||||
î |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
tэмп . = |
|
η |
|
× |
|
N -2 |
|
|
|
||||||
|
|
|
|
|
|||||||||||
|
|
|
1-η2 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
tэмп . = 0,852 × |
|
|
12 -2 |
|
|
tэмп . =5,146 |
|||||||||
1-0,852 |
2 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
Поскольку tэмп . >tтаб , постольку корреляционное отношение |
|||||||||||||||
между переменными Y и Х статистически значимо. |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
tэмп . = 0,938 × |
|
|
12 -2 |
|
|
tэмп . =8,577 |
|||||||||
1-0,938 |
2 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
Поскольку tэмп . >tтаб , постольку корреляционное отношение между переменными Х и Y статистически значимо.
Ответ: 1) Существует сильная, положительная, статистически значимая корреляционная связь между переменными Х и Y (корреляционное отношение Пирсона составило 0,938), с ошибкой в пределах 1%. То есть можно утверждать, что срок лишения свободы усиливает злокачественную агрессивность.
2. Существует сильная, положительная, статистически значимая корреляционная связь между переменными Y и Х (корреляционное отношение Пирсона составило 0,852), с ошибкой
272
в пределах 1%. То есть можно утверждать и обратное, что уровень злокачественной агрессивности увеличивает срок лишения свободы.
Если для переменных рассчитать линейный коэффициент корреляции, то он составит 0,79, а уравнение регрессии будет выглядеть следующим образом: y =6,6 +1,114 x . При этом параметры уравнения статистически значимы. Если поменять местами переменные, то результат будет: x =−1,692 +0,565 y при линейном
коэффициенте корреляции Пирсона 0,79.
§8. Множественный регрессионный (multiple regression) и корреляционный (correlation) анализ юридических процессов
Очевидно, что многообразие и сложность юридических процессов можно лишь крайне грубо и приблизительно сводить к простейшим парным линейным и нелинейным моделям. На самом деле на зависимую переменную оказывает воздействие значительное количество факторов, величина вклада которых различна и непостоянна во времени. В итоге колебания зависимой переменной – это равнодействующая всего многообразия управляющих ею сил, часть из которых мы сможем выделить и более или менее точно измерить величину их влияния. По сути, во множественном регрессионно-корреляционном анализе, мы строим модели вида:
y = f (x1, x2 ... xn ) . Они могут быть как линейными, так и нелинейными, что зависит от правила ( f ), связывающего левую и правую части уравнения. Знак над зависимой переменной указывает, что это лишь оценка зависимой переменной.
Используются |
|
|
линейные, |
степенные, |
экспоненциальные, |
|||||
гиперболические: |
|
|
|
|
|
|
||||
1) |
y = a +b1x1 +b2 x2 +... +bk xk |
+ε ; |
|
|||||||
2) |
|
b1 |
|
b2 |
|
bk |
+ ε ; |
|
|
|
y = a × x1 |
× x2 |
×...× xk |
|
|
||||||
3) |
y =e |
1 1 |
2 2 |
k |
k |
ε . |
|
|
||
|
|
a+b |
×x |
+b |
×x |
+... +b |
×x |
+ |
|
|
273
4) |
|
|
|
1 |
и другие модели множественной |
|
y = |
|
|
|
|||
a +b1 |
× x1 |
+b2 × x2 +... +bk × xk +ε |
||||
|
|
|
регрессии.
Оценки параметров соответствующих уравнений обычно находят с помощью метода наименьших квадратов (МНК)5. Как для линейных, так и для нелинейных уравнений, сводимых к линейным, строим систему нормальных уравнений:
ìå ïïíå ïïîå
y = n a+ b1å x1 + b2 å x2 + ...+bk å xk
y x1 = aå x1 + b1å x12 + b2 å x1x2 + ...+ bk å xk x1 y xk = aå xk + b1å x1xk + b2 å x2xk + ...bk å xk 2
Решение системы вручную можно выполнить с помощью метода определителей:
Проведение множественного регрессионного и корреляционного анализа гораздо более сложная задача, нежели исследование простых парных моделей. Здесь мы сталкиваемся с проблемами мультиколлинеарности и гетероскедастичности
(проблема гетероскедастичности возможна и для моделей парной регрессии).
Мультиколлинеарность – коррелированность между независимыми переменными, которая может приводить к несостоятельным оценкам. Особенно ярко проблема мультиколлинеарности проявляет себя при построении регрессий на основе временных рядов, то есть когда зависимая и независимые переменные взяты во временном развитии, например, за ряд лет, кварталов, месяцев и т.д. Важно понимать, что коррелированность между независимыми переменными не имеет отношения к точности спецификации модели. Например, коррелированные между собой независимые переменные вполне могут выступать в качестве объясняющих для исследуемой зависимой переменной, но если их использовать вместе может возникнуть проблема мультиколлинеарности, что приведет к неточным оценкам связи
5 Есть и другие методы параметризации, которые уместно использовать в определенных случаях.
274
между переменными, в частности, сила связи в модели может быть завышена. Две факторных переменных явно коллинеарны, находятся между собой в линейной зависимости, если парный коэффициент корреляции между ними не менее 0,7.
Для оценки величины мультиколлинеарности удобно использовать определитель матрицы парных коэффициентов корреляции. Если факторные переменные совершенно не коррелированны, то определитель матрицы парных коэффициентов корреляции будет равен единице. Покажем это на примере матрицы для трех независимых переменных:
|
|
|
|
|
rx1x1rx1x2rx1x3 |
|
|
100 |
|
|
|
|
|
|
|
|
|
||||
D |
|
R |
|
= |
r |
r er |
|
= |
01t0 |
=1. |
|
|
|||||||||
|
|
|
|
|
x2x1 |
x2x2x2x3 |
|
|
|
|
|
|
|
|
|
rx3x1rx2x3rx3x3 |
|
|
001 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Напротив, при функциональной связи между факторными переменными определитель матрицы будет равен нулю.
275
