model {
  xi[1:2]~dmnorm(u[1:2],latprec[1:2,1:2])

  latprec[1:2,1:2] ~ dwish(V[1:2,1:2], df)
  latcov[1:2,1:2] <- inverse(latprec[1:2,1:2])
  latcor <- latcov[1,2]/(sqrt(latcov[1,1])*sqrt(latcov[2,2]))
  df <- 4
  V[1,1] <- 1*(df-3)
  V[1,2] <- 0*(df-3)
  V[2,1] <- 0*(df-3)
  V[2,2] <- 1*(df-3)
  u[1] <- 0
  u[2] <- 0
} #end of model

