Bayesian hierarchical spatial models: Implementing the Besag York Mollié model in stan