Simulation routine for CUBE models
Generate pseudo-random observations following the given CUBE distribution.
simcube(n,m,pai,csi,phi)
n
: Number of simulated observationsm
: Number of ordinal categoriespai
: Uncertainty parametercsi
: Feeling parameterphi
: Overdispersion parametern<-300 m<-9 pai<-0.7 csi<-0.4 phi<-0.1 simulation<-simcube(n,m,pai,csi,phi) plot(table(simulation),xlab="Ordinal categories",ylab="Frequencies")
probcube
Useful links