Parallel Markov chain Monte Carlo for Bayesian hierarchical models with big data, in two stages