Saya mencoba memodelkan dan memperkirakan deret waktu yang bersifat siklik daripada musiman (yaitu ada pola yang mirip musim, tetapi tidak dengan periode tetap). Ini harus dimungkinkan untuk dilakukan dengan menggunakan model ARIMA, sebagaimana disebutkan dalam Bagian 8.5 dari Peramalan: prinsip dan praktik :
Nilai penting jika data menunjukkan siklus. Untuk mendapatkan prakiraan siklus, perlu memiliki bersama dengan beberapa kondisi tambahan pada parameter. Untuk model AR (2), perilaku siklik terjadi jika .p ≥ 2 ϕ 2 1 + 4 ϕ 2 < 0
Apa kondisi tambahan ini pada parameter dalam kasus ARIMA umum (p, d, q)? Saya belum dapat menemukannya di mana pun.
Jawaban:
Beberapa intuisi grafis
Pada model AR , perilaku siklik berasal dari akar konjugat kompleks hingga polinomial karakteristik. Untuk memberikan intuisi, saya telah merencanakan fungsi respons impuls di bawah ini ke dua contoh model AR (2).
Untuk , Root dari polinomial karakteristik adalah mana adalah nilai eigen dari matriks saya definisikan di bawah ini. Dengan nilai eigen konjugasi kompleks dan , mengontrol redaman (di mana ) dan mengontrol frekuensi gelombang kosinus.j=1…,p 1λj λ1,…,λp A λ=reiωt λ¯=re−iωt r r∈[0,1) ω
Detail AR (2) contoh
Mari kita asumsikan kita memiliki AR (2):
Anda dapat menulis AR (p) apa saja sebagai VAR (1) . Dalam hal ini, representasi VAR (1) adalah:
Perhatikan bahwa . Membentuk dekomposisi nilai eigen dan menaikkan ke kekuatan .E[Xt+k∣Xt,Xt−1,…]=AkXt A k Ak=[λ11λ21][λk100λk2]⎡⎣1λ1−λ2−1λ1−λ2−λ2λ1−λ2λ1λ1−λ2⎤⎦
Nilai eigen nyata mengarah ke peluruhan saat Anda menaikkan . Nilai eigen dengan komponen imajiner non-nol mengarah pada perilaku siklik.λ λk
Nilai eigen dengan kasus komponen imajiner:ϕ21+4ϕ2<0
Dalam konteks AR (2), kami memiliki nilai eigen yang kompleks jika . Karena nyata, mereka harus berpasangan yang merupakan konjugat kompleks satu sama lain.ϕ21+4ϕ2<0 A
Mengikuti Bab 2 Prado dan Barat (2010), marict=λλ−λ¯yt−λλ¯λ−λ¯yt−1
Anda dapat menunjukkan ramalan diberikan oleh:E[yt+k∣yt,yt−1,…]
Berbicara secara longgar, menambahkan konjugat kompleks membatalkan komponen imajiner mereka sehingga membuat Anda dengan gelombang kosinus tunggal teredam dalam ruang bilangan real. (Catatan kita harus memiliki untuk stasioneritas.)0≤r<1
Jika Anda ingin menemukan , , , , mulailah dengan menggunakan rumus Euler yang , kita dapat menulis:r ω at θt reiθ=rcosθ+rsinθ
Lampiran
Catatan Peringatan terminologi yang membingungkan! Menghubungkan polinomial karakteristik A dengan polinom karakteristik AR (p)
Trik time-series lainnya adalah menggunakan operator lag untuk menulis AR (p) sebagai:
Ganti operator lag dengan beberapa variabel dan orang sering menyebut sebagai polinom karakteristik dari model AR (p). Seperti yang dibahas dalam jawaban ini , ini adalah persis polinomial karakteristik mana . Akar adalah kebalikan dari nilai eigen. (Catatan: untuk model yang tidak bergerak yang Anda inginkan , yang ada di dalam unit cirlce, atau ekuivalen , yang berada di luar lingkaran unit.)L z 1−ϕ1z−…−ϕpzp A z = 1A z=1λ z |λ|<1 |z|>1
Referensi
Prado, Raquel dan Mike West, Time Series: Modeling, Computation, and Inference , 2010
sumber