Untuk menganalisis jumlah burung dengan inflasi nol, saya ingin menerapkan model jumlah dengan inflasi nol menggunakan paket R pscl . Namun, setelah melihat contoh yang disediakan dalam dokumentasi untuk salah satu fungsi utama ( ? Zeroinfl ), saya mulai meragukan apa keuntungan sebenarnya dari model ini. Menurut kode sampel yang diberikan di sana, saya menghitung poisson standar, quasi-poisson dan model bionomial negatif, poisson nol sederhana dan model binomial negatif dan poisson nol-meningkat dan model binomial negatif dengan regressor untuk komponen nol. Kemudian saya memeriksa histogram dari data yang diamati dan yang dipasang. (Ini kode untuk mereplikasi itu.)
library(pscl)
data("bioChemists", package = "pscl")
## standard count data models
fm_pois <- glm(art ~ ., data = bioChemists, family = poisson)
fm_qpois <- glm(art ~ ., data = bioChemists, family = quasipoisson)
fm_nb <- glm.nb(art ~ ., data = bioChemists)
## with simple inflation (no regressors for zero component)
fm_zip <- zeroinfl(art ~ . | 1, data = bioChemists)
fm_zinb <- zeroinfl(art ~ . | 1, data = bioChemists, dist = "negbin")
## inflation with regressors
fm_zip2 <- zeroinfl(art ~ fem + mar + kid5 + phd + ment | fem + mar + kid5 + phd +
ment, data = bioChemists)
fm_zinb2 <- zeroinfl(art ~ fem + mar + kid5 + phd + ment | fem + mar + kid5 + phd +
ment, data = bioChemists, dist = "negbin")
## histograms
breaks <- seq(-0.5,20.5,1)
par(mfrow=c(4,2))
hist(bioChemists$art, breaks=breaks)
hist(fitted(fm_pois), breaks=breaks)
hist(fitted(fm_qpois), breaks=breaks)
hist(fitted(fm_nb), breaks=breaks)
hist(fitted(fm_zip), breaks=breaks)
hist(fitted(fm_zinb), breaks=breaks)
hist(fitted(fm_zip2), breaks=breaks)
hist(fitted(fm_zinb2), breaks=breaks)!
Saya tidak dapat melihat perbedaan mendasar antara model-model yang berbeda (selain dari itu contoh data tidak tampak sangat "inflasi nol" bagi saya ...); sebenarnya tidak ada model yang menghasilkan estimasi setengah nol dari jumlah nol. Adakah yang bisa menjelaskan apa keuntungan dari model zero-inflated? Saya kira pasti ada alasan untuk memilih ini sebagai contoh untuk fungsi.
sumber