Тест множителей Лагранжа

Тест множителей Лагранжа (англ. Lagrange multiplier test, Score test) — статистический тест, используемый для проверки ограничений на параметры статистических моделей, оцененных на основе выборочных данных. Является одним из трёх базовых тестов проверки ограничений наряду с тестом отношения правдоподобия и тестом Вальда. Тест является асимптотическим, то есть для достоверности выводов требуется достаточно большой объем выборки.

Сущность и процедура тестаПравить

Пусть имеется эконометрическая модель с вектором параметров b. Необходимо проверить по выборочным данным гипотезу $H_0:g(b)=0$$H_{0}:~g(b)=0$ , где g — совокупность (вектор) некоторых функций параметров. Идея теста основана на применении известного метода множителей Лагранжа для оценки параметров ограниченной (короткой) модели, исходя из модели без ограничений (длинной модели). Пусть логарифмическая функция правдоподобия для длинной модели равна $l(b)$$l(b)$ . Для оценки короткой модели необходимо построить функцию Лагранжа

$F(b,\mathrm{\lambda <!--\; \lambda \; -->})=l(b)-<!--\; -\; -->{\mathrm{\lambda <!--\; \lambda \; -->}}^{T}g(b)$$F(b,\lambda )=l(b)-\lambda ^{T}g(b)$ 

Тогда условия максимума будут иметь вид:

$\frac{\mathrm{\partial <!--\; \partial \; -->}F(b,\mathrm{\lambda <!--\; \lambda \; -->})}{\mathrm{\partial <!--\; \partial \; -->}b}=\frac{\mathrm{\partial <!--\; \partial \; -->}l(b)}{\mathrm{\partial <!--\; \partial \; -->}b}-<!--\; -\; -->G(b)\mathrm{\lambda <!--\; \lambda \; -->}=0,g(b)=0,G(b)=\frac{\mathrm{\partial <!--\; \partial \; -->}g(b)}{\mathrm{\partial <!--\; \partial \; -->}b}$${\frac {\partial F(b,\lambda )}{\partial b}}={\frac {\partial l(b)}{\partial b}}-G(b)\lambda =0~,~~g(b)=0~,~G(b)={\frac {\partial g(b)}{\partial b}}$ 

Тест основан на том, что если ограничения выполняются, то множители Лагранжа должны быть равны нулю. Поскольку вместо истинных значений параметров будут использоваться оцененные, то множители Лагранжа должны быть просто максимально близки к нулю, а именно — можно показать, что оценки множителей Лагранжа имеют нормальное распределение с нулевым математическим ожиданием и ковариационной матрицей $(G{V}_{\stackrel{^<!--\; ^\; -->}{b}}{G}^T{)}^{-<!--\; -\; -->1}$$(GV_{{{\hat {b}}}}G^{T})^{{-1}}$ , зависящей от ковариационной матрицы оценок параметров длинной модели. Тогда статистика теста

$LM={\mathrm{\lambda <!--\; \lambda \; -->}}^{T}G{V}_{\stackrel{^<!--\; ^\; -->}{b}}{G}^T\mathrm{\lambda <!--\; \lambda \; -->}$$LM=\lambda ^{T}GV_{{{\hat {b}}}}G^{T}\lambda$ 

будет иметь распределение Хи-квадрат, с q степенями свободы, где q — количество ограничений.

Частные случаиПравить

При проверке линейных ограничений для линейной модели регрессии LM-статистика будет равна

Можно показать, что для классической линейной модели LM-статистика равна

$LM=\frac{ES{S}_S-<!--\; -\; -->ES{S}_L}{ES{S}_S/n}$$LM={\frac {ESS_{S}-ESS_{L}}{ESS_{S}/n}}$ 

В частности при проверке значимости регрессии в целом (то есть при проверке гипотезы о равенстве нулю всех коэффициентов при факторах, кроме константы) $ES{S}_S=TSS$$ESS_{S}=TSS$  — общей сумме квадратов (дисперсия зависимой переменной умноженной на n). Следовательно,

$LM=\frac{TSS-<!--\; -\; -->ESS}{TSS/n}=n(1-<!--\; -\; -->ESS/TSS)=n{R}^2$$LM={\frac {TSS-ESS}{TSS/n}}=n(1-ESS/TSS)=nR^{2}$ ,

где $R^2$$R^2$  — коэффициент детерминации.

Взаимосвязь с другими тестамиПравить

Доказано, что тест Вальда (W), тест отношения правдоподобия (LR) и тест множителей Лагранжа (LM) — асимптотически эквивалентные тесты (LM=LR=W). Тем не менее, для конечных выборок значения статистик не совпадают. Для линейных ограничений доказано неравенство $LM⩽<!--\; ⩽\; -->LR⩽<!--\; ⩽\; -->W$$LM\leqslant LR\leqslant W$ . Тем самым тест множителей Лагранжа будет чаще других тестов принимать нулевую гипотезу об ограничениях (реже других будет отвергать её). В случае нелинейных ограничений первая часть неравенства выполняется, а вторая — вообще говоря, нет.

Вместо LM-теста можно использовать асимптотический F-тест, статистика которого связана с LM-статистикой следующим образом:

$F=\frac{LM}{n-<!--\; -\; -->LM}\frac{n-<!--\; -\; -->k}{q}\stackrel{d}{\to }F(q,n-<!--\; -\; -->k)$$F={\frac {LM}{n-LM}}{\frac {n-k}{q}}~{\xrightarrow {d}}~F(q,n-k)$ ,

где k-количество параметров модели.

Во многих случаях на малых выборках такой тест даже предпочтительней исходного LM-теста.

См. такжеПравить

• Тест Вальда
• Тест отношения правдоподобия
• F-тест

ЛитератураПравить

• Магнус Я. Р., Катышев П. К., Пересецкий А. А. Эконометрика. — М.: Дело, 2004. — 576 с.
• William H. Greene. Econometric analysis. — New York: Pearson Education, Inc., 2003. — 1026 с.