Compute the log-likelihood function for a CUSH model with covariates for the given ordinal responses.
loglikcushcov(m, ordinal, X, omega, shelter)
Number of ordinal categories
Vector of ordinal responses
Matrix of selected covariates for explaining the shelter parameter
Vector of parameters for explaining the shelter effect, with length equal to NCOL(X)+1 to account for an intercept term (first entry of omega)
Category corresponding to the shelter choice