to jest mój pierwszy post. Jestem naprawdę wdzięczny za tę społeczność.
Usiłuję analizować dane zliczania wzdłużnego, które są obcinane przez zero (prawdopodobieństwo, że zmienna odpowiedzi = 0 wynosi 0), a średnia! = Wariancja, więc wybrano ujemny rozkład dwumianowy nad poissonem.
Funkcje / polecenia, które wykluczyłem:
R
- Funkcja gee () w R nie uwzględnia skracania zera ani ujemnego rozkładu dwumianowego (nawet przy załadowanym pakiecie MASS)
- glm.nb () w R nie pozwala na różne struktury korelacji
- vglm () z pakietu VGAM może korzystać z rodziny posnegbinomial, ale ma ten sam problem, co polecenie ztnb Staty (patrz poniżej), ponieważ nie mogę ponownie modelować przy użyciu niezależnej struktury korelacji.
Stata
- Gdyby dane nie były podłużne, mógłbym po prostu użyć pakietów Stata ztnb do uruchomienia mojej analizy, ALE to polecenie zakłada, że moje obserwacje są niezależne.
Wykluczyłem również GLMM z różnych przyczyn metodologicznych / filozoficznych.
Na razie zdecydowałem się na komendę Stata xtgee (tak, wiem, że xtnbreg robi to samo), która uwzględnia zarówno niezależne struktury korelacji, jak i ujemną dwumianową rodzinę, ale nie obcinanie zera. Dodatkową korzyścią wynikającą z używania xtgee jest to, że mogę również obliczyć wartości qic (za pomocą polecenia qic) w celu ustalenia najlepiej dopasowanych struktur korelacji dla moich zmiennych odpowiedzi.
Jeśli istnieje pakiet / polecenie w R lub Stata, które mogą wziąć pod uwagę 1) rodzinę dwumianową, 2) GEE i 3) obcinanie zera, chciałbym wiedzieć.
Byłbym bardzo wdzięczny za wszelkie twoje pomysły. Dziękuję Ci.
-Casey