multicube function

Joint plot of estimated CUBE models in the parameter space

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

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

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

Returns

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.

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)