Modelling heterogeneity in discrete spatial data models via MAP and MCMC methods