Multivariate Poisson regression with covariance structure
Understanding the Metropolis-Hastings algorithm