Monte Carlo error estimation for multivariate Markov chains