Mengukur "penyimpangan" untuk Poisson nol-meningkat atau nol-meningkat negatif?

11

Penyimpangan skala, didefinisikan sebagai D = 2 * (log-kemungkinan model jenuh dikurangi log-kemungkinan model pas), sering digunakan sebagai ukuran good-of-fit dalam model GLM. Persen penyimpangan dijelaskan, didefinisikan sebagai [D (model nol) - D (model pas)] / D (model nol), juga kadang-kadang digunakan sebagai analog GLM ke regresi linier, R-kuadrat. Selain dari fakta bahwa distribusi ZIP dan ZINB bukan bagian dari keluarga distribusi eksponensial, saya mengalami kesulitan memahami mengapa peningkatan penyimpangan dan persen penyimpangan dijelaskan tidak digunakan dalam pemodelan zero-inflated. Adakah yang bisa menjelaskan hal ini atau memberikan referensi yang bermanfaat? Terima kasih sebelumnya!

aleanjeo
sumber
pertanyaan yang sangat bagus - Saya ingin tahu ini juga
user2673238

Jawaban:

3

Penyimpangan adalah konsep GLM, model ZIP dan ZINB bukan glm tetapi diformulasikan sebagai campuran distribusi yang terbatas yang merupakan GLM dan karenanya dapat diselesaikan dengan mudah melalui algoritma EM.

Catatan-catatan ini menggambarkan teori penyimpangan secara ringkas. Jika Anda membaca catatan itu, Anda akan melihat bukti bahwa model jenuh untuk regresi Poisson memiliki kemungkinan log

(λs)=i=1,yi0n[yilog(yi)yilog(yi!)]

yang hasil dari plug-in perkiraan .yi=λ^i

Saya akan melanjutkan sekarang dengan kemungkinan ZIP karena matematika lebih sederhana, hasil yang sama berlaku untuk ZINB. Sayangnya untuk ZIP, tidak ada hubungan sederhana seperti di Poisson. The pengamatan th log-kemungkinan adalahi

i(ϕ,λ)=Zilog(ϕ+(1ϕ)eλ)+(1Zi)[λ+yilog(λ)log(yi!)].

yang tidak diamati sehingga untuk mengatasi ini Anda akan perlu mengambil derivatif parsial wrt baik λ dan φ , mengatur persamaan untuk 0 dan kemudian memecahkan λ dan φ . Kesulitan di sini adalah y i = 0 nilai-nilai, ini bisa masuk ke λ atau menjadi φ dan tidak mungkin tanpa mengamati Z i yang menempatkan y i = 0 observasi ke dalam. Namun, jika kita tahu Z iZiλϕλϕysaya=0λ^ϕ^Zsayaysaya=0Zsayanilai kami tidak akan memerlukan model ZIP karena kami tidak akan memiliki data yang hilang. Data yang diamati sesuai dengan kemungkinan "data lengkap" dalam formalisme EM.

Salah satu pendekatan yang mungkin masuk akal adalah bekerja dengan ekspektasi wrt dari kemungkinan log data lengkap, E ( i ( ϕ , λ ) ) yang menghilangkan Z i dan menggantikannya dengan ekspektasi, ini adalah bagian dari apa yang Algoritma EM menghitung (langkah E) dengan pembaruan terkini. Saya tidak tahu ada literatur yang telah mempelajari pendekatan ini untuk e x p e c t e d penyimpangan sekalipun.ZsayaE(saya(ϕ,λ))Zsayaexhalected

Juga, pertanyaan ini ditanyakan terlebih dahulu sehingga saya menjawab posting ini. Namun, ada pertanyaan lain tentang topik yang sama dengan komentar yang bagus dari Gordon Smyth di sini: penyimpangan untuk model poisson majemuk nol, data kontinu (R) di mana ia menyebutkan respons yang sama (ini adalah penjabaran dari komentar yang saya akan katakan) ditambah mereka disebutkan dalam komentar ke posting lain makalah yang mungkin ingin Anda baca. (penafian, saya belum membaca makalah yang direferensikan)

Lucas Roberts
sumber