r - Zero-inflated negative binomial non-integer error -
i fitting zinb in winbugs using r2winbug
, winbugs model code:
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 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) }
now after calling bugs function, model said correct, data loaded correctly, when initials given, says value of binomial hiv[16] must integer! how can correct problem?
Comments
Post a Comment