Modelling correlated bivariate ordinal data with CUB margins