Compute the log-likelihood function of a CUB model fitting ordinal responses with covariates for explaining the uncertainty component.
loglikcubp0(m, ordinal, Y, bbet, ccsi)
Number of ordinal categories
Vector of ordinal responses
Matrix of selected covariates for explaining the uncertainty component
Vector of parameters for the uncertainty component, with length equal to NCOL(Y)+1 to account for an intercept term (first entry of bbet)
Feeling parameter