powered by
scores.gpcm.bug(item.scores,sims, mdl = c("rasch.bug", "pcm.bug", "tpl.bug","gpcm.bug"), gibbs = c("bugs","jags"))
scores