small fix: compute the ergodic distribution at the estimates before calling markov()