multicub {CUB}

Joint plot of estimated CUB models in the parameter space

Description

Return a plot of estimated CUB models represented as points in the parameter space.

Usage

multicub(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 CUB 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. Depending on csiplot and paiplot and on desired output, x and y coordinates may be set to \pi and \xi, respectively.

Examples

data(univer)
listord<-univer[,8:12]
multicub(listord,colours=rep("red",5),cex=c(0.4,0.6,0.8,1,1.2),
  pch=c(1,2,3,4,5),xlim=c(0,0.4),ylim=c(0.75,1),pos=c(1,3,3,3,3))
###############################
m1<-5; m2<-7;  m3<-9
pai<-0.7;csi<-0.6
n1<-1000; n2<-500; n3<-1500
ord1<-simcub(n1,m1,pai,csi)
ord2<-simcub(n2,m2,pai,csi)
ord3<-simcub(n3,m3,pai,csi)
listord<-list(ord1,ord2,ord3)
multicub(listord,labels=c("m=5","m=7","m=9"),pos=c(3,1,4))