ttr.wish.speciesv.R

  # Setting up the precision matrices
  dof ~ dgamma(1.0, 1.0)
  for (m in 1:ns)
  {
    omega[m,1:nt,1:nt] ~ dwish(ttrI, nt - 1 + dof)
    # Conversion of precision matrix into var-cov matrix
    sigma[m,1:nt,1:nt] <- inverse(omega[m,1:nt,1:nt])
  }
  # monitor # sigma, dof