Bayesian bandwidth selection in discrete multivariate associated kernel estimators for probability mass functions