loglikCUB: Log-likelihood function for CUB models

Description Usage Arguments Details See Also Examples

View source: R/loglikCUB.R

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)

CUB documentation built on May 19, 2017, 4:24 p.m.

Search within the CUB package
Search all R packages, documentation and source code