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

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

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

Пример. Обнаружен труп и произведено измерение его температуры х0 = 310С в момент обнаружения. Спустя час х = 290С. Допустим, что в момент смерти температура тела в данном месте х0 = 370С и температура окружающей среды была 210С.

k =ln

31

−21 =ln1,25

=0,22314

 

29

−21

 

 

 

 

t = −ln 37 −21

=

1

 

ln1,6 = −2,10630

0,22314

 

31 −21

 

 

Таким образом, преступление было совершено за 2 часа 6 минут до момента обнаружения тела.

Как видно расчеты весьма простые потому, что температура окружающей среды (а) принимается постоянной. Такое допущение иногда вполне оправдано, поскольку резкие температурные броски окружающей среды не слишком частое явление. Однако, незначительные колебания (реже – значительные) вполне возможны. Тогда закон охлаждения трупа примет более сложную

форму: dx

+kx = ka(t) , где а(t) – температура воздуха в момент

dt

 

времени t. Выше приведенное уравнение носит название линейного неоднородного дифференциального уравнения первого порядка15.

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

Термин параметр берет свое начало от греческого parametren - отмеривающий, измеряющий что-либо». Напомню, что в

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

15 См. Амелькин В.В., Садовский А.П. Математические модели и дифференциальные уравнения. – Мн.: Высш.школа, 1982. С.20-25.

156

функции выступает показатель степени и т.п. Забегая вперед, важно отметить, что коэффициент регрессии в степенной функции одновременно выступает константой эластичности (в других видах функций коэффициент эластичности не является постоянной величиной). Регрессионный анализ или просто регрессия по существу означает подбор уравнения для каких-то количественных эмпирических данных. С точки зрения геометрии регрессия означает подбор линии (линейная регрессия) или кривой (нелинейная регрессия), которая наилучшим образом приближается к эмпирическим точкам. С точки зрения алгебры регрессия – это аппроксимирующее эмпирические данные линейное или нелинейное уравнение с конкретными параметрами. С точки зрения теории вероятностей и математической статистики – это аппроксимирующее эмпирические данные уравнение, содержащее случайный член (возмущение, шум), отражающий влияние: 1) неучтенных в модели факторов; 2) ошибок спецификации модели, например, неудачный подбор аппроксимирующей функции; 3) ошибок измерения. По этой причине параметры регрессионного уравнения всегда только оценки реальных параметров, надежность которых, во-первых, нужно проверять, а, во-вторых, всегда помнить о том, что связь носит вероятностно-статистический характер. Спецификация регрессионной модели – это формулирование гипотезы о связи между конкретными переменными, выбор зависимой и независимых (независимой) переменных, подбор правила (параметров линейного или нелинейного уравнения), связывающего переменные.

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

157

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

Параметрический корреляционный и регрессионный анализ по своей сути развивает функциональный анализ в вероятностностатистическом направлении, поскольку функциональная связь может рассматриваться как частный случай корреляционной с коэффициентом корреляции равным единице (положительная связь) или минус единице (отрицательная связь). Когда мы изучали функции, то игнорировали корреляцию, полагая её равной единице, и работали только с соответствующими линейными или нелинейными уравнениями. В параметрическом корреляционном и регрессионном анализе решаются две задачи: 1) находится коэффициент корреляции, показывающий силу связи между переменными (обычно рассчитывают так называемый линейный коэффициент корреляции Пирсона); 2) выводится регрессионное уравнение с соответствующими параметрами (константами а, b, c и т.д.), линейное или нелинейное, показывающее как объясняющая переменная (переменные) влияет на объясняемую. Кроме того, проводится работа, нацеленная на проверку качества коэффициента корреляции, параметров уравнения и всего уравнения в целом. Параметрическим регрессионный анализ называется потому, что находятся параметры уравнения, то есть здесь нужно найти неизменные константы, входящие в правую часть уравнения наряду с объясняющими (управляющими) переменными. Например, в простом линейном регрессионном уравнении вида у =а+bx, где - у оценки по переменной у, х – значения переменной икс, а (параметр) – свободный член (сдвиг, точка пересечения линии регрессионного уравнения с осью ординат), оценивающая величину игрек при икс равном нулю (иногда не интерпретируется, так как не имеет криминологического смысла), например, в случае, когда при параметре а стоит знак минус (-а). Важно заметить, что хотя в таком случае параметр а не интерпретируется, но интерпретируется знак минус. Если а<0 (очевидно, что минус а всегда меньше нуля), то изменения зависимой переменной (Δy) происходят быстрей, чем

158

изменения независимой переменной (Δx): y> x, что объясняется величиной производной (параметр b), тангенсом угла наклона (линия уходит под ноль по оси ординат, но проходит через линию абсцисс за отметкой ноль в области положительных чисел). Неслучайно, если мы возьмем коэффициент вариации по зависимой переменной и коэффициент вариации по независимой переменной

и сравним их, то получим: Vy = σyy >Vx = σxx . Соответственно и наоборот, если а>0, то изменения зависимой переменной

происходят медленнее, чем независимой y< x,

Vy =

σy

<Vx =

σx

.

 

 

 

 

y

x

Говоря о свободном члене, следует знать, что он отражает

совокупную величину влияния неучтенных в данной модели факторов.

Параметр b – коэффициент регрессии (коэффициент пропорциональности, первая производная или скорость изменения функции), показывающий, на сколько в абсолютном выражении изменяется переменная игрек при изменении переменной икс на единицу измерения; b - тангенс угла наклона касательной (tg α) в любой точке функции при х→0. То есть

b = dy

= lim

y

= f ′(x) = y′ = y = tgα.

dx

x→0

x

 

По сути, соотносятся абсолютные приросты по переменным х

и у:

yx , где у=у1-у0, х=х1-х0.

Важно заметить, что на каждый параметр при независимой

переменной в параметрическом регрессионном анализе должно приходиться не менее 7 наблюдений (параметр а во внимание не принимается) То есть в парном регрессионном анализе должно быть не менее 7 значений по исследуемым переменным икс и игрек (7 – строк таблицы). Если проводится множественный регрессионный анализ, например, с тремя независимыми переменными, то наблюдений должно быть не менее 21.

Параметризацию юридических функций можно проводить различными способами: 1) метод контрольных точек; 2) метод

средних; 3) метод сумм; 4) метод кривых Пирсона; 5) метод 159

наименьших абсолютных отклонений; 6) метод наименьших квадратов и его модифицированные варианты (дисконтированный МНК, двух шаговый МНК, косвенный МНК и т.д.)16, а также другими методами.

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

Таблица №1. Первичные статистические данные об изменениях процессов икс и игрек, а также абсолютные приросты по обеим переменным

X

Y

X

Y

23

45

 

 

34

56

11

11

43

61

9

5

52

70

9

9

64

87

12

17

76

98

12

11

83

109

7

11

96

117

13

8

111

130

15

13

9

 

 

 

åyi

=773

Х =11

Y =10,625

i=1

X =64,66667; Y =85,88889

Способ №1. Учитывая тот факт, что в сути производная (скорость изменения функции) отражает отношение абсолютных приростов Y/ X её приближенное значение можно найти, как частное от деления среднего абсолютного прироста по переменной

игрек

на средний абсолютный прирост по переменной икс: b=

 

 

 

 

10 ,625

 

Y

=

=0,965909. Свободный член найдем по формуле:

 

 

11

X

16Статистические модели в экономике/А.Ф.Гришин, С.Ф.Котов-Дарти, В.Н.Ягунов.

–Ростов н/Д: «Феникс», 2005. С. 93.

160

a=Y −bX =85,88889-0,965909∙64,66667=23,42677 Следовательно, оценочное уравнение: yˆ =23,42677+0,965909х. Проверим его качество.

 

 

X

 

 

Y

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

23

 

 

45

 

 

 

 

 

45,64268

 

 

 

 

 

 

 

34

 

 

56

 

 

 

 

 

56,26768

 

 

 

 

 

 

 

43

 

 

61

 

 

 

 

 

64,96086

 

 

 

 

 

 

 

52

 

 

70

 

 

 

 

 

73,65404

 

 

 

 

 

 

 

64

 

 

87

 

 

 

 

 

85,24495

 

 

 

 

 

 

 

76

 

 

98

 

 

 

 

 

96,83585

 

 

 

 

 

 

 

83

 

 

109

 

 

 

 

 

103,5972

 

 

 

 

 

 

 

96

 

 

117

 

 

 

 

 

116,154

 

 

 

 

 

 

 

111

 

 

130

 

 

 

 

 

130,6427

 

 

 

 

 

 

 

 

 

 

 

773

 

 

 

 

 

773

 

 

 

 

 

 

 

 

 

Проверим точность расчетов.

 

 

Видно

 

что

9

9

 

 

 

 

 

åY = åY , что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

i =1

 

говорит о хорошем

качестве

аппроксимации.

Далее

найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

среднюю ошибку аппроксимации

 

=

1

×å

 

yi - yi

×100 .

 

 

 

 

А

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

- yi

 

×100

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

45

 

 

45,64

 

 

 

 

 

 

 

 

1,422

 

 

 

 

 

 

 

 

56

 

 

56,27

 

 

 

 

 

 

 

 

0,482

 

 

 

 

 

 

 

 

61

 

 

64,96

 

 

 

 

 

 

 

 

6,492

 

 

 

 

 

 

 

 

70

 

 

73,65

 

 

 

 

 

 

 

 

5,214

 

 

 

 

 

 

 

 

87

 

 

85,24

 

 

 

 

 

 

 

 

2,023

 

 

 

 

 

 

 

 

98

 

 

96,84

 

 

 

 

 

 

 

 

1,184

 

 

 

 

 

 

 

 

109

 

 

103,597

 

 

 

 

 

 

 

 

4,957

 

 

 

 

 

 

 

 

117

 

 

116,15

 

 

 

 

 

 

 

 

0,726

 

 

 

 

 

 

 

 

130

 

 

130,64

 

 

 

 

 

 

 

 

0,492

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑=22,992

 

 

 

161

А = 22,992 = 2,55% . Средняя ошибка аппроксимации не превышает

9

3%. В свою очередь средняя ошибка аппроксимации в пределах 5- 7% считается хорошей (8-10% - допустимой). Вообще говоря, в данном случае она не превысила бы и 2%, если нами не проводилось округление до тысячных по оценочным значениям.

Способ №2. Отличается от первого только способом нахождения среднего абсолютного прироста:

X

Y

Y/ X

23

45

 

34

56

1

43

61

0,555556

52

70

1

64

87

1,416667

76

98

0,916667

83

109

1,571429

96

117

0,615385

111

130

0,866667

 

 

Среднее=0,992796

а =Y −bX =85,88889-0,992796∙64,66667=21,68808 Следовательно, оценочное уравнение: yˆ =21,68808+0,992796х. Проверим его качество.

X

Y

Y

 

 

ˆ

23

45

44,51431

34

56

55,43506

43

61

64,37023

52

70

73,30539

64

87

85,21894

76

98

97,1325

83

109

104,0821

96

117

116,9884

111

130

131,8804

 

773

772,9273

162

Способ №3. Основной. С помощью метода наименьших квадратов.

b = Kx , y

Dx

 

1

N

 

 

 

 

 

1

 

N

Dx =

å(xi - x)2

или

Dx =

 

 

 

å(xi - x)2 (несмещенная оценка

 

N

 

N -1

 

 

i=1

 

 

 

 

i =1

дисперсии).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

N

 

 

 

 

 

 

 

 

 

 

Кx, y =

å(xi - x) ×( yi - y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

Y

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

23

 

 

45

 

 

 

43,9811

 

 

 

 

34

 

 

56

 

 

 

55,0438

 

 

 

 

43

 

 

61

 

 

 

64,0951

 

 

 

 

52

 

 

70

 

 

 

73,1464

 

 

 

 

64

 

 

87

 

 

 

85,2148

 

 

 

 

76

 

 

98

 

 

 

97,2832

 

 

 

 

83

 

 

109

 

 

104,3231

 

 

 

 

96

 

 

117

 

 

117,3972

 

 

 

 

111

 

130

 

 

132,4827

 

 

 

 

 

 

 

773

 

 

772,9674

 

Рис. №10.

Можно работать и с центрированными данными, заменив исходные переменные центрированными: хцi = (xi − x) и yцi = ( yi − y) .

163

Хц

Yц

 

 

-41,7

-40,9

 

 

-30,7

-29,9

 

 

-21,7

-24,9

 

 

-12,7

-15,9

 

 

-0,7

1,1

 

 

11,3

12,1

 

 

18,3

23,1

 

 

31,3

31,1

 

 

46,3

44,1

 

 

x =64,66667 y =85,88889

Рис. №11.

Как видно из уравнения коэффициент наклона не изменился. Изменился свободный член (стал равен нулю). Чтобы прогнозировать по данному уравнению, нужно восстанавливать свободный член: a=Y −bX =85,88889-1,005∙64,66667=20,8.

Даже получив значения параметров уравнения с помощью метода наименьших квадратов, не следует обольщаться, ибо мы имеем дело лишь с оценками реальных параметров уравнения, характеризующими соответствующие генеральные совокупности.

164

Обычно а ¹ α и b ¹ β , где α – значение свободного члена в генеральной совокупности, а β – значение коэффициента регрессии (тангенса угла наклона линии) в генеральной совокупности. Существует необходимость проверки статистической надежности оценок.

Оценивается статистическая значимость: 1) уравнения в целом; 2) параметров уравнения; 3) коэффициента корреляции.

Для оценки параметров уравнения, коэффициента корреляции и коэффициента детерминации выполняется

следующий алгоритм:

1) выдвигается статистическая гипотеза H0 – показатели статистически незначимы, то есть незначимо отличаются от нуля;

2) рассчитывается фактический t-критерий Стьюдента:

а) для коэффициента корреляции по формуле: tr = mrr , где tr –

это эмпирическое значение t-статистики для коэффициента корреляции, которое сравнивается с табличным значением t-критерия Стьюдента, r – коэффициент корреляции, mr –

случайная ошибка коэффициента корреляции: mr = 1− r2 .

N − 2

Если число наблюдений превышает 70 (N≥70), а распределение случайных величин Х и Y близко к нормальному, то практически достоверно, что истинное значение коэффициента корреляции полученное по выборочным данным располагается в интервале между [r-3σr] и [r+3σr], где σr - стандартное отклонение

выборочных корреляций, вычисляется по формуле: σ = 1−r2

r N .

P.S. Тестирование коэффициента детерминации R2 осуществляется с помощью специальных таблиц критических значений данного коэффициента при различных уровнях значимости (0,05 (5%) – значимый; 0,01 (1%) – высоко значимый;

165