Яку модель можна використовувати, коли припущення про постійну дисперсію порушено?


9

Оскільки ми не можемо підходити до моделі ARIMA, коли припущення про постійну дисперсію порушується, яку модель можна використовувати для пристосування одновимірних часових рядів?


Якщо припустити, що в пристосованій моделі немає незалежних регресорів, непостійна дисперсія насправді є лише проблемою, коли дисперсія терміна помилки залежить від часу. Потім: arma + garch
user603

Відповіді:


8

Існує ряд варіантів моделювання для обліку непостійної дисперсії, наприклад, ARCH (і GARCH та їх безліч розширень) або стохастична модель мінливості.

Модель ARCH розширює моделі ARMA з додатковим рівнянням часового ряду для терміна квадратної помилки. Вони, як правило, досить легко оцінити (наприклад, пакет fGRACH R).

Моделі SV поширюють моделі ARMA додатковим рівнянням часового ряду (зазвичай AR (1)) для журналу часової дисперсії. Я знайшов, що ці моделі найкраще оцінюються за допомогою байєсівських методів (OpenBUGS працював для мене добре в минулому).


5

Можна підходити до моделі ARIMA, але спочатку потрібно стабілізувати дисперсію, застосувавши відповідне перетворення. Ви також можете використовувати перетворення Box-Cox. Це було зроблено в книзі Аналіз часових рядів: із програмами R , стор. 99, а потім вони використовують перетворення Box-Cox. Перевірте це посилання Моделювання Box-Jenkins Інше посилання - сторінка 169, Вступ до часових рядів та прогнозування, Броквелл та Девіс, «Після того, як дані були трансформовані (наприклад, за допомогою якоїсь комбінації Box-Cox та диференціальних перетворень або шляхом усунення тренду та сезонних компонентів) до точки, де трансформована серія X_t потенційно може бути пристосована до нульової середньої моделі ARMA, ми стикаємося з проблемою вибору відповідних значень для порядків p і q ». Тому вам потрібно стабілізувати дисперсію до того, як відповідати моделі ARIMA.


1
Я не бачу, як спочатку можна зробити стабілізацію дисперсії. Вам потрібно спочатку побачити залишки з моделі, щоб побачити, чи змінюється залишкова дисперсія з часом. Тоді перегляд залишків може підказати, як змінити модель або стабілізувати дисперсію.
Майкл Р. Черник

Просто побудувавши часовий ряд, ви зможете дізнатися, чи слід застосовувати стабілізацію дисперсії чи ні. Це було зроблено в книзі "Аналіз часових рядів з додатками в R", сторінка 99, а потім вони використовують перетворення Box-Cox. Ви можете перевірити це самостійно. Якщо ви підходите, не стабілізуючи дисперсію, то вона буде показана на графіку залишку. Вся справа в тому, що ми повинні спробувати виправити будь-які порушення в припущенні до моделі ARIMA, перш ніж їх встановити. Я настійно пропоную вам бути обережнішими, коли давати негативні бали на відповідь! Удачі.
Стат

Так, я був тим, хто спростував вашу відповідь. Я погоджуюся, що ви можете отримати відчуття дисперсії неоднорідності з сюжету серії. Але я все ще не вважаю, що доцільно застосувати трансформацію, стабілізуючу трансформацію, перед тим, як випробувати моделі. Усі моделі орієнтовні. Ви підходите, дивіться на залишки та змінюєте за потребою. Це підхід три кроки Box-Jenkins. Початкова ідентифікація моделі з подальшим встановленням, а потім діагностична перевірка повторюваним циклом, якщо модель не здається адекватною.
Майкл Р. Черник

Це означає, що Ви уважно не читали Box-Jenkins. Перевірте це посилання robjhyndman.com/papers/BoxJenkins.pdf Інша посилання, стор. 169, Вступ до часових рядів та прогнозування, Броквелл та Девіс, “Після того, як дані були перетворені (наприклад, за допомогою деякої комбінації Box-Cox та диференціальних перетворень або шляхом усунення трендових та сезонних компонентів) до точки, коли трансформована серія X_t потенційно може бути пристосована до нульової середньої моделі ARMA, ми стикаємося з проблемою вибору відповідних значень для замовлень p і q ». Ви можете просто визнати, що ви помилилися.
Стат

Stat і @Michael, у вас обох є дійсні моменти: Stat, тому що часто чітко вказується початкова трансформація Box-Cox - так чому б не почати процес ітеративного моделювання, попередньо застосовуючи цю трансформацію? - все ж Майкл також має право зазначити, що в центрі уваги слід звертати увагу на залишки моделі, а не на значення, що залежать від сировини (відмінність тут часто не зрозуміла в питаннях). Ні для голосування, ні для звинувачення у помилках для проведення цього обговорення не потрібно. Якщо ви збираєтесь сперечатися, зробіть це з того, про що ви обоє по-справжньому не погоджуєтесь!
шурхіт

2

Я б спершу запитав, чому залишки моделі ARIMA не мають постійної дисперсії, перш ніж я відмовляться від підходу. Чи не мають залишки в нихдерева кореляційної структури? Якщо вони є, можливо, в модель повинні бути включені деякі середні середні терміни.

Але тепер припустимо, що залишки не мають будь-якої структури автокореляції. то якими способами змінюється дисперсія з часом (збільшується, зменшується або коливається вгору і вниз)? Те, як змінюється дисперсія, може бути підказкою, що не так у існуючій моделі. Можливо, є коваріати, які взаємопов'язані з цим часовим рядом. У такому випадку коваріати можуть бути додані до моделі. Залишки можуть не проявляти незмінну дисперсію.

Ви можете сказати, що якщо ряд перехресно співвідноситься з коваріатом, який відображається в автокореляції залишків. Але це не було б, якщо кореляція здебільшого знаходиться на відставанні 0.

Якщо ні додавання ковзаючих середніх термінів, ні введення коваріатів не допомагають вирішити проблему, можна, можливо, розглянути можливість визначення функції, що змінюється за часом, для залишкової дисперсії на основі кількох параметрів. Тоді це співвідношення може бути включено у функцію вірогідності з метою зміни модельних оцінок.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.