JAGS: неверный векторный аргумент для ilogit - PullRequest
0 голосов
/ 03 мая 2018

Я пытаюсь запустить следующую модель. В нем есть несколько ссылок на логит, которые находятся в циклах.

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
...