Як інтерпретувати ці сюжети ACF та Pacf


20

Далі наведені сюжети ACF та Pacf щомісячних даних. Другий сюжет - це acf з ci.type = 'ma':

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

Наполегливість високих значень в сюжеті ACF, ймовірно, представляє довгострокову позитивну тенденцію. Питання полягає в тому, чи це означає сезонні зміни?

Я намагався бачити різні сайти на цю тему, але не впевнений, чи показують ці сюжети сезонність.

Аналіз ділянок ACF та PACF

Допоможіть інтерпретувати графіки ACF- та PACF

Допоможіть зрозуміти наступну картину АКФ

Автокореляція та часткова інтерпретація автокореляції

Редагувати: нижче наведено графік відставання до 60:

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

Далі наводяться сюжети різниці (my_series):

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

І до відставання 60:

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

Редагувати: Ці дані походять від: Чи це відповідний метод для перевірки сезонних наслідків даних про кількість самогубств? Тут дописувачі не вважали сюжет ACF і Pacf оригінальної чи різницької серії, який варто згадати (тому це не повинно бути важливо). У кількох місцях згадувались лише залишки ACF / Pacf.


1
Чи можете ви додати щось про свої дані (наприклад, основний сюжет)? Ви пробували щось подібне stl()?
gung - Відновіть Моніку

Я намагаюся зрозуміти, як визначити сезонність за сюжетами ACF та Pacf. Чи потрібен для цього перегляд базового сюжету або stl? Чи не можемо ми визначити щось із цих сюжетів?
rnso

1
Це було б добре. Для наочності ваше запитання насправді не в тому, що відбувається з вашими даними, а в тому, що можна зрозуміти з цих сюжетів поодиноко, чи не так?
gung - Відновіть Моніку

1
Так. Мені часто потрібно визначити, чи присутня сезонність у моїх даних, тому я хочу зрозуміти, яку інформацію я можу отримати з сюжетів ACF та Pacf. Сюжети функції stl досить легко зрозуміти, але не ці графіки.
rnso

Ваші дані дійсно мають деяку сезонність. Перегляньте мою відповідь на @javlacalle.
IrishStat

Відповіді:


9

дивлячись на графіки, щоб спробувати заграти дані в здогадану модель аріма, добре працює, коли 1: У даних немає атрибутів / імпульсів / зрушень рівня, місцевих тенденцій часу та жодних сезонних детермінованих імпульсів І 2), коли модель arima має постійні параметри в часі І 3), коли дисперсія помилок від моделі аріма має постійну дисперсію в часі. Коли ці три речі дотримуються .... у більшості наборів даних підручника, що представляють простоту моделювання арими. Коли 1 чи більше з трьох не втримуються .... у кожному реальному наборі даних, який я коли-небудь бачив. Проста відповідь на ваше запитання вимагає доступу до оригінальних фактів (історичних даних), а не до другої описової інформації у ваших сюжетах. Але це лише моя думка!

ВИДАЛЕНО ПІСЛЯ ПРИЙМАННЯ ДАНИХ:

Я був у грецькій відпустці (фактично робив щось інше, ніж аналіз часових рядів) і не зміг проаналізувати ДАНІ про суїцид, але спільно з цим повідомленням. Тепер я підходящий і правильний, що я подаю аналіз, щоб продовжити / довести на прикладі мої коментарі щодо стратегій багатостадійної ідентифікації моделі та недоліків простого візуального аналізу простих кореляційних сюжетів, як "доказ в пудингу".

Ось ACF вихідних даних введіть тут опис зображенняPACF оригінальної серії введіть тут опис зображення. AUTOBOX http://www.autobox.com/cms/ програмне забезпечення, яке я допомогла розробити, використовує евристику для виявлення стартової моделі введіть тут опис зображення. Діагностична перевірка залишків з цієї моделі запропонувала деяке доповнення моделі за допомогою зсуву рівня, імпульсів та сезонного імпульсу. Зауважте, що Зміна рівня виявляється приблизно на періоді 164, що майже ідентично попередньому висновку про період 176 від @forecaster. Всі дороги не ведуть до Риму, але деякі можуть вас зблизити!введіть тут опис зображення. Тестування на сталість параметрів, що відхиляються, змінюється з часом. Перевіривши детерміновані зміни дисперсії помилок, дійшов висновку, що детермінованих змін у дисперсії помилки не виявлено. введіть тут опис зображення. Випробування Box-Cox щодо необхідності перетворення потужності було позитивним, зробивши висновок про необхідність логарифмічного перетворення. введіть тут опис зображення. Остаточна модель тут введіть тут опис зображення. Залишки остаточної моделі не містять автокореляції введіть тут опис зображення. Сюжет остаточних моделей, здається, не містить Гауссових порушень введіть тут опис зображення. Сюжет Актуального / Підходящого / Прогнозування тут введіть тут опис зображенняз прогнозами тутвведіть тут опис зображення


Дякую за вашу відповідь. Чи такі припущення є настільки важливими і завжди настільки сформованими у даних реального світу, що сюжети ACF та Pacf майже ніколи не можуть трактуватися поодиноко?
rnso

Я ненавиджу ніколи не говорити, АЛЕ припущення, які я виклав, суттєво ускладнить процес візуальної ідентифікації, якщо буде порушено. Ваш набір даних чітко (на мої старі очі) є прикладом цього. Визначення початкової моделі, оцінка та повторна ідентифікація на основі залишкової діагностики є багатоетапним процесом, який не є єдиним, і робиться ВСІМО в тривіальних випадках.
IrishStat

Ще раз повторюю, слідуючи за моїм другом stats.stackexchange.com/users/48766/javlacalle : Перевірка наявності імпульсів та зрушень рівня І І імпульсів сезонного періоду, а також місцевих тенденцій часу та постійності дисперсії помилок також необхідна.
IrishStat

(+1) Хороший аналіз даних. Однак як щодо оригінального питання? чи можна визначити сезонність у даних? Можливо, це можна зробити з висновку, який ви показуєте, але я не міг цього зрозуміти.
javlacalle

сезонність присутня в терміні AR (12) в моделі ARIMA та в сезонному пульсі, починаючи з періоду 98 (2003/2)
IrishStat

10

Інтерпретація ACF та PACF

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

ACF та PACF порядку 12 виходять за рамки достовірності. Однак це не обов'язково означає наявність ідентифікованого сезонного малюнка. ACF та PACF інших сезонних замовлень (24, 36, 48, 60) знаходяться в межах довіри. З графічної картини неможливо зробити висновок, чи обумовлена ​​значимість ACF та PACF порядку 12 через сезонність або тимчасові коливання.

Згадане раніше збереження ACF свідчить про те, що спочатку можуть знадобитися розбіжності, щоб зробити дані нерухомими. Однак ACF / PACF з розрізнених серій виглядають підозрілими, негативна кореляція, можливо, була викликана диференціальним фільтром і може бути фактично не підходящою. Дивіться цю публікацію для отримання детальної інформації.

Визначте, чи є сезонність

Аналіз ACF та PACF слід доповнити іншими інструментами, наприклад:

  • Спектр (погляд на АЧХ у частотній області), може виявити періодичність циклів, що пояснюють більшість змінних даних.
  • Підійдіть до базової моделі структурних часових рядів і перевірте, чи дисперсія сезонного компонента близька до нуля відносно інших параметрів (у функції R stats::StructTSта пакеті stsm ).
  • Тести на сезонність на основі сезонних муляжів, сезонних циклів або тих, що описані та впроваджені в X-12 .
  • Перевірка наявності імпульсів і зсувів рівня, про яку згадує IrishStat, також необхідна, оскільки вони можуть спотворювати висновки з попередніх методів (у R пакет для цього можуть бути корисними для цього.

Я додав сюжет до відставання 60. Якою була б команда R для отримання "відмінної серії"? Я додам сюжети для diff (my_series).
rnso

@mso Я додав значні зміни до своєї попередньої відповіді. Команда для відмінного ряду - це функція, diffяку ви використовували.
javlacalle

javlacalle - було два дуже схожі абзаци, і @rnso намагався допомогти, видаливши один. Я видалив те, що, на мою думку, те, що ви хотіли замінити. Не могли б ви перевірити, чи було вилучено правильний абзац?
Glen_b -Встановіть Моніку

@Glen_b дякую за редагування, я вніс деякі зміни.
javlacalle

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