data(relgoods)
m<-10
ordinal<-relgoods[,29]
gender<-relgoods[,2]
data<-na.omit(cbind(ordinal,gender))
ordinalnew<-data[,1]
W<-data[,2]
pai<-0.44
gama<-c(-0.91,-0.7)
pr<-probcub0q(m,ordinalnew,W,pai,gama)
Run the code above in your browser using DataLab