Integrated nested Laplace approximation method for hierarchical Bayesian inference of spatial model with application to functional magnetic resonance imaging data