Ольков_С_Г_Аналитическая юриспруденция
.pdf
курсором, и они отобразятся в формуле автоматически. Фигурные скобки отобразятся в командной строке (строке формул) после выполнения 4-го действия.
4) Встаем курсором в первую ячейку столбца Х и набираем комбинацию Ctrl+Shift+Enter (одновременно нажимаем три клавиши). Если нажать просто Enter, то появится только одно значение – величина параметра а.
Таблица. Выполнение расчетов в Excel.
|
А |
В |
С |
D |
|
|
|
|
|
1 |
|
А |
|
В |
|
|
|
|
|
2 |
10 |
20339 |
|
18365 |
|
|
|
|
|
3 |
|
4231828 |
|
|
|
20339 |
1 |
|
38686988 |
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
5 |
Х |
|
Параметр |
|
|
-1018,04 |
«а» |
|
|
|
|
|
|
|
6 |
|
|
Параметр |
|
|
|
1,403481 |
«b» |
|
|
|
|
|
|
Не стоит говорить, что система линейных уравнений любой размерности решается с той же степенью легкости.
Среди функций нелинейных по объясняющим переменным для юридических исследований достаточно большое значение имеют квадратическая (парабола второго порядка), кубическая и гиперболическая функции. Так, квадратическая функция используется для описания зависимости уровня свободы от уровня государственного принуждения или общественного насилия. То есть представляет двумерную модель политических режимов. При вычислении коэффициента локализации удобно проводить параметризацию кривой Лоренца с помощью квадратической функции, а дальше находить разницу между этой квадратической функцией и линейной функцией, отражающей абсолютное равенство. Экономические процессы тесно связанные с криминологическими также иногда имеют вид параболы второго порядка.
186
Кубические функции неплохо моделируют развитие эпидемиологических и демографических процессов, в определенной мере связанных с криминологическими процессами, а кроме того сами криминологические процессы могут развиваться по подобному сценарию. Например, производственная функция определенного преступного предприятия, зависимость выпуска продукции от переменного ресурса, потребляемого за период времени может развиваться по кубическому правилу: 1) выработка растет с большой скоростью, (увеличивающаяся отдача); 2) выработка растет меньшими темпами (уменьшающаяся отдача); 3) отрицательная отдача, когда выработка начинает уменьшаться с увеличением затрат. Производная от этой функции, показывающая скорость изменения кубической функции опять же будет квадратической функцией.
Гиперболические функции неплохо моделируют реальное течение юридических процессов. Например, снижение уровня какого-либо вида преступности от угнетающего его фактора даст скорее гиперболическую, чем линейную зависимость.
Модели нелинейные по оцениваемым параметрам а, b, c и т.д. (модели второго класса), как мы выяснили, делятся на два подкласса: внутренне линейные (1) и нелинейные (2). Модели первого подкласса (внутренне линейные) сводятся к линейному виду и решаются с помощью МНК. Модели же второго класса решаются итерационными способами. Например, модель, содержащая степенную функцию: y = a × xb ×u приводится к линейному виду путем логарифмирования: ln y = ln a +ln x +ln u , в то время как модель: внутренне нелинейная20, и решать её придется итерационным способом.
Корреляция для нелинейных регрессионных моделей.
20 Подробнее: Эконометрика: Учебник/Под ред. И.И.Елисеевой. – М.: Финансы и статистика, 2001. С.70.
187
Нелинейный регрессионный анализ, как и линейный всегда дополняется корреляционным анализом, проверкой надежности подобранного уравнения, коэффициента корреляции и параметров уравнения, а также находим среднюю ошибку аппроксимации, с помощью которой судим об эффективности подобранного уравнения для характеристики эмпирических данных.
В случае парной линейной регрессии, как мы знаем, используется линейный коэффициент корреляции Пирсона. Для нелинейных же моделей применяется индекс корреляции, вычисляемый по формуле:
|
|
|
|
|
|
|
|
2 |
|
|
|
Dост |
|
|
|
|
|
||
R = 1− |
|
или R = 1 |
− |
å( yi − yi ) |
|
, где R - индекс корреляции, |
|||
|
|
|
|
||||||
Dy |
|
å( y − y) |
2 |
||||||
|
|
|
|
|
|
|
|
||
Dост – остаточная дисперсия, Dy – общая дисперсия управляемой переменной. Обратите внимание, что во второй по порядку формуле в числителе и знаменателе стоят суммы квадратов не приведенные на число степеней свободы (то есть не дисперсии, а суммы квадратов отклонений). В числителе стоит остаточная сумма квадратов отклонений, а в знаменателе общая сумма квадратов отклонений. Напомним, что
|
1 |
|
N |
|
2 |
|
|
1 |
|
N |
2 |
|
|
Dост = |
|
|
|
×å( yi - yi ) |
|
, |
Dy = |
|
|
×å( yi - yi ) |
|
. |
|
N - 2 |
|
N -1 |
|
||||||||||
|
i=1 |
|
|
|
|
i=1 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||
N |
|
2 |
- |
остаточная сумма |
квадратов (если её привести на |
||||||||
å( yi - yi ) |
|
||||||||||||
i=1
одну степень свободы, то как раз и получится остаточная дисперсия. Число степеней свободы для остаточной дисперсии зависит от числа независимых переменных. Для случая парной регрессии число степеней свободы остаточной дисперсии равно N-
N
2, поскольку имеется одна независимая переменная); å( yi - yi )2 -
i=1
общая сумма квадратов отклонений (число степеней свободы для нее N-1), которая складывается из факторной суммы квадратов или суммы квадратов, объясненных регрессионным уравнением:
N 2 (число степеней свободы для факторной (объясненной)
å( yi - yi )
i=1
188
суммы квадратов определяется по числу независимых переменных
– для парной регрессии df=1) и остаточной суммы квадратов отклонений:
N |
2 |
N |
|
2 |
N |
|
2 |
. |
å( yi − yi ) |
|
= å(yi − y) |
|
+ å( yi − yi ) |
|
|||
i=1 |
|
i=1 |
|
|
i=1 |
|
|
|
Индекс корреляции изменяется в пределах от нуля до единицы. Чем ближе индекс корреляции к единице, тем более точно кривая аппроксимирует эмпирические значения. Надежность индекса корреляции проверяется способом аналогичным с проверкой надежности линейного коэффициента корреляции: 1) сначала находим ошибку индекса корреляции по формуле:
mR = |
1− R2 |
; |
2) |
определяем |
фактическое |
значение |
t-критерия |
||||||
|
N −2 |
||||||||||||
|
|
|
|
|
|
|
|
R |
|
|
|
|
|
Стьюдента |
по |
формуле: |
tR = |
|
|
|
× |
N - 2 |
; 3) |
сравниваем |
|||
|
|
|
|
||||||||||
|
- R2 |
||||||||||||
1 |
|
|
|
||||||||||
фактическое значение t-критерия Стьюдента с табличным.
Кроме индекса корреляции для нелинейных регрессионных моделей используется индекс детерминации R2 и F-критерий
Фишера: Fэмп |
= |
|
|
R2 |
|
× |
N - m -1 |
, |
где m – |
число |
параметров |
при |
1 |
- R |
2 |
m |
|||||||||
|
|
|
|
|
|
|
|
|
||||
объясняющих |
|
|
|
переменных |
(число |
степеней |
свободы |
для |
||||
объясненной суммы квадратов); (N-m-1) – число степеней свободы для остаточной суммы квадратов.
На практике часто возникает вопрос о выборе между линейной и нелинейной моделью регрессии. Ответить на него помогает индекс детерминации, который сравнивается с
коэффициентом |
детерминации. |
Вычисляется |
оценка |
|||||||||||
существенности |
различия |
с помощью |
t-критерия |
Стьюдента: |
||||||||||
tR2 −r 2 |
= |
Rxy |
2 −rxy |
2 |
, |
где |
|
|
ошибка |
разности: |
||||
m 2 |
−r |
2 |
|
|
|
|||||||||
|
|
|
R |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||||
mR2 −r 2 |
= 2× |
|
(R2 |
- r 2 ) - (R2 - r 2 )2 ×(2×(R2 + r2 ) |
|
. Как и обычно, если фактическое |
||||||||
|
|
|
|
|
|
N |
|
|
||||||
189
значение t-критерия Стьюдента превышает табличное, то различия между индексом корреляции и коэффициентом корреляции существенны, и необходимо использовать нелинейную модель. Напротив, если фактическое значение t-критерия Стьюдента меньше табличного, то не стоит «ломать копий» и пользоваться нелинейной моделью. В данном случае нужно остановиться на линейной модели. На практике, если tфакт<2, то следует выбрать линейную модель. Кроме того, если разница между индексом корреляции и коэффициентом корреляции не превышает 0,1 (Rxy 2 −rxy 2 ) < 0,1 , то можно сразу отдать предпочтение линейной модели.
Для оценки качества модели нелинейной регрессии, как и для линейной, наряду с другими показателями используется средняя ошибка аппроксимации:
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
2 |
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
||||||||
|
|
= |
1 |
× å |
( yi - yi |
|
|
×100% |
|
|
|
100 |
|
å(yi - yi ) |
|
|
|
|
A |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
N |
yi |
|
или A = |
× |
. |
||||||||||||
|
|
|
i 1 |
|
|
|
i =1 |
|
||||||||||
|
|
|
|
= |
|
|
|
|
|
|
|
y |
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Проведем параметризацию некоторых нелинейных функций, оцениваемых по параметрам.
Степенная функция вида y = a × xb ×u приводится к следующей системе нормальных уравнений:
ì |
|
N |
N |
|
ï N × lna |
+ b× å lnx |
= å lny |
|
|
ï |
|
i=1 |
i=1 |
|
í |
|
|
|
|
ïïlna× åN |
lnx + b× åN |
(lnx)2 = åN |
lny× lnx |
|
î |
i=1 |
i=1 |
i=1 |
|
Коэффициент регрессии (b) и оценку свободного члена (а) находим, решив систему нормальных уравнений и выполнив операцию потенцирования. Важно отметить, что логарифм уместно
190
брать по любому основанию (это будет отражаться только на процедуре потенцирования, но не на результате).
Дано: таблица первичных выборочных данных, где Х – число пробелов в нормативно-правовых актах, регламентирующих данную деятельность, за фиксированный период времени на 7 объектах, Y – число нарушений данных нормативно-правовых актов на указанных объектах, связанных с выявленными пробелами.
Объект |
Y |
X |
log Y |
log X |
logY∙logX |
( yi − y)2 |
|
ы |
|
|
|
|
|
|
|
№ п/п |
|
|
|
|
|
|
|
1 |
746 |
53 |
2,87 |
1,72 |
4,94 |
2809 |
|
2 |
852 |
67 |
2,93 |
1,83 |
5,36 |
||
2809 |
|||||||
3 |
710 |
48 |
2,85 |
1,68 |
4,78 |
||
7921 |
|||||||
4 |
811 |
61 |
2,9 |
1,79 |
5,19 |
||
144 |
|||||||
5 |
914 |
72 |
2,96 |
1,86 |
5,51 |
||
13225 |
|||||||
6 |
1043 |
81 |
3,01 |
1,91 |
5,75 |
||
59536 |
|||||||
7 |
517 |
43 |
2,71 |
1,63 |
4,42 |
||
79524 |
|||||||
|
|
|
|
|
|
||
Сумма |
5593 |
425 |
20,23 |
12,42 |
35,95 |
165968 |
|
Средне |
799 |
60,7 |
2,89 |
1,77 |
|
|
|
е |
|
|
|
|
|
|
|
D |
27661, |
185,5 |
0,00923 |
0,0102 |
|
|
|
|
3 |
7 |
3 |
9 |
|
|
|
D |
166,31 |
13,6 |
0,09609 |
0,1014 |
|
|
|
|
|
|
|
|
|
||
|
7 |
|
|
65 |
|
|
191
Продолжение таблицы |
|
|
|
|
|
|
|
|
|
|
|
(logХ)2 |
yi |
yi − yi |
( yi − yi ) |
|
Ai |
|
|
|
|
2 |
|
2,96 |
701 |
45 |
2025 |
(45/746)∙100% |
|
|
|
|
=6,03% |
3,35 |
863 |
-11 |
121 |
(11/852)∙100% |
|
|
|
|
=1,29 |
2,82 |
642 |
68 |
4624 |
(68/710)∙100% |
|
|
|
|
=9,58 |
3,20 |
794 |
17 |
289 |
(17/811)∙100% |
|
|
|
|
=2,096 |
3,46 |
920 |
-6 |
36 |
(6/914)∙100% |
|
|
|
|
=0,656 |
3,65 |
1021 |
22 |
484 |
(22/1043)∙100% |
|
|
|
|
=2,1 |
2,66 |
583 |
-66 |
4356 |
(66/517)∙100% |
|
|
|
|
=12,77 |
Итого |
5524 |
69 |
11935 |
|
22,10 |
|
|
|
|
Среднее |
|
|
|
4,92 |
D |
|
|
11935/5= |
|
|
|
|
2387 |
|
192
Для того, чтобы заполнить столбец с оценками игрек, необходимо решить систему нормальных уравнений:
ì |
|
|
N |
|
|
N |
|
|
ï N × lna |
+ b× å lnx |
= å lny |
|
|||||
ï |
|
i=1 |
|
|
i=1 |
|
|
|
í |
|
lnx + b× åN |
( lnx)2 = åN |
|
||||
ïïlna× åN |
lny× lnx |
|||||||
î |
i=1 |
|
|
i=1 |
|
|
i=1 |
|
ì |
7× l na+ b×1 |
2,4 |
2= |
2 |
0,2 |
3 |
|
|
í |
l na×1 |
2,4 2+ |
b× |
2 2,1 |
0= 3 |
5,9 |
|
|
î |
|
|||||||
Составим матрицы:
A= |
é |
7 |
1 |
ù2 |
В= |
é |
2 |
ù |
ê |
1 |
22 |
ú |
ê |
3 |
ú |
||
|
ë |
,û2 |
|
ë |
û |
Находим матрицу Х с параметрами уравнения:
Х= |
é |
1,3 |
1ù. |
|
|
ê |
0 |
, |
ú |
|
ë |
û |
||
Получилось оценочное линейное уравнение регрессии:
193
Y =1,3165 +0,887 × X .
Выполнив потенцирование, получим окончательное оценочное уравнение по степенной функции:
|
=10 |
1,3165 |
× x |
0,887 |
= 20,725 × x |
0,887 |
|
y |
|
|
|
|
|||
|
|
× x |
0,887 |
. |
|
||
y |
=20,725 |
|
|
||||
Подставляя значения независимой переменной (х) в полученное уравнение степенной функции, заполним столбец с
оценками зависимой переменной (игрек): 1) |
y1 = 20,725 ×53 |
0,887 |
= 701 |
… |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
= |
|
1 |
|
|
|
|
|
yi |
- yi |
|
×100% = 4,92 |
|
|
|
|
|
|
|
|
|
|
|||||||||
A |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
N |
|
y |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
i=1 |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
100 |
|
|
å(yi - yi ) |
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
11935 |
|
|
|
|
|
|
|||||||||||||||||
A = |
× |
|
|
= |
× |
= 5,168% |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
y |
|
|
|
|
|
N |
|
|
|
|
799 |
|
|
7 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Найдем индекс корреляции Rxy : |
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Dост |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Rxy |
= |
|
1− |
= |
|
|
1− |
|
å( уi − yi ) |
|
|
|
|
|
|
|
|
|||||||||||||||
|
Dy |
|
|
|
å( yi − y)2 |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Rxy |
= 1- |
11935 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
165968 |
=0,96. |
|
|
Проверим индекс |
корреляции |
на |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
статистическую значимость с помощью t-критерия Стьюдента: 1) сначала найдем индекс детерминации – возведем в квадрат индекс корреляции: R2=0,92.
|
2) Находим ошибку индекса |
корреляции по |
формуле: |
||
|
|
|
|
|
|
|
1- R2 |
. В нашем случае: mR = |
|
. Отсюда |
получим: |
mR = |
1 -0,92 |
||||
|
N -2 |
|
7 -2 |
|
|
mR = 0,126 .
3) Определяем фактическое значение t-критерия Стьюдента по
формуле: |
tR = |
|
|
R |
|
× |
|
или t R= |
R |
. |
|
|
|
N - 2 |
|||||||
|
|
|
|
|
||||||
|
- R2 |
|||||||||
|
1 |
|
|
|
|
mR |
||||
194
tR = |
0,96 |
|
|
|
|
|
|
|
|
|
|
|||
|
× 7 |
- 2 = 7,58 |
или с использованием ошибки индекса |
|||||||||||
|
|
|
|
|||||||||||
1 |
- 0,92 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||
корреляции:t R= |
R |
= |
0,96 |
= 7,61 . |
||||||||||
mR |
|
0,126 |
||||||||||||
3) сравниваем фактическое значение t-критерия Стьюдента с табличным при числе степеней свободы в остаточной дисперсии (в данном случае парной регрессии N-2=7-2=5):
ì |
2,5 7 д л αя= 0 , 0 |
tт а=бí |
6 , 8 5д9л αя= 0 , 0 |
î |
Как видно эмпирическое значение t-критерия составило 7,61, что выше критического при любом уровне значимости.
Проверим целесообразно ли заменять линейную регрессионную модель моделью степенной. Ответить на данный вопрос помогает индекс детерминации, который сравнивается с коэффициентом детерминации. Вычисляется оценка существенности различия с помощью t-критерия Стьюдента:
|
|
R |
2 |
−r |
2 |
|
|
|
|
|
|
|
|||
tR2 −r 2 |
= |
|
xy |
|
|
|
xy |
|
, |
где |
ошибка |
|
разности: |
||
|
m |
2 |
−r |
2 |
|
|
|||||||||
|
|
|
R |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||||
mR 2 −r 2 |
= 2 × |
|
(R2 - r2 ) - (R2 - r2 )2 × (2 - (R2 + r2 )) |
. |
Для |
нашего |
примера |
||||||||
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
||
коэффициент детерминации, полученный из линейного коэффициента корреляции Пирсона составит 0,931 (r2=0,931).
Отсюда: tR 2 −r 2 = |
0,92 −0,931 |
= -0,139 . |
|||||||||
|
0,079 |
|
|||||||||
|
|
|
|
|
(0,9 2- 0,9 3)1- (0,9 2- 0,9 3)12 × |
|
|
|
|
||
m 2 |
|
|
= 2× |
|
(2- (0,9 2+ 0,9 321)) |
|
= |
0,0 7 |
. |
||
|
2 |
|
7 |
|
|
|
|||||
R |
− r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Видно, что эмпирическое t меньше двух (всего -0,139). Следовательно, не имеет смысла строить степенную функцию.
195
