Hi all,
We have a model that compiles and runs in JAGS, but not BUGS. In BUGS,
we get the, "expected multivariate node" error. We've been careful with
the data format, so we don't think that is the problem. The model is
below. We suspect it is something to do with the assignment of Tau:
model {
for (k in 1:Ntree) {
p[k]<-1/Ntree
}
for (i in 1:Ndata){
Y[i,1:Nspec]~dmnorm(mu[1:Nspec],Tau[1:Nspec,1:Nspec])
}
for (j in 1:Nspec) {
mu[j] ~ dnorm(0,1.0E-06)
}
K~dcat(p[])
Tau[1:Nspec,1:Nspec] <- invA[1:Nspec,1:Nspec,K]
}
Example data sets for the model can be found at:
http://www.uq.edu.au/~uqsblomb/jagsdata.txt
http://www.uq.edu.au/~uqsblomb/bugsdata.txt
We are happy to use JAGS rather than BUGS, but we would like end-users
to use BUGS, as it is the standard software for Bayesian MCMC stuff, and
more user-friendly.
Thanks,
Simon.
--
Simon Blomberg, BSc (Hons), PhD, MAppStat.
Lecturer and Consultant Statistician
School of Biological Sciences
The University of Queensland
St. Lucia Queensland 4072
Australia
Room 320 Goddard Building (8)
T: +61 7 3365 2506
http://www.uq.edu.au/~uqsblomb
email: S.Blomberg1_at_uq.edu.au
Policies:
1. I will NOT analyse your data for you.
2. Your deadline is your problem.
Statistics is the grammar of science - Karl Pearson
-------------------------------------------------------------------
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
|