«Ошибка в handleRes (res): NA» из базовой полиномиальной логистической регрессии с использованием BRugs

Question

Я получаю ошибку «Ошибка в handleRes (res): NA» для модели многочленной регрессии в R с использованием BRugs.Приведенный ниже код представляет собой аналогичный пример, который намного проще и выдает ту же ошибку.

Модель была скомпилирована, но она не работает на этапе инициализации.

Кажется, у dmulti () есть некоторые проблемыпоскольку код будет работать, если я изменю результат на одну категориальную переменную и использую dcat ().

library(BRugs) 
library(Hmisc)

n<-100
K<-4
dat<-t(rmultinom(n,2,rep(1/K,K)))
datlist<-list(y=dat,N=n)
bugsData(datlist, file="./test.txt")

modelString = "
model
{# likelihood
for (i in 1:N) {
  denom[i]<-1+exp(beta[1])+exp(beta[2])+exp(beta[3])
  p[i,1]<-exp(beta[1])/denom[i]
  p[i,2]<-exp(beta[2])/denom[i]
  p[i,3]<-exp(beta[3])/denom[i]
  p[i,4]<-1/denom[i]
  y[i,1:4]~dmulti(p[i,1:4],2) 

}
# priors
for(k in 1:3) {
beta[k]~dnorm(0,0.0001)
}
}
"

writeLines(modelString,con="./model.txt")


modelFile="./model.txt"
mydata="./test.txt"
modelCheck(modelFile)
modelData(mydata)
modelCompile(numChains=1)
modelGenInits()

modelCheck (modelFile)

модель синтаксически верна

modelData (mydata)

данные загружены

modelCompile (numChains = 1)

скомпилированная модель

modelGenInits ()

Ошибка в handleRes (res): NA