A comparison of a pseudo-likelihood estimation and full information maximum likelihood estimation for fitting multivariate longitudinal ordinal data