Saya sedang mengerjakan model data hitung nol-naik menggunakan pscl
paket. 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 glm
with family=poisson
atau 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.
8
Jawaban:
Model Poisson satu-inflasi untuk hitunganYi adalah
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.μi π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.
sumber
score function
danhessian matrix
. Bisakah Anda merekomendasikan saya teks / artikel yang dapat membantu saya untuk belajar lebih banyak tentang itu?zeroinfl
kode harus melakukannya - mengubah kemungkinan Poisson & skor nol untuk mencocokkan model di atas (atau coba hanya mengubah kemungkinan & tidak melewati skor keoptim
). Tentu saja Anda juga bisa bertanya di sini atau di SO jika sesuai untuk referensi atau membantu hal-hal yang membuat Anda macet.Paket R
VGAM
memiliki fungsivglm
yang dapat digunakan untuk memenuhi semua jenis model Poisson-esque. Anda dapat menggunakannya untuk menentukan model satu-inflasi, jadi kira-kira seperti ituvglm(Y~X,family=oipospoisson(),data=data)
. Lihat di sini untuk detail lebih lanjut.sumber