Hi all,
thanks to those who responded. Martyn Plummer gave me the solution:
Don't calculate Tau at all and just include the invA matrix directly in
the dmnorm call.
Thanks again,
Simon.
On Wed, 2010-02-17 at 10:36 +1000, Simon Blomberg wrote:
> 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
|