multicube {CUB}

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,...)

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

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.

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)