Berechnung des ICC für die logistische Regression mit zufälligen Effekten

Ich führe ein logistisches Regressionsmodell in der folgenden Form aus:

lmer(response~1+(1|site), family=binomial, REML = FALSE)

Normalerweise würde ich den ICC aus dem Achsenabschnitt und den Restvarianzen berechnen, aber die Zusammenfassung des Modells enthält keine Restvarianzen. Wie berechne ich das?

logistic lme4-nlme glmm Megan
quelle

Warum berechnen Sie den ICC?

AdamO

Um die Annahme zu überprüfen, dass die normale logistische Regression für diese Daten nicht gültig ist, als Beweis, dass ich GLMM verwenden sollte. Ich habe eine Gleichung gefunden: ICClogit = Achsenabschnittsvarianz ^ 2 / (Achsenabschnittsvarianz ^ 2 + pi ^ 2/3). Scheint das vernünftig?

Megan

Sie verwenden den vollständigen Maximum-Likelihood-Ansatz. Können Sie keinen Likelihood-Ratio-Test mit 1 Freiheitsgrad gegen das Fixed-Effects-Modell durchführen?

AdamO

@ Megan: Du hast recht. In der Praxis haben Zeger et al. (1988) schlägt vor ,

funktioniert besser als

als Restvarianz für logistische Regressionsmodelle, obwohl die zwei sehr nahe sind. Siehe SL Zeger, KY Liang und PS Albert. Modelle für longitudinale Daten: ein verallgemeinerter Schätzgleichungsansatz. Biometrics, 44: 1049–1060, 1988.

(15 / 16)^{2} π^{2} / 3

$(15/16)^2\pi^2/3$

π^{2} / 3

$\pi^2/3$

Randel,

@ Megan: Es ist intercept_variance / (intercept_variance + pi^2/3)- also nicht die Varianz quadrieren .

Wolfgang

Antworten:

Sie können die icc()-Funktion aus dem sjstats-Paket verwenden .

In der Hilfedatei ?sjstats::iccfinden Sie einen Verweis auf die Formel für gemischte Modelle mit Binärantwort:

Wus, Crespi CM, Wong WK. 2012. Vergleich von Methoden zur Schätzung des Intraclass-Korrelationskoeffizienten für binäre Antworten in randomisierten Krebspräventionscluster-Studien. Contempory Clinical Trials 33: 869-880 (doi: 10.1016 / j.cct.2012.05.004)

Die verbleibende Abweichung in der logistischen Regression ist auf festgelegt (pi ^ 2) / 3.

Daniel
quelle

Haben Sie eine Referenz für diese Formel?

Jeanine

Meinst du mich? War Ihr Kommentar anfangs nicht beim OP-Posting?

Daniel

@ Jeanine-ICC-Zitat: R. Moineddin, FI Matheson & RH Glazier (2007). Eine Simulationsstudie zur Stichprobengröße für logistische Regressionsmodelle auf mehreren Ebenen. BMC Medical Research Methodology, 7, 34. doi.org/10.1186/1471-2288-7-34

murrt am