Ich versuche, Zähldaten in R zu modellieren , die scheinbar unterdispers sind (Dispersionsparameter ~ .40). Dies ist wahrscheinlich der Grund, warum ein glm
with- family = poisson
oder ein negatives binomial ( glm.nb
) -Modell nicht signifikant sind. Wenn ich mir die Beschreibungen meiner Daten anschaue, habe ich nicht den typischen Versatz der Zähldaten und die Residuen in meinen beiden Versuchsbedingungen sind auch homogen.
Meine Fragen sind also:
Muss ich für meine Zähldaten sogar spezielle Regressionsanalysen verwenden, wenn sich meine Zähldaten nicht wirklich wie Zähldaten verhalten? Manchmal stelle ich eine Nichtnormalität fest (normalerweise aufgrund der Kurtosis), aber ich habe die Perzentil-Bootstrap-Methode zum Vergleichen der Mittelwerte verwendet (Wilcox, 2012), um die Nichtnormalität zu berücksichtigen. Können Methoden zur Zählung von Daten durch eine von Wilcox vorgeschlagene und im WRS-Paket umgesetzte robuste Methode ersetzt werden?
Wie kann ich die Unterdispersion berücksichtigen, wenn ich Regressionsanalysen für Zählungsdaten verwenden muss? Das Poisson und die negative Binomialverteilung nehmen eine höhere Streuung an, also sollte das nicht angemessen sein, oder? Ich habe darüber nachgedacht, die Quasi-Poisson- Verteilung anzuwenden , aber dies wird normalerweise für eine Überdispersion empfohlen. Ich habe etwas über Beta-Binomial- Modelle gelesen , die sowohl Über- als auch Unterdispersionen zu berücksichtigen
VGAM
scheinen. Die Autoren scheinen jedoch eine getitelte Poisson-Distribution zu empfehlen , aber ich kann sie nicht im Paket finden .
Kann jemand eine Prozedur für unterdisperse Daten empfehlen und vielleicht ein Beispiel für R-Code dafür bereitstellen?