Ольков_С_Г_Аналитическая юриспруденция
.pdf
0,001 (0,1%) – особо высоко значимый; 0,1 (10%) – средне значимый) 17. В верхней строке таблицы указывается число независимых переменных, а в левом крайнем столбце - количество наблюдений (N). Например, если мы работали с парной линейной регрессией, то выбираем критическое значение коэффициента детерминации в столбце, соответствующем одной переменной. Если у нас имелось 10 наблюдений, то критическое значение R2 составит 0,399 при α=0,05. Если имелось 3 независимых переменных (модель множественной регрессии) при 40 наблюдениях, то критическое значение R2 составит 0,193 при α=0,05. При любом числе переменных, когда число наблюдений превышает 60, критическое значение коэффициента детерминации вычисляется по специальной формуле и множителям, указанным в таблице критических значений:
Крит.знач. = |
множитель 1 |
× |
множитель 2 |
. |
|
N |
|
N 2 |
б) Для коэффициента регрессии по формуле:
это эмпирическое (фактическое) значение t-критерия Стьюдента для коэффициента регрессии, b – коэффициент
регрессии, mb – случайная ошибка коэффициента регрессии:
|
|
N |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ N − 2 |
|
|
|
|
|
|
|
|
Dостаточная |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
å(yi − yi ) |
|
|
|
|
mb = |
|
|
|
|||||||
mb = |
i=1 |
|
|
|
|
или что равнозначно: |
|
N |
, |
||||||||
|
|
|
|
|
|
|
|
å(xi − x)2 |
|||||||||
|
N |
|
2 |
|
|
|
|||||||||||
|
|
|
å(xi − x) |
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
σ y |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
или через стандартную ошибку оценки: mb = |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||
σx × N -1 |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
в) для свободного члена по формуле: ta = |
a |
, где tа – это t- |
|||||||||||||||
m |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
критерий Стьюдента для коэффициента сдвига (свободного
17 Сигел Э.Ф. Практическая бизнес-статистика: Пер. с англ. – М.: Издательский дом «Вильямс», 2002. С. 632-641.
166
члена), который также будет сравниваться с табличным, а
– коэффициент сдвига, mа – случайная ошибка свободного
|
|
N |
2 |
|
N |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
å(yi - yi ) |
|
å x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|||
члена: ma = |
|
|
|
или ma =σy × |
+ |
X |
||||||||||
|
|
× N |
|
|
N |
D ×(N -1) . |
||||||||||
|
|
i=1 |
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N - 2 |
|
N × å(xi |
- x)2 |
|
|
x |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
3)Находим табличное значение t-критерия Стьюдента, которое зависит от: 1) числа степеней свободы; 2) числа
определяемых по ним констант. В случае парной регрессии число степеней свободы df=N-2.
4) Сравниваются фактическое значения t-критерия Стьюдента, полученные по вышеприведенным формулам с табличными (теоретическими) значениями t-критерия Стьюдента. Если tтеор. < tфакт., то принимается статистическая гипотеза H1, утверждающая статистическую значимость соответствующего коэффициента, его значимое отличие от нуля. В противном случае, если tтаб. > tфакт., то принимается статистическая гипотеза H0, указывающая на незначимое отличие оценочного значения параметра от нуля.
Рассчитаем |
|
t-критерий Стьюдента |
для |
нашего |
уравнения: |
||||||||||
y =20 ,85 +1,005 x . При R2=0,992, r= |
|
|
|
=0,995. |
|
||||||||||
0,992 |
|
|
|||||||||||||
Проверим |
|
статистическую |
значимость коэффициента |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
1- r2 |
= |
|
=0,034; |
|
|
|
|
|
||||
корреляции: mr = |
|
1- 0,992 |
|
|
|
|
|
||||||||
|
|
|
|
N - 2 |
9 - 2 |
|
|
|
|
|
|
|
|||
tr = |
r |
= 0,995 |
= 29,26 |
. Находим табличное |
значение |
t-критерия |
|||||||||
mr |
|||||||||||||||
|
0,034 |
|
|
|
|||||||||||
Стьюдента при 7 степенях свободы (9-2) для уровня значимости α=0,01 (то есть вероятность ошибки составляет только 1%): tтабл.=3,49. Поскольку 3,49<29,26 примем статистическую гипотезу H1 о статистически значимом отличии коэффициента корреляции от нуля.
167
Проверим |
статистическую значимость коэффициента |
|||||||
регрессии: |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
N |
|
2 |
|
|
|
|
mb = |
å(yi − yi ) / N − 2 |
Dост |
|
|||||
i=1 |
|
|
= |
|
||||
|
N |
|
N |
|||||
|
|
|
å(xi |
− x)2 |
å(xi − x)2 . |
|||
|
|
|
i=1 |
|
|
|
i=1 |
|
Для оценки значимости уравнения в целом используется нижеследующий алгоритм:
1). Находится эмпирическое значение F-критерия Фишера по
формуле: F = |
Dфакт |
= |
r2 |
×(N - 2) , где F – это F-критерий Фишера для |
|
D |
|
1- r2 |
|
|
ост |
|
|
|
проверки статистической гипотезы H0: Dфакт = Dост . Чтобы доказать
альтернативную статистическую гипотезу H1 нужно, чтобы факторная дисперсия (объясненная регрессией) в несколько раз превышала остаточную дисперсию; Dфакт – факторная дисперсия (сумма квадратов отклонений зависимой переменной объясненная регрессией приведенная на число степеней свободы); Dост – остаточная дисперсия (остаточная сумма квадратов отклонений приведенная на число степеней свободы); r2 - коэффициент детерминации, характеризует долю дисперсии зависимой переменной (у) объясняемую регрессионным уравнением в общей дисперсии:
r2 = Dфакт
Dобщ
=æç Kx,y çè Dx
ö2 |
æ |
ö2 |
æ |
ö2 |
|
|
|
÷ |
= çb × |
σx |
÷ |
= ç |
ср xy - y × x ÷ |
где |
ср ху – среднее от |
|
|||||||
÷ |
ç |
÷ |
ç |
÷ |
|||
ø |
è |
σ y ø |
è |
σ x ×σ y ø |
|
|
|
произведения значений х помноженных на значения у. Соответственно коэффициент корреляции r находится по формуле:
|
|
|
|
Kx,y |
|
σx |
|
cр xy - y × x |
|
r = r |
2 |
= |
= b × |
= |
. |
||||
|
Dx |
σ y |
σx ×σ y |
||||||
168
Важно отметить, что величина 1-r2 показывает долю факторной дисперсии объясненную независимыми переменными, не включенными в данную объяснительную модель.
2)Находится табличное (теоретическое) значение F-критерия Фишера.
3)Проводится сравнение табличного значения со значением эмпирическим. При этом, если эмпирическое значение больше теоретического, то уравнение статистически значимо.
F-критерий Фишера, как и другие критерии находится для различных доверительных уровней 0,1; 0,05; 0,01; 0,001. Величина
вероятной ошибки зависит от величины доверительного уровня и оценивается как 1-р, где р - уровень значимости18. Например, для р=0,05 вероятность ошибки составляет 5%. Очевидно, чем меньше ошибка, тем выше значение F-критерия Фишера.
Учитывая тот факт, что для нахождения F-критерия Фишера и понимания сути регрессионно-корреляционного анализа ключевое значение имеет разложение общей суммы квадратов отклонений зависимой переменной от среднего (TSS – total sum of squares) на
объясненную (ESS – explain sum of squares) и необъясненную
(RSS-residual sum of squares) составляющие покажем это:
å( yi − yi ) |
|
= å(yi − y) |
|
+ å( yi − yi ) |
, где |
||||
N |
|
|
2 |
N |
|
2 |
N |
|
2 |
|
|
|
|
|
|||||
i=1 |
|
|
|
i=1 |
|
|
i=1 |
|
|
N |
|
|
|
- общая сумма квадратов отклонений; |
|||||
å( yi − yi ) |
2 |
||||||||
i=1 |
|
|
|
|
|
|
|
|
|
N |
|
2 |
- факторная сумма квадратов отклонений или сумма |
||||||
å( yi − yi ) |
|
|
|||||||
i=1 |
|
|
|
|
|
|
|
|
|
квадратов отклонений, объясненная регрессией; |
|||||||||
N |
|
2 |
- остаточная сумма квадратов отклонений или сумма |
||||||
å( yi − yi ) |
|
|
|||||||
i=1
квадратов не объясненная регрессией.
18 Также обозначают буквой α.
169
Поиск сумм квадратов отклонений имеет важное значение, поскольку обеспечивает: 1) получение коэффициента корреляции, например, коэффициент корреляции для нелинейной регрессии рассчитывается по общему принципу:
|
|
|
|
|
2 |
|
|
|
|
Dост |
|
|
|
||
r = |
|
|
å( уi − yi ) |
|
, где |
r - коэффициент корреляции, |
|
1− |
Dy |
= |
1− å( yi − y)2 |
||||
Dост – остаточная дисперсия, Dy – общая дисперсия управляемой переменной; 2) проверку значимости: а) уравнения в целом (находим F-критерий Фишера); б) параметров уравнения; в) коэффициента корреляции.
Общая сумма квадратов отклонений управляемой (объясняемой)
переменной игрек (у) |
N |
по |
от среднего значения игрек å( yi - yi )2 |
||
|
i=1 |
|
существу и определяет величину силы связи между переменными
икс и игрек. |
В |
случае, |
когда |
остаточная |
сумма квадратов |
|
N |
|
2 |
равна |
нулю, |
перед нами |
функциональная |
отклоненийå( yi - yi ) |
|
|||||
i=1
связь между переменными (все точки легли на аппроксимирующую линию), поскольку общая сумма квадратов отклонений в данном случае равна сумме квадратов отклонений объясненных регрессий
(или |
факторной сумме |
квадратов |
отклонений). Наоборот, когда |
y = y , |
то объясненная |
регрессией |
сумма квадратов отклонений |
равна нулю (общая сумма квадратов отклонений равна остаточной сумме квадратов отклонений), связь между переменными модели отсутствует, и вариация игрек обусловлена воздействием иных факторов, не включенных в данную модель.
Приведение сумм квадратов на соответствующее им число степеней свободы дает величины средних квадратов отклонений или конкретных дисперсий – общей, факторной и остаточной, а также приводит их к сопоставимому виду.
Dy = N1-1×åN ( yi - yi )2 общая дисперсия,
i=1
|
N |
|
2 |
|
|
|
|
|
å( yi − y) |
|
N |
|
|
- факторная дисперсия; |
|
Dфакт = i=1 |
|
|
= å( yi − y) |
2 |
|||
|
|
|
|
|
|
|
|
|
|
1 |
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
170
|
1 |
N |
|
2 |
|
Dост = |
|
×å( yi - yi ) |
|
- остаточная дисперсия. |
|
N - 2 |
|
||||
|
i=1 |
|
|
|
|
Число степеней свободы зависит от числа наблюдений N и числа определяемых параметров. Продемонстрируем это на примере - определим число степеней свободы для общей суммы квадратов отклонений. Пусть у нас имеется вариационный ряд переменной игрек: 10; 12; 17; 25; 31; 40. Среднее арифметическое составило 22,5. Находим отклонения 10-22,5=-12,5; 12-22,5=-10,5;
17-22,5=-5,5; 25-22,5=2,5; 31-22,5=8,5; 40-22,5=17,5. åN ( yi - yi ) =0.
i=1
Спрашивается, сколько независимых отклонений из N возможных существует в данном случае? Ответ: N-1 или для нашего примера 6- 1=5. То есть пять единиц выборки из шести варьируют свободно, а шестое может быть определено, если пять других значений известно. Для нахождения суммы квадратов отклонений (равна нулю) нам достаточно знать величины пяти отклонений, а шестое будет найдено автоматически. Пусть нам неизвестна величина четвертого отклонения (2,5), тогда, зная пять величин, находим: (- 12,5)+(-10,5)+(-5,5)+8,5+17,5=-2,5, а среднее равно нулю. Следовательно, недостающее число равно 2,5. Можно сказать и иначе, что в данном случае несвободно лишь одно отклонение из шести.
Для решения задачи парного регрессионно-корреляционного анализа в программе Excel достаточно ввести исходные данные.
X |
Y |
23 |
45 |
34 |
56 |
43 |
61 |
52 |
70 |
64 |
87 |
76 |
98 |
83 |
109 |
96 |
117 |
111 |
130 |
171
Далее последовательно выбрать «Данные», «Анализ данных», «Регрессия», заполнить окно для икс и для игрек, указать, куда вывести итоги регрессионной статистики и получить ответ.
Таблица №. Итоги регрессионной статистики.
Множественн |
|
|
|
|
|
|
ый R |
0,996201 |
|
|
|
|
|
|
|
|
|
|
|
|
R-квадрат |
0,992417 |
|
|
|
|
|
|
|
|
|
|
|
|
Нормированн |
|
|
|
|
|
|
ый R-квадрат |
0,991334 |
|
|
|
|
|
|
|
|
|
|
|
|
Стандартная |
|
|
|
|
|
|
ошибка |
|
|
|
|
|
|
оценки |
2,760112 |
|
|
|
|
|
|
|
|
|
|
|
|
Наблюдения |
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Дисперсионный анализ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Значим |
|
|
df |
SS |
MS |
F |
ость F |
|
|
|
|
|
|
|
|
|
|
|
|
|
1,11E- |
|
Регрессия |
1 |
6979,56136 |
6979,56 |
916,17 |
08 |
|
|
|
|
|
|
|
|
Остаток |
7 |
53,3275201 |
7,61822 |
|
|
|
|
|
|
|
|
|
|
Итого |
8 |
7032,88888 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
Верхн |
|
Коэффициен |
Стандартн |
-статист |
P- |
Нижни |
ие |
|
ты |
ая ошибка |
ика |
Значение |
е 95% |
95% |
|
|
|
|
|
|
|
|
|
|
|
|
15,3233 |
|
Y-пересечение |
20,85047 |
2,33741850 |
8,9203 |
5E-05 |
5 |
26,377 |
|
|
|
|
|
|
|
Переменная |
|
|
|
|
|
|
X 1 |
1,005749 |
0,03322784 |
30,2683 |
1E-08 |
0,92717 |
1,0843 |
|
|
|
|
|
|
|
SS – сумма квадратов отклонений.
MS – средний квадрат отклонений (сумма квадратов отклонений приведенная на число степеней свободы - df).
Dфакт =6979,56 (df=1).
Dост =7,618 (df=7) (53,3275201/7=7,6182).
Dобщ =879 (df=8) (7032,88888/8=879,1)
F = |
Dфакт |
= |
|
|
|
r2 |
|
×(N - 2) |
|
D |
1 |
|
- r2 |
||||||
|
|
|
|
||||||
|
ост |
|
|
|
|
|
|
|
|
F = |
6979 ,56 |
=916 ,1 |
|||||||
|
7,618 |
|
|
|
|
|
|
||
F = |
0,992417 |
|
×(9 -2) =916 ,1 |
||||||
1-0,992417 |
|||||||||
|
|
|
|||||||
172
Fтабл=5,59 при α=0,05 (то есть при вероятности ошибки 5%). Табличное значение находим так: 1) по горизонтальной оси берем значение 1, что соответствует числу степеней свободы в факторной дисперсии (определялось числом независимых переменных – 1 независимая переменная); 2) по вертикали берем 7, что соответствует числу степеней свободы в остаточной дисперсии.
Поскольку Fтабл< Fфакт , постольку уравнение в целом статистически значимо.
0
P-значение (p-value) – вероятность наступления гипотезы H ,
доверительная вероятность (один из совокупности диагностических
0
показателей) – результат проверки статистической гипотезы H ,
показывающий вероятность того, что данные соответствуют
нулевой гипотезе. Следовательно, чем меньше р-значение, тем
0
ниже вероятность гипотезы H и выше вероятность альтернативной
1
гипотезы H . Например, в нашем случае р-значение для сдвига=
0,00005 или 0,005%, а для коэффициента наклона 0,00000001 или
0,000001%. |
То |
есть и для |
сдвига и |
для наклона |
вероятность |
|
0 |
|
|
|
|
гипотезы |
H |
меньше |
1%, что |
подчеркивает |
значимость |
полученных коэффициентов. Сходный результат о значимости
параметров уравнения будет получен и при исследовании t-
статистики Стьюдента.
173
Найдем стандартные ошибки, а также t-статистики Стьюдента
для коэффициента регрессии и свободного члена:
|
|
|
|
|
|
N |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
mb = |
å(yi - yi ) / N - 2 |
|
|
Dост |
|
=0,03322784 |
||||||||||||||||
|
|
i=1 |
|
|
|
|
|
|
= |
|
|
|
||||||||||
|
|
|
|
N |
|
|
|
|
|
|
N |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
å(xi |
- x) |
2 |
|
|
|
å(xi - x)2 |
||||||||
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
2 |
|
|
|
|
N |
|
2 |
|
|
|
|
|
||
ma = |
|
å(yi - yi ) |
|
å x |
|
|
|
|
=2,33741850 |
|||||||||||||
|
|
i=1 |
|
× |
|
|
N |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
N - 2 |
|
|
|
N × å(xi - x)2 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
||
ta = |
a |
|
|
= |
|
20,85047 |
|
=8,92 |
|
|
|
|
|
|
|
|||||||
ma |
|
|
2,3374185 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
tb = |
|
b |
|
= |
|
1,005749 |
|
|
= 30,268 |
|
|
|
|
|
||||||||
|
mb |
|
|
0,03322784 |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Табличное значение t-статистики Стьюдента для 7 степеней свободы при α=0,05 (то есть при вероятности ошибки 5%) соответствует значению 2,3646. Как видно фактические значения, как для свободного члена, так и для коэффициента регрессии больше табличного, что подтверждает их статистическую значимость и статистическую гипотезу H1.
Доверительный интервал для коэффициента регрессии находим по формуле: b ±t ×mb . Соответственно верхняя граница доверительного интервала:
b +t ×mb =1,005749 + 2,3646 ×0,03322784 =1,0843 . Обратите внимание, что
вместо t подставляется значение t табличного (2,3646), а не фактическое значение t-статистики для данного параметра (параметра регрессии) (30,268).
Нижняя граница доверительного интервала:
b -t ×mb =1,005749 -2,3646 ×0,03322784 = 0,92717 .
Отсюда делаем вывод, что в генеральной совокупности реальный коэффициент наклона (β), вероятно, находится в пределах от 0,927 до 1,084.
174
Доверительный интервал для коэффициента сдвига находим по формуле: a ±t ×ma . Соответственно верхняя граница доверительного интервала:
a +t ×ma = 20,85047 + 2,3646 ×2,33741850 = 26,377 .
Нижняя граница доверительного интервала:
a -t ×ma = 20,85047 -2,3646 ×2,33741850 =15,323 .
Отсюда делаем вывод, что в генеральной совокупности реальный коэффициент сдвига (α), вероятно, находится в пределах от 15,323 до 26,377.
Нормированный R-квадрат (иногда обозначают R 2 или adjusted R2 ) приведенный в таблице «итоги регрессионной статистики» (в литературе называют также: скорректированный R2 ; исправленный R2 ) используется для уточнения коэффициента детерминации при включении в модель дополнительных независимых переменных (x). Это обусловлено тем, что добавление новой независимой переменной к уравнению регрессии повышает значение коэффициента детерминации R2. Нормированный R- квадрат сдерживает рост коэффициента детерминации при включении в модель дополнительных независимых переменных.
Нормированный |
|
|
R-квадрат |
|
|
|
– |
|
это |
показатель, |
который |
||||||||||||
рассчитывается по формуле: |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
2 =1- (1- R2 ) × |
|
|
n −1 |
|
= R2 |
- |
|
k |
|
|
×(1- R2 ) или |
что |
также |
||||||||
R |
|
|
|
||||||||||||||||||||
|
n - k -1 |
n - k |
-1 |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
равнозначно: |
|
|
2 = |
|
n −1 |
|
|
× R2 |
- |
|
k |
|
|
|
, |
где k |
– число независимых |
||||||
R |
|
|
|
|
|
||||||||||||||||||
|
|
n - k -1 |
|
n - k -1 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
переменных. |
С |
|
ростом |
числа независимых переменных (k) |
|||||||||||||||||||
увеличивается значение дроби: |
|
|
|
k |
|
, и возрастает корректировка |
|||||||||||||||||
n -k -1 |
|||||||||||||||||||||||
коэффициента детерминации R2 в сторону уменьшения.
Вообще говоря, R2, как и нормированный R-квадрат рассматриваются, как диагностические критерии качества регрессионной модели наряду с другими диагностическими показателями (проверка статистической значимости уравнения в целом, статистической значимости параметров и коэффициента
корреляции). |
|
2 |
|
|
|
2 |
Du |
|
æ |
|
Kx, y |
ö2 |
, где |
- дисперсия в ряду |
||||
|
|
Dy |
|
|
ç |
|
÷ |
|||||||||||
|
R |
|
= |
|
|
=1- |
|
|
= |
|
|
|
|
|
|
|
|
Dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
D |
y |
D |
y |
ç |
|
D |
× D |
÷ |
|
|||||||
|
|
|
|
|
|
|
|
è |
|
y |
|
x ø |
|
|
||||
175
