Как получить и подобрать прогнозы от glmer?

Я попробовал несколько предложений здесь и на других сайтах, но ни одно из них не работает.

Моя модель выглядит так

model = glmer(cbind(live, dead)~cont1*cont2+ (1|random),
                family = binomial, data = all) 

cont1 и cont2 — непрерывные независимые переменные, random — случайный фактор примерно с 8 уровнями.

Для стандартного glm (т.е.)

model = glm(cbind(live, dead)~cont1*cont2, 
                   family = binomial, data = all) 

Я использовал этот код для получения прогнозов

plot(-100,-100,xlim=c(0,30),ylim=c(0,100),xlab='cont1',
   ylab='alive(%)',col.axis='black',col.lab='black')
   preds = predict(model,data.frame(cont1=0:30,cont2=0),type='response')*100
   lines(0:30,preds,lty=5,lwd=3,col='magenta')   
   preds = predict(model,data.frame(cont1=0:30,cont2=-0.25),type='response')*100
   lines(0:30,preds,lty=3,lwd=3,col='purple')

и т. д.

Мне удалось получить прогнозы с помощью glmer, но я не могу получить прогнозы для каждого уровня cont2, например, в стандартном glm.

Я попытался скопировать код, предложенный здесь cbind количество данных)

predframe <- data.frame(cont1=(all$cont1), cont2=(all$cont2))
predframe$AlivePercent= predict(m1, newdata= predframe,type="response",REform=NULL)

Однако я получаю сообщение об ошибке

Error in `[[<-.data.frame`(`*tmp*`, i, value = integer(0)) : 
  replacement has 0 rows, data has 224

Если бы кто-нибудь мог предложить, как я могу решить эту проблему, это было бы очень признательно

Привет, Дилан


person user3358662    schedule 28.02.2014    source источник
comment
Пожалуйста, опубликуйте минимальный воспроизводимый пример.   -  person fabians    schedule 28.02.2014
comment
дикая догадка: вы хотите REform=NA, а не REform=NULL?   -  person Ben Bolker    schedule 01.03.2014