Я создаю GLMM (используя glmer() из пакета R "lme4"), и иногда я получаю сообщение об ошибке при оценке значений R2 (используя r.squaredGLMM() из пакета "MuMIn").
Модель, которую я пытаюсь подогнать, похожа на эту:
library(lme4)
lmA <- glmer(x~y+(1|w)+(1|w/k), data = data1, family = binomial(link="logit"))
Затем, чтобы оценить R2, я использую:
library(MuMIn)
r.squaredGLMM(lmA)
И я получаю это:
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]
У вас есть идеи, почему появляется эта ошибка? Например, если я использую только один случайный фактор (в данном случае (1|w)
), эта ошибка не появляется.
Вот мой набор данных:
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))
есть идеи?
lme4
вопросы...) - person Ben Bolker   schedule 16.10.2016