Парный линейный регрессионный анализ.

Лекция 9. Элементы регрессионного анализа

Парный линейный регрессионный анализ.

Пусть проводится наблюдение над двумерной г.с. (Х, Y).

Если фиксировать значение x случайной величины X, то можно рассмотреть условное математическое ожидание с.в. Y при X=x: M[Y/X=x]. Таким образом, M[Y/X=x] является некоторой детерминированной функцией от x: M[Y/X=x] =j(x). Эта функция называется функцией регрессии Y на X, а график функции y=j(x) кривой регрессии Y на X. Если наблюдаетя с.в. Y при определенных значениях x, то случайную величину Y можно представить в виде Y=j(x)+e, где e – с.в. Пусть наблюдения проводятся при фиксированных значениях x1, x2, …,xn. При этом случайная величина Y приняла соответственно значения y1, y2, …, yn. Тогда можно считать, что имеет место выборка yi = j(xi) + ei, i=1, …, n. В дальнейшем будем считать, что случайные величины ei, i=1, …, n, удовлетворяют следующим условиям.

1) ei (i=1, …, n) распределены по нормальному закону с математическим ожиданием 0 и дисперсией s2;

2) они попарно некоррелированны.

Если функция регрессии j(x) линейна, то говорят, что имеет место линейная регрессионная модель. Рассмотрим подробно эту модель.

Пусть Y= b0+b1x+e, yi= b0+b1xi+ei, (i=1, …, n) и выполняются условия 1) – 2).

Задача корреляционного и регрессионного анализа состоит в следующем.

1) Получить наилучшие точечные и интервальные оценки параметров b0, b1, s линейной модели;

2) Проверить значимость модели;

3) Проверить адекватность модели наблюдаемым данным.

Для нахождения точечных оценок применяется метод наименьших квадратов (сокращенно – МНК).

Обозначим искомое уравнение . По МНК коэффициенты ищут, такие, чтобы принимала минимальное значение сумма

,

где обозначено .

По необходимому условию экстремума частные производные функции S по переменным должны обратиться в нуль в точке минимума.

Итак, решаем систему

Преобразовав систему, получаем

(9.1)

По методу Крамера, получим

 

,

где – выборочный корреляционный момент, – выборочная дисперсия с.в. X.

Преобразовав далее, получим

,

где – выборочный коэффициент корреляции.

В итоге получаем формулы для оценок коэффициентов уравнения линейной регрессии:

, . (9.2)

Коэффициент называется выборочным коэффициентом регрессии. Выборочное уравнение регрессии имеет вид

Введем в рассмотрение следующие суммы.

– сумма квадратов отклонений,

– остаточная сумма квадратов,

– сумма квадратов, обусловленная регрессией,

называется остаточной дисперсией.

Теорема 9.1 .

Если параметр = 0, то линейная модель называется незначимой. Для проверки значимости линейной модели выдвигается основная гипотеза H0: = 0 при альтернативной гипотезе H1: ¹ 0.

Статистика имеет распределение Фишера с 1 и n–2 степенями свободы, если основная гипотеза верна. Таким образом, если выборочное значение Fв больше квантили распределения Фишера, то основная гипотеза отвергается с вероятностью a, то есть на уровне значимости a линейная модель статистически значима.

Доверительные интервалы для коэффициентов , c доверительной вероятностью 1– a имеют вид:

< < ,

< < ,

где – квантиль распределения Стьюдента порядка 1–a./2 со степенью свободы n–2.

Коэффициентом детерминации называется величина . Чем ближе значение коэффициента детерминации к 1, тем лучше линейная модель описывает наблюдаемые данные. Если имеет место линейная регрессионная модель, то выборочный коэффициент корреляции между X, Y .