On the estimation of the polychoric correlation coefficient via Markov Chain Monte Carlo methods