Compute the log-likelihood function of a CUB model fitting ordinal data with covariates for explaining both the feeling and the uncertainty components.
loglikcubpq(m, ordinal, Y, W, bet, gama)
Number of ordinal categories
Vector of ordinal responses
Matrix of selected covariates for explaining the uncertainty component
Matrix of selected covariates for explaining the feeling 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)
Vector of parameters for the feeling component, whose length equals NCOL(W) + 1 to account for an intercept term (first entry of gama)