R 零膨胀负二项式非整数误差

我正在使用

R2Winbug

将ZINB安装到WinBugs中，这是WinBugs模型代码：

Model{
  for(i in 1:n){
    injections[i] ~ dnegbin(p[i],r)
    p[i] <- r/r+ramda[i]*(1-u[i])
    u[i] ~ dbern(p0[i])
    logit(p0[i])  <- alpha0 + alpha1*hiv[i] + alpha2*parteners[i] + 
                     alpha3*married[i] + alpha4*divorced[i] + alpha5*mage[i]
    log(ramda[i]) <- beta0 + beta1*hiv[i] + beta2*parteners[i] + 
                     beta3*married[i] + beta4*divorced[i] + beta5*mage[i]
  }
  # priors for missing covariates
  for(i in 1:n){
    parteners[i] ~ dpois(mu)
    hiv[i] ~ dbin(q,1)
  }

  # priors
  q ~ dunif(0.0,1.0)
  r <- exp(logalpha)
  logalpha ~ dnorm(0,1.0E-8)
  beta0 ~ dnorm(0,1.0E-8); beta1 ~ dnorm(0,1.0E-8) 
  beta2 ~ dnorm(0,1.0E-8); beta3 ~ dnorm(0,1.0E-8)
  beta4 ~ dnorm(0,1.0E-8); beta5 ~ dnorm(0,1.0E-8) 
  alpha0 ~ dnorm(0,1.0E-8); alpha1 ~ dnorm(0,1.0E-8); alpha2 ~ dnorm(0,1.0E-8)
  alpha3 ~ dnorm(0,1.0E-8); alpha4 ~ dnorm(0,1.0E-8); alpha5 ~ dnorm(0,1.0E-8)
  mu ~ dgamma(0.01, 0.01)
}

模型{
for（1:n中的i）{
注射[i]~dnegbin（p[i]，r）
p[i]你能发布准确的错误消息吗？好吧，二项式家族确实期望值是……二项式的。这似乎是一个应该首先向包维护者提出的问题，只有在失败之后，才可以向其他论坛提出。你能发布准确的错误消息吗？好吧，二项式家族确实期望值是。。。。这似乎是一个应该首先向包维护者提出的问题，只有在这个问题失败之后，才应该向其他论坛提出。