### Real data example: Asahi-Todai survey (not run)
## Collapses 5-category ordinal survey items into 3 categories for estimation
data(AsahiTodai)
out.varinf <- ordIRT(.rc = AsahiTodai$dat.all, .starts = AsahiTodai$start.values,
.priors = AsahiTodai$priors, .D = 1,
.control = {list(verbose = TRUE,
thresh = 1e-6, maxit = 500)})
## Compare against MCMC estimates using 3 and 5 categories
cor(ideal3, out.varinf$means$x)
cor(ideal5, out.varinf$means$x)
Run the code above in your browser using DataLab