On bivariate Poisson regression models