Log-likelihood function for CUB models

Description

Compute the log-likelihood value of a CUB model fitting given data, with or without covariates to explain the feeling and uncertainty components, or for extended CUB models with shelter effect.

Usage

1
loglikCUB(ordinal,m,param,Y=0,W=0,X=0,shelter=0)

Arguments

ordinal

Vector of ordinal responses (factor type)

m

Number of ordinal categories

param

Vector of parameters for the specified CUB model

Y

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

W

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

X

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

shelter

Category corresponding to the shelter choice (default: no shelter effect is included in the model)

Details

If no covariate is included in the model, then param has the form (π,ξ). More generally, it has the form (\bold{β,γ)} where, respectively, \bold{β} and \bold{γ} are the vectors of coefficients explaining the uncertainty and the feeling components, with length NCOL(Y)+1 and NCOL(W)+1 to account for an intercept term in the first entry. When shelter effect is considered, param corresponds to the first possibile parameterization and hence should be given as (pai1,pai2,csi)

See Also

logLik

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
## Log-likelihood of a CUB model with no covariate
m<-9; n<-300
pai<-0.6; csi<-0.4
ordinal<-simcub(n,m,pai,csi)
param<-c(pai,csi)
loglikcub<-loglikCUB(ordinal,m,param)
##################################
## Log-likelihood of a CUB model with covariate for uncertainty

data(relgoods)
attach(relgoods)
m<-10
naord<-which(is.na(Physician))
nacov<-which(is.na(Gender))
na<-union(naord,nacov)
ordinal<-Physician[-na]; Y<-Gender[-na]
bbet<-c(-0.81,0.93); ccsi<-0.2
param<-c(bbet,ccsi)
loglikcubp0<-loglikCUB(ordinal,m,param,Y=Y)
#######################
## Log-likelihood of a CUB model with covariate for feeling
data(relgoods)
attach(relgoods)
m<-10
naord<-which(is.na(Physician))
nacov<-which(is.na(Gender))
na<-union(naord,nacov)
ordinal<-Physician[-na]; W<-Gender[-na]
pai<-0.44; gama<-c(-0.91,-0.7)
param<-c(pai,gama)
loglikcub0q<-loglikCUB(ordinal,m,param,W=W)
#######################
## Log-likelihood of a CUB model with covariates for both parameters
data(relgoods)
attach(relgoods)
m<-10
naord<-which(is.na(Walking))
nacovpai<-which(is.na(Gender))
nacovcsi<-which(is.na(Smoking))
na<-union(naord,union(nacovpai,nacovcsi))
ordinal<-Walking[-na]
Y<-Gender[-na]; W<-Smoking[-na]
bet<-c(-0.45,-0.48); gama<-c(-0.55,-0.43)
param<-c(bet,gama)
loglikcubpq<-loglikCUB(ordinal,m,param,Y=Y,W=W)
### Log-likelihood of a CUB model with shelter effect
m<-7; n<-400
pai<-0.7; csi<-0.16; delta<-0.15
shelter<-5
ordinal<-simcubshe(n,m,pai,csi,delta,shelter)
pai1<- pai*(1-delta); pai2<-1-pai1-delta
param<-c(pai1,pai2,csi)
loglik<-loglikCUB(ordinal,m,param,shelter=shelter)
##############
### Log-likelihood of a GeCUB
data(univer)
attach(univer)
shelter<-7
modelgecub<-GEM(Formula(officeho~gender|gender|gender),family="cub",shelter=shelter,maxiter=100)
logLik(modelgecub)
param=rep(0.1,6)
loglik<-loglikCUB(officeho,m=7,param=param,shelter=shelter,Y=gender,W=gender,X=gender)

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.