# NOT RUN {
## ordinal data
f <- gjamSimData(S = 15, Q = 3, typeNames = 'OC')
ml <- list(ng = 1500, burnin = 500, typeNames = f$typeNames, holdoutN = 10)
out <- gjam(f$formula, f$xdata, f$ydata, modelList = ml)
# repeat with ng = 2000, burnin = 500, then plot data here:
pl <- list(trueValues = f$trueValues, width=3, height=2)
fit <- gjamPlot(output = out, plotPars = pl)
# }
Run the code above in your browser using DataLab