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

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

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

курсором, и они отобразятся в формуле автоматически. Фигурные скобки отобразятся в командной строке (строке формул) после выполнения 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

y = a +b × xc +u

Кубические функции неплохо моделируют развитие эпидемиологических и демографических процессов, в определенной мере связанных с криминологическими процессами, а кроме того сами криминологические процессы могут развиваться по подобному сценарию. Например, производственная функция определенного преступного предприятия, зависимость выпуска продукции от переменного ресурса, потребляемого за период времени может развиваться по кубическому правилу: 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