Сравнение AIC модели и ее лог-преобразованной версии


17

Суть моего вопроса заключается в следующем:

Пусть YRn быть многомерной нормальной случайной величиной со средним и ковариационной матрицей . Пусть , то есть . Как сравнить AIC модели, подходящей для наблюдаемых реализаций с моделью, подходящей для наблюдаемых реализаций Z ?μΣZ:=log(Y)Zi=log(Yi),i{1,,n}YZ



Мой начальный и чуть более длинный вопрос:

Пусть YN(μ,Σ) - многомерная нормальная случайная величина. Если я хочу сравнить модель, подходящую для Y с моделью, подходящей для log(Y) , я мог бы взглянуть на их лог-вероятности. Однако, поскольку эти модели не являются вложенными, я не могу сравнивать правдоподобие правдоподобия (и тому подобное, например, AIC и т. Д.), Но мне нужно их преобразовать.

Я знаю, что если X1,,Xn являются случайными переменными с совместным pdf g(x1,,xn) и если Yi=ti(X1,,Xn) для взаимно-однозначных преобразований ti и i{1,,n} , затем pdf-файл Y1,,Yn задается как

f(y1,,yn)=g(t11(y),,tn1(y))det(J)
где J - якобиан, связанный с преобразованием.

Нужно ли просто использовать правило преобразования для сравнения

l(Y)=log(i=1nϕ(yi;μ,Σ))
до
l(log(Y))=log(i=1nϕ(log(yi);μ,Σ))

или я могу что-то еще сделать?


[править] Забыл поместить логарифмы в двух последних выражениях.


Вы также, кажется, потеряли якобиана в последнем выражении.
whuber

2
Я не понимаю преобразование . Как вы можете взять когда отрицательный? loglogYY
полубрюин

Ответы:


6

Вы не можете сравнить AIC или BIC при установке двух различных наборов данных , то есть и . Вы можете сравнивать только две модели, основанные на AIC или BIC, только при подгонке к одному и тому же набору данных. Взгляните на выбор модели и вывод на основе нескольких моделей: практический информационно-теоретический подход (Burnham and Anderson, 2004). Они упомянули мой ответ на странице 81 (раздел 2.11.3 Преобразования переменной ответа):YZ

Исследователи должны быть уверены, что все гипотезы смоделированы с использованием одной и той же переменной ответа (например, если бы весь набор моделей основывался на log (y), проблем не возникло бы; это неверное смешивание переменных ответа).

И, между прочим, для использования критериев AIC или BIC ваши модели не обязательно должны быть вложенными (то же самое, стр. 88, раздел 2.12.4 Модели без вложенных данных), и на самом деле это одно из преимуществ использования BIC.


5

Akaike (1978, стр. 224) описывает, как AIC может быть скорректирована в присутствии преобразованной исходной переменной для обеспечения возможности сравнения моделей. Он заявляет: «Эффект преобразования переменной представлен просто умножением вероятности на соответствующий якобиан в AIC ... для случая оно равно −2 , где суммирование распространяется на ».log{y(n)+1}log{y(n)+1}n=1,2,...,N

Akaike, H. 1978. «О вероятности модели временных рядов», Журнал Королевского статистического общества, Серия D (The Statistician), 27 (3/4), с. 217–235.


1
есть ли подход в R для этого?
лесничий
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.