Poisson log-linear models are ubiquitous in many applications, and one of the most popular approaches for parametric count regression. In Bayesian context, however, there no sufficient specific computational tools efficient sampling from posterior distribution parameters, standard algorithms, such as random walk Metropolis-Hastings or Hamiltonian Monte Carlo typically used. Herein, we developed...