Cara menginterpretasikan plot ACF dan pacf ini

20

Berikut ini adalah plot acf dan pacf dari seri data bulanan. Plot kedua adalah acf dengan ci.type = 'ma':

masukkan deskripsi gambar di sini

Bertahannya nilai tinggi dalam petak acf mungkin mewakili tren positif jangka panjang. Pertanyaannya adalah apakah ini mewakili variasi musiman?

Saya mencoba melihat berbagai situs tentang topik ini, tetapi saya tidak yakin apakah plot ini menunjukkan musim.

Analisis plot ACF dan PACF

Membantu menafsirkan plot ACF- dan PACF

Bantu memahami gambar ACF berikut

Autokorelasi dan interpretasi autokorelasi parsial

Sunting: berikut ini adalah grafik untuk jeda hingga 60:

masukkan deskripsi gambar di sini

Berikut ini adalah plot diff (my_series):

masukkan deskripsi gambar di sini

Dan hingga 60:

masukkan deskripsi gambar di sini

Sunting: Data ini dari: Apakah ini metode yang tepat untuk menguji efek musiman dalam data jumlah bunuh diri? Di sini kontributor tidak mempertimbangkan plot acf dan pacf dari seri asli atau berbeda yang perlu disebutkan (jadi itu tidak boleh penting). Hanya plot acf / pacf dari residu yang dirujuk di beberapa tempat.

juga
sumber
1
Bisakah Anda menambahkan sesuatu tentang data Anda (misalnya, plot dasar)? Apakah Anda mencoba yang seperti itu stl()?
gung - Reinstate Monica
Saya mencoba memahami cara menentukan musim dari plot acf dan pacf. Apakah peninjauan plot dasar atau stl diperlukan untuk ini? Bisakah kita tidak menentukan sesuatu dari plot ini?
rnso
1
Itu akan baik-baik saja. Untuk lebih jelasnya, pertanyaan Anda sebenarnya bukan tentang apa yang terjadi dengan data Anda, tetapi apakah tentang apa yang dapat dipahami dari plot ini secara terpisah, apakah itu benar?
gung - Reinstate Monica
1
Iya nih. Saya sering perlu menentukan apakah ada data musiman pada musim saya, jadi saya ingin memahami informasi apa yang dapat saya peroleh dari plot acf dan pacf. Plot fungsi stl cukup mudah dimengerti tetapi tidak plot ini.
rnso
Data Anda memang memiliki musiman. Silakan lihat tanggapan saya terhadap @javlacalle.
IrishStat

Jawaban:

9

melihat plot untuk mencoba memasukkan data menjadi model arima yang ditebak bekerja dengan baik ketika 1: Tidak ada perubahan outlier / pulsa / level, tren waktu lokal dan tidak ada pulsa deterministik musiman dalam data DAN 2) ketika model arima memiliki parameter konstan dari waktu ke waktu DAN 3) ketika varians kesalahan dari model arima memiliki varians konstan dari waktu ke waktu. Kapan ketiga hal ini berlaku .... di sebagian besar set data buku teks menyajikan kemudahan pemodelan arima. Kapan 1 atau lebih dari 3 tidak tahan .... di setiap set data dunia nyata yang pernah saya lihat. Jawaban sederhana untuk pertanyaan Anda memerlukan akses ke fakta asli (data historis) bukan informasi deskriptif sekunder di plot Anda. Tapi ini hanya pendapat saya!

DIedit SETELAH PENERIMAAN DATA:

Saya sedang berlibur di Yunani (sebenarnya melakukan sesuatu selain analisis deret waktu) dan tidak dapat menganalisis SUICIDE DATA tetapi dalam hubungannya dengan posting ini. Sekarang sudah tepat dan benar bahwa saya menyerahkan analisis untuk menindaklanjuti / membuktikan dengan contoh komentar saya tentang strategi identifikasi model multi-tahap dan kegagalan analisis visual sederhana dari plot korelasi sederhana sebagai "buktinya ada di puding".

Berikut adalah ACF dari data asli masukkan deskripsi gambar di siniThe PACF dari seri asli masukkan deskripsi gambar di sini. AUTOBOX http://www.autobox.com/cms/ perangkat lunak yang saya bantu kembangkan menggunakan heuristik untuk mengidentifikasi model awal. Dalam hal ini ditemukan model yang diidentifikasi awalnya masukkan deskripsi gambar di sini. Pemeriksaan diagnostik residu dari model ini menyarankan beberapa model augmentasi menggunakan pergeseran level, pulsa dan pulsa musiman. Perhatikan bahwa Level Shift terdeteksi pada atau sekitar periode 164 yang hampir identik dengan kesimpulan sebelumnya tentang periode 176 dari @forecaster. Semua jalan tidak mengarah ke Roma tetapi beberapa bisa membuat Anda dekat!masukkan deskripsi gambar di sini. Pengujian untuk parameter konstan menolak perubahan parameter dari waktu ke waktu. Memeriksa perubahan deterministik dalam varian kesalahan menyimpulkan bahwa tidak ada perubahan deterministik terdeteksi dalam varian kesalahan. masukkan deskripsi gambar di sini. Tes Box-Cox untuk kebutuhan transformasi daya adalah positif dengan kesimpulan bahwa transformasi logaritmik diperlukan. masukkan deskripsi gambar di sini. Model terakhir ada di sini masukkan deskripsi gambar di sini. Sisa dari model akhir tampaknya bebas dari autokorelasi apa pun masukkan deskripsi gambar di sini. Plot residual model akhir tampaknya bebas dari Pelanggaran Gaussian masukkan deskripsi gambar di sini. Plot Aktual / Fit / Prakiraan ada di sini masukkan deskripsi gambar di sinidengan perkiraan di sinimasukkan deskripsi gambar di sini

IrishStat
sumber
Terima kasih atas jawaban anda. Apakah asumsi ini begitu penting dan selalu begitu dirusak dalam data dunia nyata sehingga plot acf dan pacf hampir tidak pernah bisa ditafsirkan secara terpisah?
rnso
Saya benci untuk mengatakan TAPI TETAPI asumsi yang saya paparkan akan sangat menyulitkan proses identifikasi visual jika dilanggar. Kumpulan data Anda dengan jelas (di mata saya yang lama) adalah contohnya. Mengidentifikasi model awal, memperkirakan dan mengidentifikasi kembali berdasarkan sisa diagnostik adalah proses multi-tahap bukan satu dan dilakukan KECUALI dalam kasus-kasus sepele.
IrishStat
Untuk mengulangi mengikuti stats.stackexchange.com/users/48766/javlacalle teman saya : Memeriksa keberadaan pulsa dan pergeseran level DAN pulsa musiman DAN tren waktu lokal DAN konstanta varian kesalahan juga diperlukan.
IrishStat
(+1) Analisis data yang bagus. Namun, bagaimana dengan pertanyaan aslinya? dapat musiman diidentifikasi dalam data? Mungkin bisa disimpulkan dari output yang Anda tunjukkan, tapi saya tidak bisa mengetahuinya.
javlacalle
musiman hadir dalam istilah AR (12) dalam model ARIMA dan dalam pulsa musiman mulai dari periode 98 (2003/2)
IrishStat
10

Interpretasi dari ACF dan PACF

Peluruhan lambat fungsi autokorelasi menunjukkan data mengikuti proses memori panjang. Durasi guncangan relatif persisten dan memengaruhi data beberapa pengamatan di depan. Ini mungkin tercermin oleh pola tren yang halus dalam data.

ACF dan PACF dari pesanan 12 berada di luar batas kepercayaan yang signifikan. Namun, ini tidak selalu berarti adanya pola musiman yang dapat diidentifikasi. ACF dan PACF dari pesanan musiman lainnya (24, 36, 48, 60) termasuk dalam kelompok kepercayaan. Dari grafik, tidak mungkin untuk menyimpulkan apakah signifikansi ACF dan PACF pesanan 12 adalah karena fluktuasi musiman atau sementara.

Kegigihan ACF yang disebutkan sebelumnya menunjukkan bahwa perbedaan pertama mungkin diperlukan untuk membuat data stasioner. Namun, ACF / PACF dari seri yang dibedakan terlihat mencurigakan, korelasi negatif mungkin disebabkan oleh filter pembeda dan mungkin sebenarnya tidak sesuai. Lihat posting ini untuk beberapa detail.

Tentukan apakah ada musim

Analisis ACF dan PACF harus dilengkapi dengan alat-alat lain, misalnya:

  • Spectrum (pandangan ke ACF dalam domain frekuensi), dapat mengungkapkan periodisitas siklus yang menjelaskan sebagian besar variabilitas dalam data.
  • Sesuaikan model deret waktu struktural dasar dan periksa apakah varians komponen musiman mendekati nol relatif terhadap parameter lainnya (dalam fungsi R stats::StructTSdan paket stsm ).
  • Tes untuk musiman, berdasarkan pada boneka musiman, siklus musiman atau yang dijelaskan dan diimplementasikan dalam X-12 .
  • Memeriksa keberadaan pulsa dan pergeseran level sebagaimana disebutkan oleh IrishStat juga diperlukan karena mereka dapat mengubah kesimpulan dari metode sebelumnya (dalam R paket tsoutliers dapat berguna untuk tujuan ini).
javlacalle
sumber
Saya telah menambahkan plot hingga lag 60. Apa yang akan menjadi perintah R untuk mendapatkan "seri berbeda"? Saya akan menambahkan plot untuk diff (my_series).
rnso
@ juga saya telah menambahkan perubahan besar pada jawaban saya sebelumnya. Perintah untuk seri berbeda adalah fungsi diffyang Anda gunakan.
javlacalle
javlacalle - ada dua paragraf yang sangat mirip dan @rnso berusaha membantu dengan menghapus satu paragraf. Saya telah menghapus apa yang saya pikir adalah yang ingin Anda ganti. Bisakah Anda memeriksa apakah paragraf yang benar telah dihapus?
Glen_b -Reinstate Monica
@Glen_b terima kasih untuk pengeditannya, saya telah membuat beberapa perubahan.
javlacalle
@ javlacalle, @IrishStat: silakan lihat edit di pertanyaan saya kembali data asli.
rnso