CUB (version 0.1)

betabinomialcsi: Beta-Binomial probabilities of ordinal responses, given feeling parameter for each observation

Description

Compute the Beta-Binomial probabilities of the given ordinal responses, with feeling parameter specified for each observation, and with the same overdispersion parameter for all the responses.

Usage

betabinomialcsi(m, ordinal, csivett, phi)

Arguments

m

Number of ordinal categories

ordinal

Vector of ordinal responses

csivett

Vector of feeling parameters of the Beta-Binomial distribution for the given ordinal responses

phi

Overdispersion parameter of the Beta-Binomial distribution

Value

A vector of the same length as ordinal, containing the Beta-Binomial probability of each observation for the corresponding feeling parameter and for the specified overdispersion parameter

References

Iannario, M. (2014). Modelling Uncertainty and Overdispersion in Ordinal Data, Communications in Statistics - Theory and Methods, 43, 771--786

See Also

betar, betabinomial

Examples

data(relgoods)
m<-10
ordinal<-relgoods[,37]
age<-2014-relgoods[,4]
lage<-log(age)-mean(log(age))
nona<-na.omit(cbind(ordinal,lage))
ordinal<-nona[,1]
W<-nona[,2]
gama<-c(-0.61,-0.31)
phi<-0.16 
csivett<-logis(W,gama)
pr<-betabinomialcsi(m,ordinal,csivett,phi)