Efficient Markov chain Monte Carlo sampling for hierarchical hidden Markov models