CUB (version 0.1)

loglikIHG: Log-likelihood function for IHG models

Description

Compute the log-likelihood function for IHG models with or without covariates to explain the preference parameter.

Usage

loglikIHG(ordinal,m,param,U=0)

Arguments

ordinal

Vector of ordinal responses

m

Number of ordinal categories

param

Vector of parameters for the specified IHG model

U

Matrix of selected covariates to explain the preference parameter (default: no covariate is included in the model)

Details

If no covariate is included in the model, then "param" is the estimate of the preference parameter (theta), otherwise "param" has length equal to NCOL(U) + 1 to account for an intercept term (first entry)

See Also

IHG

Examples

#### Log-likelihood of an IHG model with no covariate
m<-10; theta<-0.14; n<-300
ordinal<-simihg(n,m,theta)
loglik<-loglikIHG(ordinal,m,param=theta)
##################################
#### Log-likelihood of a IHG model with covariate 
data(relgoods)
m<-10
ordinal<-relgoods[,41]
gender<-relgoods[,9]
nona<-na.omit(cbind(ordinal,gender))
ordinal<-nona[,1]
gender<-nona[,2]
nu<-c(-1.55,-0.11) # first entry: intercept term
loglik<-loglikIHG(ordinal, m, param=nu, U=gender)