3
Ich baue GLMMs (mit glmer() von "lme4" R-Paket) und manchmal bekomme ich einen Fehler bei der Schätzung von R2-Werten (mit r.squaredGLMM() aus "MuMIn" -Paket).Fehler in r.squaredGLMM()
Das Modell I ist simmilar zu diesem zu passen versuchen:
library(lme4)
lmA <- glmer(x~y+(1|w)+(1|w/k), data = data1, family = binomial(link="logit"))
Dann estime R2, verwende ich:
library(MuMIn)
r.squaredGLMM(lmA)
Und ich bekomme diese:
The result is correct only if all data used by the model has not changed since model was fitted. Error in .rsqGLMM(fam = family(x),
varFx = var(fxpred), varRe = varRe, : 'names' attribute [2] must be the same length as the vector [0]
Haben Sie eine Idee, warum dieser Fehler auftritt? Wenn ich zum Beispiel nur einen zufälligen Faktor verwende (in diesem Fall (1|w)), wird dieser Fehler nicht angezeigt.
Hier ist meine Daten-Set:
data1 <-
structure(list(w = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L,
1L, 2L, 1L), .Label = c("CA", "CB"), class = "factor"), k = structure(c(4L,
4L, 3L, 3L, 3L, 4L, 1L, 3L, 2L, 3L, 2L), .Label = c("CAF01-CAM01",
"CAM01", "CBF01-CBM01", "CBM01"), class = "factor"), x = c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L), y = c(-0.034973549,
0.671720643, 4.557044729, 5.347170897, 2.634240583, -0.555740207,
4.118277809, 2.599825716, 0.95853864, 4.327804344, 0.057331718
)), .Names = c("w", "k", "x", "y"), class = "data.frame", row.names = c(NA,
-11L))
Irgendwelche Gedanken?
Vielen Dank für die Arbeit, um Ihre Frage in Form zu bringen, indem Sie ein reproduzierbares Beispiel hinzufügen. Leider weiß ich nicht, wie ich diesen Fehler beheben kann, aber er ist nun reproduzierbar, basierend auf den Informationen in der Frage, so dass hoffentlich jemand anderer in der Community helfen kann. – josliber
Es ist jetzt auf R-Schmiede behoben: Update auf Version 1.15.8 –
@ KamilBartoń, können Sie Ihren Kommentar als Antwort posten? (Nur alte "lme4" Fragen aufräumen ...) –