An approximate maximum likelihood procedure for parameter estimation in multivariate discrete data regression models