# \donttest{
data("cesc")
fit<-ordinalbayes(Stage~1, data=cesc, x=cesc[,5:45],
model="regressvi",gamma.ind="fixed",
pi.fixed=0.99, adaptSteps=1000, burnInSteps=1000, nChains=2,
numSavedSteps=2000, thinSteps=2, seed=26)
plot(fit)
# }
Run the code above in your browser using DataLab