Simulation routine for IHG models
Generate pseudo-random observations following the given IHG distribution.
simihg(n,m,theta)
n
: Number of simulated observationsm
: Number of ordinal categoriestheta
: Preference parametern<-300 m<-9 theta<-0.4 simulation<-simihg(n,m,theta) plot(table(simulation),xlab="Number of categories",ylab="Frequencies")
probihg
Useful links