Simulation routine for CUSH models
Generate pseudo-random observations following the distribution of a CUSH model without covariates.
simcush(n,m,delta,shelter)
n
: Number of simulated observationsm
: Number of ordinal categoriesdelta
: Shelter parametershelter
: Category corresponding to the shelter choicen<-200 m<-7 delta<-0.3 shelter<-3 simulation<-simcush(n,m,delta,shelter) plot(table(simulation),xlab="Ordinal categories",ylab="Frequencies")
probcush
Useful links