A Hamiltonian Monte Carlo EM algorithm for generalized linear mixed models with spatial skew latent variables