Критерий Бройша-Пагана
Материал из MachineLearning.
м |
м |
||
Строка 27: | Строка 27: | ||
==Пример== | ==Пример== | ||
- | [[Изображение:het_plot.png|thumb]] | + | [[Изображение:het_plot.png|thumb|Сгенерированный отклик <tex>y_1</tex>]] |
- | [[Изображение:hom_plot.png|thumb]] | + | [[Изображение:hom_plot.png|thumb|Сгенерированный отклик <tex>y_2</tex>]] |
Рассмотрим пример с использованием системы R: | Рассмотрим пример с использованием системы R: |
Версия 14:42, 28 декабря 2013
Критерий Бройша-Пагана (также Бреуша-Пагана, англ. Breusch-Pagan test) — один из статистических критериев для проверки наличия гетероскедастичности (то есть непостоянной дисперсии) случайных ошибок модели линейной регрессии. Применяется, если есть основания полагать, что дисперсия ошибок может зависеть от некоторой совокупности наблюдаемых переменных:
где
Содержание |
Определение
Формулировки проверяемой и альтернативной гипотез выглядят следующим образом:
- остатки гомоскедастичны;
- неверна.
Статистика критерия может быть получена на основе метода множителей Лагранжа и имеет следующий вид:
- .
Вычисление статистики сводится к следующей процедуре[1].
- Шаг 1: Исходная модель оценивается обычным методом наименьших квадратов, вычисляются остатки .
- Шаг 2: Дисперсия ошибки модели (в предположении её гомоскедастичности) оценивается как .
- Шаг 3: Вычисляются стандартизированные остатки .
- Шаг 4: Строится дополнительная регрессия квадратов стандартизированных ошибок на исходные наблюдаемые переменные:
- .
- Шаг 5: , где — коэффициент детерминации построенной на предыдущем шаге регрессии.
При справедливости нулевой гипотезы о гомоскедастичности остатков статистика критерия имеет распределение хи-квадрат с степенями свободы.
Пример
Рассмотрим пример с использованием системы R:
> ## моделируем наблюдаемые переменные > x <- rep(c(-1,1), 50) > ## генерируем гетероскедастичные ошибки > err1 <- rnorm(100, sd=rep(c(1,2), 50)) > ## генерируем гомоскедастичные ошибки > err2 <- rnorm(100) > ## генерируем отклик > y1 <- 1 + x + err1 > y2 <- 1 + x + err2 > ## проводим тест Бройша-Пагана > bptest(y1 ~ x)$p.value BP 0.0007141008 > bptest(y2 ~ x)$p.value BP 0.9464273
Реализации
- MatLab: встроенной реализации нет, есть реализации на File Exchange.
- R: функция
bptest
в стандартном пакетеlmtest
иncvtest
в пакетеcar
.
Ссылки
- Breusch T.S., Pagan A.R. (1979). "Simple test for heteroscedasticity and random coefficient variation". Econometrica, 47(5), 1287-1294.
- EnWiki: Breusch–Pagan test
- Магнус Я.Р., Катышев П.К., Пересецкий А.А. (2007) "Эконометрика. Начальный курс". М.:Дело, стр. 179-183.
- YouTube: The Breusch Pagan test for heteroscedasticity