Improved multivariate prediction regions for Markov process models