model {

for (i in 1:n) {
  Y[i,1:2] ~ dmnorm(mu.Y[],Prec.Y[,]);
}

mu.Y[1:2] ~ dmnorm(mean.mu[],prec.mu[,])
for (k in 1:q) {
  mean.mu[k] <- 0.0;
  prec.mu[k,k] <- 0.001;
  for (l in (k+1):q) {
    prec.mu[l,k] <- 0.0;
    prec.mu[k,l] <- prec.mu[l,k]
  }
}

Prec.Y[1:q,1:q] ~ dwish(B[,],3)
for (k in 1:q) {
  B[k,k] <- 1;
  for (l in (k+1):q) {
    B[k,l] <- 0;
    B[l,k] <- B[k,l];
  }
}

}