Thanks for your reply. I am trying to get fit statistics with the M2 command. M2 requires complete dataset with no missing. Therefore, I tried the following command but failed:
> M2(modeling.irt,impute=10,QMC=TRUE)
Error in rep.int(rep.int(seq_len(nx), rep.int(rep.fac, nx)), orep) :
vector is too large
In addition: Warning messages:
1: High-dimensional models should use quasi-Monte Carlo integration. Pass QMC=TRUE
2: High-dimensional models should use quasi-Monte Carlo integration. Pass QMC=TRUE
3: High-dimensional models should use quasi-Monte Carlo integration. Pass QMC=TRUE
I also tried imputation but it was said that the information matrix was not positive definite:
> careful.scores <- mirt::fscores(modeling.irt, method='MAP', MI=100)
Error: Information matrix is not positive definite
Each of my items has 5 ordinal categories (1=Strongly Disagree; 2 = Disagree; 3= Neutral; 4 = Agree; 5 = Strongly Agree). Some items have very sparse response in the Strongly Disagree category (less than 5), even that my sample size is close to 500. The number of items is 55. I wonder if it helps to collapse categories into 3 categories (e.g., 1 = Disagree; 2 = Neutral; 3 = Agree)?
Thanks again!
Chester Kam