Published in Journal of Modern Applied Statistical Methods, 2009, Vol. 8, No. 2, pp. 646-658.


A Fortran 77 subroutine is provided for implementing the Gibbs sampling procedure to a multiunidimensional IRT model for binary item response data with the choice of uniform and normal prior distributions for item parameters. In addition to posterior estimates of the model parameters and their Monte Carlo standard errors, the algorithm also estimates the correlations between distinct latent traits.The subroutine requires the user to have access to the IMSL library.