Я пытаюсь запустить следующую модель. В нем есть несколько ссылок на логит, которые находятся в циклах.
model {
for(i in 1:8){
year.eff[i] ~ dnorm(0,tau.yr)
annual.controlrenest[i] <- pow(ilogit(beta0+year.eff[i]+beta.renest),37)
annual.treatmentrenest[i] <-
pow(ilogit(beta0+year.eff[i]+beta.treatment+beta.renest),37)
annual.controlnest[i] <- pow(ilogit(beta0+year.eff[i]),37)
annual.treatmentnest[i] <- pow(ilogit(beta0+year.eff[i]+beta.treatment),37)
annual.controlrenest[i] <- pow(ilogit(beta0+year.eff[i]+beta.renest),37)
annual.treatmentrenest[i] <-
pow(ilogit(beta0+year.eff[i]+beta.treatment+beta.renest),37)
}
# Likelihood
for(i in 1:539){
for(j in (enter.nest[i]+1):left.nest[i]){
logit(phinest[i,j]) <- beta0 + year.eff*year.nest[i] +
beta.treatment*Treatment[i] + beta.renest*Renest[i]
mu[i,j] <- phinest[i,j] * n[i,j-1]
n[i,j] ~ dbern(mu[i,j])
}
}
}
Однако я получаю следующую ошибку и не могу понять, в чем проблема. Я получил эту ошибку раньше, когда моя ссылка на логит не была в цикле (вы не можете векторизовать ссылку на логит в JAGS).
Error in checkForRemoteErrors(lapply(cl, recvResult)) :
3 nodes produced errors; first error: RUNTIME ERROR:
Invalid vector argument to ilogit