- •Парная регрессия. Точечные оценки параметров в условиях классической нормальной модели.
- •Интервальные оценки параметров линейной регрессии.
- •Интервальная оценка средних и индивидуальных значений объясняемой переменной при заданном значении объясняющей переменной.
- •Проверка гипотезы о значимости коэффициентов линейной регрессии. Оценивание качества приближения с помощью коэффициента детерминации.
- •Коэффициент детерминации.
- •Проверка гипотезы о значимости уравнения регрессии в целом.
Коэффициент детерминации.
Введем показатель,
характеризующий качество
аппроксимации
экспериментальных данных с помощью
функции линейной регрессии. Отметим,
что среди всевозможных наблюдаемых пар
данных (xi,yi)
наблюдаемый в статистическом эксперименте
разброс значений переменной y
обязан двум факторам. Во–первых,
переменная y
зависит от x
и, следовательно, меняется потому, что
в разных экспериментах наблюдались
различные значения переменной x.
Во–вторых, значения y
меняются от эксперимента к эксперименту
потому, что y
зависит от случайной величины ,
которая будет принимать различные
значения даже при постоянном значении
x.
Мерой
разброса
(вариацией) значений переменной y
естественно считать выборочную дисперсию
.
Определим, какая доля дисперсии приходится
на детерминированную зависимость y
от x,
а какая – на стохастическую зависимость
yот.
По определению выборочной дисперсии, получаем
Покажем, что второе слагаемое в приведенной сумме можно отбросить (то есть сумма этих слагаемых равна 0). Вспомним, что параметры b0 и b1 были найдены как решения линейной системы уравнений
С учетом этих уравнений получаем:
Следовательно, дисперсия распадается на два слагаемых: остаточную сумму квадратов
и регрессионную дисперсию, связанную с детерминированной частью зависимости,
.
Коэффициентом
детерминации R2
называется отношение этой последней
суммы к суммарной дисперсии
,
то есть
.
С учетом соотношения
,
получаем
.
Окончательная формула:
.
Таким образом, коэффициент детерминации R2 есть доля в вариации переменной Y, за которую отвечает детерминированная часть регрессии. По определению, величина R2 принимает значения от 0 до 1. Если R2=1, то стохастическая часть не оказала влияния на значения Y в статистическом опыте (то есть можно считать, что =0). Наоборот, R2=0 только в случае b0=b1=0.
Проверка гипотезы о значимости уравнения регрессии в целом.
При отсутствии
линейной зависимости y
от x
(то есть в случае y=)
статистики
и
имеют 2–распределение
с числом степеней свободы N2
и 1 соответственно. Действительно, в
этом случае
,
гдеei
– нормально распределенные случайные
величины с нулевым математическим
ожиданием и средним квадратичным
уклонением
.
При этом эти случайные величины ei
подчиняются двум линейно независимым
соотношением, следовательно, при
достаточно большихN,
статистика
распределена по закону 2
с N2
степенями свободы.
Аналогично, для статистики DR имеем:
,
где
,
.
Подставляя в сумму для DR, получим
.
При этом b1
нормально распределена с математическим
ожиданием a1
= 0 и средним квадратичным
уклонением
.Следовательно,
статистика
распределена по закону 2
с 1 степенью свободы.
Теперь мы можем найти закон распределения коэффициента детерминации и сформулировать критерий незначимости уравнения регрессии с заданным уровнем значимости. Для этого нам потребуется знание еще одного распределения, а именно, распределение Фишера–Снедекора, или F–распределения.
Рассмотрим случайную величину
,
где
– независимые случайные величины с
числом степеней свободы ki.
Случайная величина
называется
F–распределением
с числом степеней свободы
.
Функция плотности F–распределения
имеет вид
.
Можно показать, что статистики и являются независимыми, поэтому случайная величина
имеет F–распределение с числом степеней свободы 1, N2. Это дает следующий критерий проверки гипотезы о незначимости уравнения регрессии.
Зададимся уровнем
значимости
и найдем критическую точку
F–распределения
с числом степеней свободы 1, N2
и данным уровнем значимости . Для этого
следует решить уравнение
.
Для вычисления можно использовать следующие команды программ как EXCEL, MATHCAD и MAPLE.
>with(stats):statevalf[icdf,fratio[k1,k2]](1-); (Maple)
qF(1,k1,k2) (MathCad)
=FРАСПОБР(;k1;k2) (Excel)
В
частности,
.
Уравнение
регрессии признается незначимым на
уровне значимости ,
если
.
Иначе, гипотеза о незначимости уравнения регрессии отвергается.
