Пример. Обнаружен труп и произведено измерение его температуры х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.
функции выступает показатель степени и т.п. Забегая вперед, важно отметить, что коэффициент регрессии в степенной функции одновременно выступает константой эластичности (в других видах функций коэффициент эластичности не является постоянной величиной). Регрессионный анализ или просто регрессия по существу означает подбор уравнения для каких-то количественных эмпирических данных. С точки зрения геометрии регрессия означает подбор линии (линейная регрессия) или кривой (нелинейная регрессия), которая наилучшим образом приближается к эмпирическим точкам. С точки зрения алгебры регрессия – это аппроксимирующее эмпирические данные линейное или нелинейное уравнение с конкретными параметрами. С точки зрения теории вероятностей и математической статистики – это аппроксимирующее эмпирические данные уравнение, содержащее случайный член (возмущение, шум), отражающий влияние: 1) неучтенных в модели факторов; 2) ошибок спецификации модели, например, неудачный подбор аппроксимирующей функции; 3) ошибок измерения. По этой причине параметры регрессионного уравнения всегда только оценки реальных параметров, надежность которых, во-первых, нужно проверять, а, во-вторых, всегда помнить о том, что связь носит вероятностно-статистический характер. Спецификация регрессионной модели – это формулирование гипотезы о связи между конкретными переменными, выбор зависимой и независимых (независимой) переменных, подбор правила (параметров линейного или нелинейного уравнения), связывающего переменные.
Регрессионный анализ обычно дополняется корреляционным, но даже коэффициент корреляции близкий к единице не гарантирует того, что одна переменная непременно детерминирует другую. То есть реально объясняет её поведение. Например, существует сильная корреляционная связь между умышленными убийствами и суицидами в нашей стране. Однако это вовсе не означает, что убийства объясняют суициды, или наоборот. Просто обе переменные характеризуют социально-патологическое поведение, которое детерминируется во многом схожим
причинным комплексом. Полезным из этой связи будет то, что по одной переменной можно будет прогнозировать другую, но не более того.
Параметрический корреляционный и регрессионный анализ по своей сути развивает функциональный анализ в вероятностностатистическом направлении, поскольку функциональная связь может рассматриваться как частный случай корреляционной с коэффициентом корреляции равным единице (положительная связь) или минус единице (отрицательная связь). Когда мы изучали функции, то игнорировали корреляцию, полагая её равной единице, и работали только с соответствующими линейными или нелинейными уравнениями. В параметрическом корреляционном и регрессионном анализе решаются две задачи: 1) находится коэффициент корреляции, показывающий силу связи между переменными (обычно рассчитывают так называемый линейный коэффициент корреляции Пирсона); 2) выводится регрессионное уравнение с соответствующими параметрами (константами а, b, c и т.д.), линейное или нелинейное, показывающее как объясняющая переменная (переменные) влияет на объясняемую. Кроме того, проводится работа, нацеленная на проверку качества коэффициента корреляции, параметров уравнения и всего уравнения в целом. Параметрическим регрессионный анализ называется потому, что находятся параметры уравнения, то есть здесь нужно найти неизменные константы, входящие в правую часть уравнения наряду с объясняющими (управляющими) переменными. Например, в простом линейном регрессионном уравнении вида у =а+bx, где - у оценки по переменной у, х – значения переменной икс, а (параметр) – свободный член (сдвиг, точка пересечения линии регрессионного уравнения с осью ординат), оценивающая величину игрек при икс равном нулю (иногда не интерпретируется, так как не имеет криминологического смысла), например, в случае, когда при параметре а стоит знак минус (-а). Важно заметить, что хотя в таком случае параметр а не интерпретируется, но интерпретируется знак минус. Если а<0 (очевидно, что минус а всегда меньше нуля), то изменения зависимой переменной (Δy) происходят быстрей, чем
изменения независимой переменной (Δ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.
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 |
|
|
|
А = 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 |
Способ №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) .
Хц |
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.
Даже получив значения параметров уравнения с помощью метода наименьших квадратов, не следует обольщаться, ибо мы имеем дело лишь с оценками реальных параметров уравнения, характеризующими соответствующие генеральные совокупности.
Обычно а ¹ α и 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%) – высоко значимый;