Modelling Heterogeneity in Discrete Spatial Data Models via MAP and MCMC Methods