Saya ingin menghasilkan deret waktu sintetis. Rangkaian waktu harus berupa rantai markov dengan distribusi marginal gamma dan parameter AR (1) dari . Dapatkah saya melakukan ini hanya dengan menggunakan distribusi gamma sebagai istilah derau dalam model AR (1), atau apakah saya perlu menggunakan pendekatan yang lebih canggih?
time-series
markov-process
gamma-distribution
autoregressive
ahli hidrologi
sumber
sumber
Jawaban:
Orang mungkin menebak (begitu juga saya awalnya) bahwa ya, tetapi proses AR (1) akan memiliki parameter baru. Untuk bentuk dan skala , biarkan . Tulis .a s gt∼Γ(a,s) g~t=gt−E(gt)
Kemudian, suatu AR (1) dalam , juga dapat ditulis sebagai Ingat kembali dan . Dengan properti dari AR (1) -proses, dan Memecahkan sistem persamaan dari dua momen pertama distribusi gamma untuk dua parameternya menghasilkan parameter bentuk baru , dan .gt yt=ρyt−1+gt
Namun argumen ini tidak lengkap karena tidak menunjukkan bahwa memang . Pada dasarnya, tuliskan representasi sehingga dapat dilihat sebagai seri tertimbang gamma diremehkan. Pembacaan saya terhadap posting seperti ini (lihat juga jawaban yang lebih baru lainnya) menunjukkan bahwa ini bukan varian gamma.yt Γ M.A ( ∞ )
Yang mengatakan, sedikit simulasi menunjukkan bahwa pendekatan tersebut menghasilkan perkiraan yang cukup baik:
sumber
plot(grid,dgamma(grid,shape=shape.y,scale=scale.y), lwd=2, col="red", type = "l")
danlines(density(y), type="l", col="lightblue", lwd=2)
bagaimanapun memang menunjukkan bahwa ada perbedaan bahkan sangat besarn
, ketika penaksir kepadatan kernel daridensity
harus OK.Saya sekarang memiliki jawaban untuk pertanyaan yang saya ajukan, tetapi itu membawa saya ke pertanyaan lebih lanjut.
Jadi, pertama, solusinya adalah sebagai berikut:
Untuk Rantai Markov stasioner dengan aΓ[α,p] distribusi marjinal, fungsi kepadatan probabilitas Pt di x diberikan oleh:
kemudian pdf bersyarat dariPt+1 di x diberikan $ P_t = u adalah:
dimanaIν menunjukkan fungsi Bessel yang dimodifikasi. Ini memberikan Rantai Markov dengan distribusi marginal gamma, dan struktur korelasi AR di manaρ(1) adalah ρ .
Rincian lebih lanjut dari ini diberikan dalam makalah yang sangat baik oleh David Warren, yang diterbitkan pada tahun 1986 di Journal of Hydrology, "Skewness Outflow dalam reservoir linier non-musiman dengan aliran yang terdistribusi gamma" (Volume 85, pp127-137; http: // www.sciencedirect.com/science/article/pii/0022169486900806# ).
Ini bagus, karena menjawab pertanyaan awal saya, namun, sistem yang ingin saya wakili dengan PDF ini memerlukan pembuatan seri sintetik. Jika parameter bentuk dan skala distribusi besar, maka ini mudah. Namun, jika saya ingin parameternya kecil maka saya tidak dapat menghasilkan seri dengan karakteristik yang sesuai. Saya menggunakan MATLAB untuk melakukan ini dan kodenya adalah sebagai berikut:
Jika saya menggunakan angka yang jauh lebih besar untuk parameter distribusi gamma maka marginal keluar tepat, tapi saya perlu menggunakan nilai kecil. Adakah pemikiran tentang bagaimana saya bisa melakukan ini?
sumber
Ada sejumlah cara untuk mendapatkan proses Markov urutan pertama dengan margin gamma. Referensi yang sangat baik tentang topik ini adalah makalah oleh GK Grunwald, RJ Hyndman dan LM Tedesko: Pandangan terpadu model AR (1) .
Seperti yang akan Anda lihat, "bentuk inovasi" klasikyt=ϕyt−1+εt bukan cara termudah untuk menentukan transisi Markov p(yt|yt−1) , kecuali kalau ϕ diambil secara acak. Menggunakan distribusi yang dipilih dengan baik; Beta untukϕ dan Gamma untuk
εt , seseorang dapat memperoleh margin gamma.
Proses AR (1) waktu kontinu yang terkenal dengan margin Gamma adalah proses derau-tembakan dengan langkah-langkah eksponensial, yang banyak digunakan misalnya dalam hidrologi dan berkaitan dengan proses Poisson. Ini dapat digunakan dengan sampling waktu diskrit juga, kemudian muncul sebagai koefisien acak AR (1) dengan distribusi tipe campuran untuk inovasi.
sumber
Ide yang diilhami kopula adalah mengubah proses Gaussian AR (1), katakanlah
Sedangkan proses yang dihasilkanyt akan memiliki properti Markov, tidak akan menjadi AR (1), namun, karena fungsi autokorelasi parsialnya tidak terputus untuk kelambatan yang lebih besar dari 1 seperti yang terlihat dalam simulasi berikut:
Jika sebaliknya membiarkanxt menjadi AR (p) dengan koefisien yang sesuai, maka mungkin yt dapat dibuat kira-kira AR (1), yaitu, pilih urutannya p dan ϕ1,…,ϕp sedemikian rupa sehingga pacf dari yt menjadi cukup kecil untuk semua keterlambatan lebih tinggi dari 1. Tapi sekarang prosesnya yt tidak akan lagi memiliki properti Markov.
sumber