CUB (version 0.1)

ellecub: Log-likelihood function of a CUB model without covariates

Description

Compute the log-likelihood function of a CUB model without covariates fitting ordinal responses, possibly with different vector of parameters for each observation.

Usage

ellecub(m, ordinal, assepai, assecsi)

Arguments

m

Number of ordinal categories

ordinal

Vector of ordinal responses

assepai

Vector of uncertainty parameters for given observations (with the same length as ordinal)

assecsi

Vector of feeling parameters for given observations (with the same length as ordinal)

See Also

loglikCUB

Examples

m<-7
n0<-230
n1<-270
bet<-c(-1.5,1.2)
gama<-c(0.5,-1.2)
pai0<-1/(1+exp(-bet[1])); csi0<-1/(1+exp(-gama[1]))
pai1<-1/(1+exp(-sum(bet))); csi1<-1/(1+exp(-sum(gama)))
ordinal0<-simcub(n0,m,pai0,csi0)
ordinal1<-simcub(n1,m,pai1,csi1)
ordinal<-c(ordinal0,ordinal1)
assepai<-c(rep(pai0,n0),rep(pai1,n1))
assecsi<-c(rep(csi0,n0),rep(csi1,n1))
lli<-ellecub(m,ordinal,assepai,assecsi)