Ошибка в r.squaredGLMM()

Я создаю 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))

есть идеи?


person mtao    schedule 02.07.2016    source источник
comment
Спасибо за всю работу, чтобы привести ваш вопрос в форму, добавив воспроизводимый пример. К сожалению, я не знаю, как решить эту ошибку, но теперь она воспроизводима на основе информации в вопросе, поэтому, надеюсь, кто-то еще в сообществе сможет помочь.   -  person josliber♦    schedule 03.07.2016
comment
На R-forge теперь исправлено: обновление до версии 1.15.8   -  person Kamil Bartoń    schedule 05.07.2016
comment
@KamilBartoń, можете ли вы опубликовать свой комментарий в качестве ответа? (Просто очищаю старые lme4 вопросы...)   -  person Ben Bolker    schedule 16.10.2016
comment
@BenBolker сделал. Как только новая версия появится в CRAN, этот вопрос можно будет удалить.   -  person Kamil Bartoń    schedule 16.10.2016
comment
можно оставить на случай, если люди используют устаревшие версии (такое бывает...)   -  person Ben Bolker    schedule 16.10.2016


Ответы (1)


Это была ошибка, которая была исправлено в версии >= 1.15.8 (скоро в CRAN, сейчас в R-Forge).

person Kamil Bartoń    schedule 16.10.2016