powered by
probihg(m, theta)
IHG
m<-10 theta<-0.30 pr<-probihg(m, theta) plot(1:m,pr,type="h",xlab="Ordinal categories") points(1:m,pr,pch=19)
Run the code above in your browser using DataLab