Проверка эквивалентности не вложенных моделей

12

Скажем, - линейная функция от и фиктивная . Моя гипотеза состоит в том, что само по себе , как гедонистическому индексу вектора других переменных, . У меня есть поддержка для этого в из (т.е. , , ..., ) на . Есть ли способ проверить эквивалентность этих двух моделей:x d d Z M A N O V A Z z 1 z 2 z n dyxddZMANOVAZz1z2znd

Модель 1:y=b0+b1x+b2d+e1

Модель 2:y=g0+ZG+e2

где - вектор-столбец параметров.G

user3671
источник

Ответы:

8

Для начала вы должны определить понятие эквивалентности . Кто-то может подумать, что две модели эквивалентны, когда они дают почти одинаковую точность прогнозирования (эта будет иметь отношение к временным рядам и панельным данным), а другая может быть заинтересована, если совпадения с моделью близки . Первый является объектом для различной перекрестной проверки (обычно с помощью домкратного ножа или некоторых тестов вне выборки, Роб accuracy()делает это хорошо), второй идет для минимизации некоторого информационного критерия.

В микроэкономике выбор , хотя вы также можете рассмотреть если вы работаете с небольшими размерами выборки. Обратите внимание, что выбор, основанный на минимизации информационного критерия, также актуален для вложенных моделей.BICAIC

Приятное обсуждение дано в должн--то заказать Камерон и Триведите (Глава 8.5 обеспечивает отличный обзор методов), более конкретные теоретические детали находятся в Гонконге и Престоне здесь .

Грубо говоря, выбор из двух моделей более экономный (с меньшим количеством параметров для оценки, следовательно, с большей степенью свободы) будет предложен в качестве предпочтительного. Информационный критерий вводит специальную функцию штрафа, которая ограничивает включение дополнительных объясняющих переменных в линейную модель, концептуально аналогичную ограничениям, введенным скорректированным .R2

Однако вас может не заинтересовать выбор модели, которая минимизирует выбранный информационный критерий. Концепция эквивалентности подразумевает, что некоторые тестовые статистические данные должны быть сформулированы. Таким образом , вы можете пойти для отношения правдоподобия испытаний либо Кокса или Voung тесты, Дэвидсон-МакКиннон теста. JLRJ

Наконец, согласно тегам, вас могут заинтересовать только Rфункции:

library(lmtest)
coxtest(fit1, fit2)
jtest(fit1, fit2)

Где fit1и fit2два невложенных подогнанные линейные модели регрессии, coxtestявляется Cox тест, и Дэвидсон-МакКиннон тест.JLRjtestJ

Дмитрий Челов
источник
Спасибо, Дмитрий. Если я правильно понимаю, и coxtest, и jtest являются существенно измененными вложенными тестами. Шаг 1: Запустите модель с объединенным пулом регрессоров из модели 1 и модели 2. Шаг 2: Проверьте каждую модель 1 и модель 2 отдельно как подмножества «супермодели». Я прав? Кроме того, в примечании мер IC, есть ли способ статистически сравнить различия AIC / BIC между моделями 1 и 2? Примечание: я НЕ пытаюсь выбрать «лучшую» модель, но вы правы в том, что я пытаюсь проверить, имеют ли две модели одинаковые соответствия.
user3671
@user, вам не нужно добавлять в супермодель что-либо, просто поставьте jtestили coxtestне вложенные подгонки из Step1. Информационный критерий для не вложенных элементов будет хорошим руководством для определения того, какая модель является более статистически подходящей (экономной), но для проверки гипотез я бы просто пошел на любой (на самом деле логарифмическая вероятность является частью любого информационного критерия). , Выводы будут несколько близки, но, поскольку существуют две детерминистически заданные штрафные функции, их немного сложно сравнить статистически. LR
Дмитрий Челов