data(relgoods)
m<-10
ordinal<-relgoods[,29]
gender<-relgoods[,2]
nona<-na.omit(cbind(ordinal,gender))
ordinalnew<-nona[,1]
Y<-nona[,2]
bet<-c(-0.81, 0.93)
csi<-0.20
probi<-probcubp0(m,ordinalnew,Y,bet,csi)
Run the code above in your browser using DataLab