Mengapa tidak ada model data hitungan satu kali lipat?

8

Saya sedang mengerjakan model data hitung nol-naik menggunakan psclpaket. Saya hanya ingin tahu mengapa tidak ada pengembangan model untuk model data hitungan satu kali lipat! Juga mengapa tidak ada pengembangan model data bimodal, katakanlah zero-and-2-inflated! Setelah saya menghasilkan data Poisson satu-inflasi dan menemukan bahwa baik model glmwith family=poissonatau binomial negatif ( glm.nb) tidak cukup baik untuk menyesuaikan data dengan baik. Jika ada yang bisa menjelaskan pemikiran saya, meskipun eksentrik mungkin, itu akan sangat membantu bagi saya.

r generalized-linear-model zero-inflation poisson-regression

— kewalahan
sumber

3

Saya menulis proposal hibah sekali untuk mengerjakan perangkat lunak untuk ini, tetapi ditolak.

— Peter Flom

3

Saya berharap ada sedikit untuk dikatakan daripada ada sedikit panggilan untuk itu sehingga orang-orang tidak banyak menyebutkannya, tetapi itu tidak seperti itu tidak akan pernah dilakukan; mungkin belum banyak dibahas. Ide-ide tidak benar-benar memperkenalkan banyak yang secara substansial berbeda dari kasus-0 umum.

— Glen_b -Reinstate Monica

12

Model Poisson satu-inflasi untuk hitungan $Y_i$ adalah

\begin{aligned} Pr (Y_{i} = 1) & = π_{i} + (1 - π_{i}) \cdot μ_{i} e^{- μ_{i}} \\ Pr (Y_{i} = y_{i}) & = (1 - π_{i}) \cdot \frac{μ_{i}^{y_{i}} e^{- μ_{i}}}{y_{i}!} when y_{i} \neq 1 \end{aligned}

$\begin{align}\Pr(Y_i = 1) &= \pi_i +(1-\pi_i)\cdot\mu_i\mathrm{e}^{-\mu_i}\\ \Pr(Y_i = y_i) &= (1-\pi_i)\cdot\frac{\mu_i^{y_i}\mathrm{e}^{-\mu_i}}{y_i!} \qquad \text{when } y_i\neq 1 \end{align}$

di mana Poisson berarti & probabilitas Bernoulli terkait dengan prediktor melalui fungsi tautan yang sesuai. Anda dapat menentukan model yang sama untuk mengembang probabilitas untuk nilai apa pun yang Anda pilih. $\mu_i$ $\pi_i$

Namun, nol memiliki tempat khusus (& sekali kontroversial) di antara angka-angka penghitungan — dalam arti mewakili tidak adanya apa pun untuk dihitung. Dan itu adalah perbedaan "tidak ada" vs "sesuatu", daripada perbedaan "satu" vs "hitungan lain" yang cenderung relevan di berbagai fenomena yang ingin kita model: ada satu proses yang memberikan sia-sia, satu , dua, ... hitung & lainnya yang tidak memberi nilai sama sekali.

— Scortchi - Reinstate Monica
sumber

Terima kasih atas tanggapan Anda. Saya setuju dengan bagaimana Anda menjelaskan mengapa jumlah nol menarik begitu banyak bunga sehingga jumlah lainnya dihitung. Namun, saya ingin belajar lebih banyak tentang model data hitungan satu kali . Ketika saya mencoba untuk menulis kode dalam R untuk satu-meningkat model regresi dengan kovariat maka perkiraan yang salah. Saya yakin saya telah membuat kesalahan dalam score functiondan hessian matrix. Bisakah Anda merekomendasikan saya teks / artikel yang dapat membantu saya untuk belajar lebih banyak tentang itu?

— kewalahan

Tidak tahu apa-apa tentang pemasangan model Poisson satu-inflasi secara khusus. Karena satu-inflasi hanya sedikit modifikasi dari nol-inflasi, membaca yang terakhir akan membantu. Saya pikir beberapa perubahan pada zeroinflkode harus melakukannya - mengubah kemungkinan Poisson & skor nol untuk mencocokkan model di atas (atau coba hanya mengubah kemungkinan & tidak melewati skor ke optim). Tentu saja Anda juga bisa bertanya di sini atau di SO jika sesuai untuk referensi atau membantu hal-hal yang membuat Anda macet.

— Scortchi

3

Paket R VGAMmemiliki fungsi vglmyang dapat digunakan untuk memenuhi semua jenis model Poisson-esque. Anda dapat menggunakannya untuk menentukan model satu-inflasi, jadi kira-kira seperti itu vglm(Y~X,family=oipospoisson(),data=data). Lihat di sini untuk detail lebih lanjut.

— Justin
sumber