Saya memiliki database peristiwa (yaitu variabel tanggal) dan kovariat terkait.
Peristiwa ini dihasilkan oleh proses Poisson non-stasioner dengan parameter menjadi fungsi yang tidak diketahui (tetapi mungkin linier) dari beberapa kovariat.
Saya pikir paket NHPoisson ada hanya untuk tujuan ini; tetapi setelah 15 jam penelitian yang gagal, saya masih belum tahu cara menggunakannya.
Heck, saya bahkan mencoba membaca kedua buku yang dirujuk: Coles, S. (2001). Pengantar pemodelan statistik nilai ekstrem. Peloncat. Casella, G. dan Berger, RL, (2002). Kesimpulan statistik. Brooks / Cole.
Satu-satunya contoh dalam dokumentasi fitPP.fun tampaknya tidak sesuai dengan pengaturan saya; Saya tidak memiliki nilai ekstrim! Saya hanya punya acara kosong.
Adakah yang bisa membantu saya dengan contoh sederhana pas proses Poisson dengan parameter dengan kovariat tunggal , dan asumsi, bahwa ? Saya tertarik pada estimasi dan . Saya memberikan kumpulan data dua kolom dengan waktu kejadian (misalkan, diukur dalam detik setelah beberapa waktu sewenang-wenang ) dan kolom lain dengan nilai kovariat ?
sumber
Jawaban:
Pas proses Poisson stasioner
Pertama-tama, penting untuk disadari, seperti apa input data yang dibutuhkan NHPoisson.
Terutama, NHPoisson membutuhkan daftar indeks momen acara. Jika kita mencatat interval waktu dan jumlah peristiwa dalam interval waktu, daripada I kita harus menerjemahkannya ke dalam satu kolom tanggal, mungkin "mengolesi" tanggal selama interval yang direkam.
Untuk kesederhanaan saya akan menganggap, bahwa kita menggunakan waktu yang diukur dalam detik, dan bahwa "kedua" adalah satuan alami .λ
Mari kita mensimulasikan data untuk proses Poisson stasioner sederhana, yang memiliki peristiwa per menit:λ=1
The
simNHP.fun
membuat simulasi. Kami menggunakan untuk mendapatkanaux$posNH
, variabel dengan indeks saat-saat acara menembakkan simulasi. Kita dapat melihat bahwa kita memiliki sekitar 60 * 24 = 1440 acara, dengan memeriksa `length (aux $ posNH).Sekarang mari kita rekayasa balik dengan :λ
fitPP.fun
Karena fungsi hanya mengambil indeks peristiwa, maka perlu juga mengukur berapa banyak indeks yang mungkin dimungkinkan. Dan ini adalah bagian yang sangat membingungkan , karena dalam proses Poisson yang sebenarnya dimungkinkan untuk memiliki jumlah kejadian yang tak terbatas (jika hanya ). Tetapi dari perspektif kita perlu memilih unit waktu yang cukup kecil. Kami memilihnya sangat kecil, sehingga kami dapat mengasumsikan maksimum satu peristiwa per unit waktu.λ>0
fitPP
Jadi apa yang kita lakukan sebenarnya adalah bahwa kita mendekati proses Poisson dengan urutan granular peristiwa binomial, setiap peristiwa berlangsung tepat satu unit waktu, analogi dengan mekanisme di mana distribusi Poisson dapat dilihat sebagai batas distribusi binomial dalam hukum peristiwa langka .
Setelah kami memahaminya, sisanya jauh lebih sederhana (setidaknya untuk saya).
Untuk mendapatkan perkiraan kami saya perlu mengambil eksponen parameter dipasang , . Dan di sini datang sepotong informasi penting, yang hilang di manual : dengan kami cocok dengan logaritma dari λ , tidak λ itu sendiri.λ λ λ
beta
exp(coef(out)[1])
NHPoisson
Pas proses Poisson non-stasioner
NHPoisson
cocok dengan model berikut:Sekarang mari kita siapkan proses Poisson non-stasioner.
Sama seperti sebelumnya,
aux$posNH
akan memberi kita indeks peristiwa, tetapi kali ini kita akan melihat, bahwa intensitas peristiwa berkurang secara eksponensial dengan waktu. Dan laju penurunan ini adalah parameter yang ingin kami perkirakan.Penting untuk dicatat, bahwa kita perlu menempatkan
all.seconds
sebagai kovariat, bukanlambdas
. Eksponensial / logaritisasi dilakukan secara internal olehfitPP.fun
. BTW, terlepas dari nilai yang diprediksi, fungsi membuat dua grafik secara default.Bagian terakhir adalah fungsi pisau swiss untuk validasi model
globalval.fun
,.Antara lain, fungsi membagi waktu menjadi interval, masing-masing
lint
sampel panjang, sehingga dimungkinkan untuk membuat grafik kasar yang membandingkan intensitas yang diprediksi dengan intensitas yang diamati.sumber