Multivariate Poisson Regression with Full Covariance Structure