R/loglikcub0q.R

In CUB: A Class of Mixture Models for Ordinal Data

#' @title Log-likelihood function of a CUB model with covariates for the feeling component
#' @description Compute the log-likelihood function of a CUB model fitting ordinal data, with \eqn{q} 
#' covariates for explaining the feeling component.
#' @aliases loglikcub0q
#' @usage loglikcub0q(m, ordinal, W, pai, gama)
#' @param m Number of ordinal categories
#' @param ordinal Vector of ordinal responses
#' @param W Matrix of selected covariates for explaining the feeling component
#' @param pai Uncertainty parameter
#' @param gama Vector of parameters for the feeling component, with length NCOL(W) + 1 to account for 
#' an intercept term (first entry of gama)
#' @keywords internal
#' 


loglikcub0q <-
function(m,ordinal,W,pai,gama){
  
  if (is.factor(ordinal)){
    ordinal<-unclass(ordinal)
  }
  W<-as.matrix(W)
  probn<-probcub0q(m,ordinal,W,pai,gama)
  sum(log(probn))
}