Предварительные условия для сравнения моделей AIC

26

Какие именно предварительные условия необходимо выполнить для сравнения моделей AIC для работы?

Я только пришел к этому вопросу, когда я сделал сравнение, как это:

> uu0 = lm(log(usili) ~ rok)
> uu1 = lm(usili ~ rok)
> AIC(uu0)
[1] 3192.14
> AIC(uu1)
[1] 14277.29

Таким образом, я оправдал logпреобразование переменной usili. Но я не знаю, могу ли я сравнить модели AIC, когда, например, зависимая переменная отличается?

Идеальный ответ будет включать в себя список предпосылок (математические предположения).

любознательный
источник

Ответы:

29

Вы не можете сравнить две модели, так как они не моделируют одну и ту же переменную (как вы правильно узнаете себя). Тем не менее AIC должен работать при сравнении как вложенных, так и не вложенных моделей.

Просто напоминание, прежде чем мы продолжим: логарифмическая вероятность Гаусса определяется

журнал(L(θ))знак равно-|D|2журнал(2π)-12журнал(|К|)-12(Икс-μ)TК-1(Икс-μ),

- ковариационная структура вашей модели, | D | количество точек в ваших наборах данных, μ средний ответ и x ваша зависимая переменная.К|D|μИкс

В частности, AIC рассчитывается так, чтобы он равнялся , где k - число фиксированных эффектов в вашей модели, а L - функция вероятности [1]. Он практически сравнивает компромисс между дисперсией ( 2 k ) и смещением ( 2 log ( L ) ) в ваших предположениях моделирования. Таким образом, в вашем случае он будет сравнивать две разные логарифмические структуры правдоподобия, когда дело доходит до смещения. Это потому, что когда вы вычисляете свою логарифмическую вероятность, вы фактически смотрите на два термина: подходящий термин, обозначаемый - 12k2log(L)kL2К2журнал(L), и штраф за сложность, обозначаемый как-1-12(Икс-μ)TК-1(Икс-μ). Таким образом, вы видите, что ваш подходящий термин для двух моделей совершенно различен; в первом случае вы сравниваете остатки из необработанных данных, а в другом случае - остатки зарегистрированных данных.-12журнал(|К|)

Помимо Wikipedia, AIC также определен, чтобы приравнять: [3]; эта форма делает еще более очевидным, почему разные модели с разными зависимыми переменными не сравнимы. RSS-это два случая просто несопоставимы между ними.|D|журнал(рSS|D|)+2К

Оригинальная статья Акаике [4] на самом деле довольно трудна для понимания (я думаю). Он основан на дивергенции KL (грубо говоря, различие между двумя распределениями) и доказывает, как вы можете аппроксимировать неизвестное истинное распределение ваших данных и сравнивать его с распределением данных, которое предполагает ваша модель. Вот почему «чем меньше оценка AIC, тем лучше» ; вы ближе к приблизительному истинному распределению ваших данных.

Итак, чтобы свести воедино все, что нужно помнить при использовании AIC, это три [2,5]:

  1. Вы не можете использовать его для сравнения моделей разных наборов данных.

  2. Вы должны использовать одинаковые переменные ответа для всех моделей-кандидатов.

  3. Вы должны иметь , так как в противном случае вы не получите хорошую асимптотическую последовательность.|D|>>К

Извините, что сообщаю вам плохие новости, но использование AIC, чтобы показать, что вы выбираете одну зависимую переменную вместо другой, не является статистически обоснованной вещью. Проверьте распределение ваших остатков в обеих моделях. Если зарегистрированный случай данных имеет нормально распределенные остатки, а случай необработанных данных - нет, у вас есть все обоснование, которое вам может понадобиться. Вы также можете проверить, соответствуют ли ваши необработанные данные логнормальным нормам, что также может быть достаточным оправданием.

Для строгих математических предположений игра KL расхождение и теория информации ...

Ах, и некоторые ссылки:

  1. http://en.wikipedia.org/wiki/Akaike_information_criterion
  2. Информационный критерий Акаике, Шухуа Ху (презентация, стр. 17-18)
  3. Прикладной многомерный статистический анализ, Johnson & Wichern, 6th Ed. (с. 386-387)
  4. Новый взгляд на идентификацию статистической модели, Х. Акаике, IEEE Транзакции по автоматическому управлению 19 (6): 716–723 (1974)
  5. Урок выбора модели № 1: Информационный критерий Акаике, Д. Шмидт и Э. Макалич, (презентация с.39)
usεr11852 говорит восстановить Monic
источник
Благодарность! Я не понял математику, но я понял суть сообщения. Однако не могли бы вы перечислить все предпосылки, необходимые для сравнения моделей AIC? Просто чтобы быть уверенным, я не сделаю еще одну ошибку в следующий раз. Я пойду и проверю их по одному.
Любопытно
1
|D|пL(θ)θп(Икс|θ) а ваши оценки ML были непротиворечивыми, но я думаю, что эти предположения являются излишним показателем в статье. ..
1
спасибо за добавление списка этих 3 предположений в ответ! Это то, что мне было нужно.
Любопытно
1
Еще раз посмотрим на ваш ответ: ваш пункт 1. «Вы не можете использовать его для сравнения моделей разных наборов данных» . Что вы подразумеваете под «набором данных»? Что если я изменил набор зависимых переменных? Я думаю, что в этом случае AIC должен быть еще сопоставимым? Не могли бы вы обновить свой ответ, чтобы уточнить это?
Любопытно
1
рSSμИкс
11

uu0ΠяYя-12Σяжурнал(Yя)uu0AIC (uu0)+2*sum (log (usili))AIC (uu1)

probabilityislogic
источник
Я не понимаю, что вы следуете за вашей попыткой как-то «исправить» AIC и что вы на самом деле получили (как интерпретировать свой результат). В любом случае, не углубляйтесь в это, это не имеет значения, потому что мой вопрос был о чем-то совершенно ином: каковы общие предпосылки для того, чтобы AIC (фактическая, неоткорректированная) была разумно сопоставимой. Не зацикливайтесь на этом конкретном примере, это всего лишь пример общего.
Любопытно
1
@curious - моя точка зрения такова, что мой «исправленный AIC» - это фактический AIC, а то, что вы получаете от функции AIC, неверно, когда вы сравниваете преобразования «зависимой переменной». Дело в том-2журнал(п(Y|θ)) изменения при трансформации, Иксзнак равног(Y) (например,Иксзнак равноLог(Y)). Вы должны учитывать якобиан этого изменения при использовании AIC. Используемая AIC()вами функция не учитывает это.
вероятностная
@probabilityislogic: Есть ли у вас какие-либо академические ссылки на ваше предложение (AIC (uu0) + 2 * sum (log (usili)))), чтобы я мог цитировать их в научных трудах? Спасибо.
KuJ
3

Взято из Акаике 1978

Этот отрывок из Akaike 1978 содержит цитату в поддержку решения @probabilityislogic.

Akaike, H. 1978. О вероятности модели временных рядов. Журнал Королевского статистического общества. Серия D (Статистик) 27: 217-235.

BJD
источник
1
извините, я не понимаю, что такое «преобразование переменной» и как это связано с моим вопросом. Пожалуйста, объясните, спасибо
Любопытный