An R routine (sim150_bugs.R) reads the data (tsim150.txt), calls bugs commands and then plot the bayesian confidence intervals for ordinal response variable using cumulative logit model.

The following files are needed to execute the BUGS software:
sim150_bugs.bug
to run the ordinal logistic model;
sim150_bugs.cmd
for burn-in and to get summary iterations;
sim150_bugs.in
have the initial values for parameters to run bayesian model.