bu benim ilk yazım. Bu topluluk için gerçekten minnettarım.
Sıfır kesilmiş (yanıt değişkeni = 0 olması olasılığı) uzunlamasına sayım verilerini ve ortalama! = Varyansı analiz etmeye çalışıyorum, bu nedenle bir poisson üzerinde negatif bir binom dağılımı seçildi.
Dışarıda bıraktığım işlevler / komutlar:
R,
- R'deki gee () işlevi sıfır kesmeyi veya negatif binom dağılımını (MASS paketi yüklü olsa bile) açıklamaz
- R'deki glm.nb () farklı korelasyon yapılarına izin vermez
- VGAM paketinden vglm () posnegbinomial ailesini kullanabilir, ancak bağımsız olmayan bir korelasyon yapısı kullanarak modelleri yeniden yerleştiremediğim için Stata'nın ztnb komutu (aşağıya bakınız) ile aynı soruna sahiptir.
Stata
- Veriler boylamsal değilse, analizimi çalıştırmak için sadece Stata paketleri ztnb'yi kullanabilirim, ancak bu komut gözlemlerimin bağımsız olduğunu varsayar.
Ayrıca çeşitli metodolojik / felsefi nedenlerle GLMM'yi de dışladım.
Şimdilik, Stata'nın xtgee komutuna karar verdim (evet, xtnbreg'in de aynı şeyi yaptığını biliyorum), hem bağımsız olmayan korelasyon yapılarını hem de negatif binom ailesini hesaba katar, ancak sıfır kesmeyi değil. Xtgee kullanmanın ek yararı, yanıt değişkenlerim için en uygun korelasyon yapılarını belirlemek için qic değerlerini de (qic komutunu kullanarak) hesaplayabilmemdir.
R veya Stata'da 1) nbinomiyal aile, 2) GEE ve 3) sıfır kesmeyi dikkate alabilecek bir paket / komut varsa, bilmek için ölmek isterim.
Sahip olabileceğiniz herhangi bir fikri çok takdir ediyorum. Teşekkür ederim.
-Casey