Last week I've been trying to fit some nominal response models (NRM, Bock, 1972) with both MIRT & Mplus. Both do this well, with the MIRT-package having the advantage in producing nice plots of the estimated response probabilities conditional on the latent variable. However, I seem to be on the wrong track in interpreting the ak0, ak1, ak2, ... and d0, d1, d2,... parameter values in MIRT. I was initially under the impression that the models I was fitting had the following form:
So I used the estimated parameter values from MIRT as this:
However, when calculating the probabilities using these parameter values they are quite a bit off from what I see on the MIRT-plots using itemfit and what I calculate with the Mplus parameter values. The MIRT-plots using itemfit and what I calculate with the estimated Mplus parameter values are in agreement though. This makes me suspect that the model used for NRM in MIRT may not be what I thought it was. For there certainly are different ways to define such a model. Unfortunately, I cannot find any documentation of the exact NRM used within MIRT? Is there any information somewhere on the NRM used in MIRT?