data(univer)
listord<-as.list(univer[,8:12])
mlist<-as.list(rep(7,length(listord)))
multicub(listord, mlist, labelpoints=as.character(1:length(mlist)),
caption="CUB models", colours=rep("black",length(mlist)), symbols=19,
thickness=1.5, xwidth=c(0,1), ywidth=c(0,1), pos=rep(3,length(mlist)))
###############################
m1<-5; m2<-7; m3<-9
mlist<-list(m1,m2,m3)
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,mlist,labelpoints=c("m=5","m=7","m=9"),pos=c(3,1,4))
Run the code above in your browser using DataLab