CUB (version 0.1)

bitgama: Shifted Binomial distribution with covariates

Description

Return the shifted Binomial probabilities of ordinal responses where the feeling component is explained by covariates via a logistic link.

Usage

bitgama(m, ordinal, W, gama)

Arguments

m

Number of ordinal categories

ordinal

Vector of ordinal responses

W

Matrix of covariates for the feeling component

gama

Vector of parameters for the feeling component, with length equal to NCOL(W)+1 to account for an intercept term (first entry of gama)

Value

A vector of the same length as ordinal, where each entry is the shifted Binomial probability of the corresponding observation and feeling parameter

See Also

logis, probcub0q, probcubpq

Examples

n<-100
m<-7
W<-sample(c(0,1),n,replace=TRUE)
gama<-c(0.2,-0.2)
csivett<-logis(W,gama)
ordinal<-rbinom(n,m-1,csivett)+1
pr<-bitgama(m,ordinal,W,gama)