Simulation routine for CUB models with shelter effect
Generate pseudo-random observations following the given CUB distribution with shelter effect.
simcubshe(n,m,pai,csi,delta,shelter)
n
: Number of simulated observationsm
: Number of ordinal categoriespai
: Uncertainty parametercsi
: Feeling parameterdelta
: Shelter parametershelter
: Category corresponding to the shelter choicen<-300 m<-9 pai<-0.7 csi<-0.3 delta<-0.2 shelter<-3 simulation<-simcubshe(n,m,pai,csi,delta,shelter) plot(table(simulation),xlab="Ordinal categories",ylab="Frequencies")
probcubshe1
, probcubshe2
, probcubshe3
Useful links