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.

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

`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) |

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)`

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? Tweet to @rdrrHQ 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.