library(rjags)
mvn.data <- list(q=2,n=1)

m <- jags.model("/Users/dbmad/Documents/G6102/jags/mvn.bug", 
                data=mvn.data,								
                nchain=1
               )
m$update(1000)
parameters <- c("Prec.Y")
x <- jags.samples(m,parameters,n.iter=10000)

samples <- coda.samples(m, parameters, 10000)
summary(samples)
foo<-samples[[1]]

par(mfrow=c(2,2))
plot(density(foo[,1]),main="1,1")
plot(density(foo[,2]),main="2,1")
plot(density(foo[,3]),main="1,2")
plot(density(foo[,4]),main="2,2")

#foo<-summary(x$Prec.Y,mean)$stat
#faa<-summary(x$Prec.Y,sd)$stat