Dear all,
we would like to calculate the DIC in the following models using Winbugs
(1.4) output if possible. Does anybody know how to do this?
model 1;
{
#priors
M<-100*maxC
lambda ~ dunif(0,M)
p ~ dbeta(1,1)
for( i in 1 : 31 ) {
abundance[i] ~ dpois(lambda)
catch[i] ~ dbin(p,abundance[i])
}
}
model 2;
{
#priors
alpha ~ dunif(0.01,1)
beta~dunif(0.01 ,100)
maxi<-100*maxC
lambda ~ dgamma(0.001,0.001)
for( i in 1 : 31 ) {
abundance[i] ~ dpois(lambda)
p[i] ~ dbeta(alpha,beta)
catch[i] ~ dbin(p[i],abundance[i])
}
}
Thank you very much for any replies.
Regards,
Verena
-------------------------------------------------------------------
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
|