I calculating the graduation of mortality taxes and I am with problems in
the calculation of the DIC in winbugs. In the model with function indicated
in priori I do not obtain to calculate the DIC, then I do not obtain to
compare models. Any suggestions please? Thanks
model{
alpha[1] ~ dgamma(0.1, 0.1)
beta[1] ~ dgamma(0.1, 0.1)
theta[1] ~ dgamma(alpha, beta)I(0, theta[2])
lambda[1] <- theta[1] * e[1]
d[1] ~ dpois(lambda[1])}
for (i in 2 : 66) {
theta[i] ~ dgamma(alpha, beta)I(theta[i-1], theta[i+1])
lambda[i] <- theta[i] * e[i]
d[i] ~ dpois(lambda[i])
alpha[i] ~ dgamma(0.1, 0.1)
beta[i] ~ dgamma(0.1, 0.1)
}
alpha[66] ~ dgamma(0.1, 0.1)
beta[66] ~ dgamma(0.1, 0.1)
theta[66] ~ dgamma(alpha, beta)I(theta[65],1)
lambda[66] <- theta[66] * e[66]
d[66] ~ dpois(lambda[66])
for (i in 1 : 66) {
lambdarep[i] <- theta[i] * e[i]
drep[i] ~ dpois(lambda[i])
thetapred[i]<-drep[i]/e[i]}}
César Neves - Brazil
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|