Another thing is that the polychor function produces maximum likelihood estimates of the pairwise correlation between two variables, but the overall correlation matrix among K variables may not be positive definite.
To do the Bayesian analogue, you basically need the stuff at
and to estimate the correlation matrix, the latent z variables, and the cutpoints simultaneously. Hopefully, that is enough to get started.