Я попробовал несколько предложений здесь и на других сайтах, но ни одно из них не работает.
Моя модель выглядит так
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
Если бы кто-нибудь мог предложить, как я могу решить эту проблему, это было бы очень признательно
Привет, Дилан
REform=NA
, а неREform=NULL
? - person Ben Bolker   schedule 01.03.2014