Title: On model fitting for multivariate polytomous response data