CUB (version 1.1.5)

multicube: Joint plot of estimated CUBE models in the parameter space

Description

Return a plot of estimated CUBE models represented as points in the parameter space, where the overdispersion is labeled.

Usage

multicube(listord,mvett,csiplot=FALSE,paiplot=FALSE,...)

Value

Fit a CUBE model to list elements, and then by default it returns a plot of the estimated \((1-\pi, 1-\xi)\) as points in the parameter space, labeled with the estimated overdispersion. Depending on csiplot and paiplot and on desired output, \(x\) and \(y\)

coordinates may be set to \(\pi\) and \(\xi\), respectively.

Arguments

listord

A data matrix, data frame, or list of vectors of ordinal observations (for variables with different number of observations)

mvett

Vector of number of categories for ordinal variables in listord (optional: if missing, the number of categories is retrieved from data: it is advisable to specify it in case some category has zero frequency)

csiplot

Logical: should \(\xi\) or \(1-\xi\) be the \(y\) coordinate

paiplot

Logical: should \(\pi\) or \(1-\pi\) be the \(x\) coordinate

...

Additional arguments to be passed to plot, text, and GEM

Examples

m1<-5; m2<-7;  m3<-9
pai<-0.7;csi<-0.6;phi=0.1
n1<-1000; n2<-500; n3<-1500
ord1<-simcube(n1,m1,pai,csi,phi)
ord2<-simcube(n2,m2,pai,csi,phi)
ord3<-simcube(n3,m3,pai,csi,phi)
listord<-list(ord1,ord2,ord3)
multicube(listord,labels=c("m=5","m=7","m=9"),pos=c(3,1,4),expinform=TRUE)