I have the same problem with Hend. I need to compare the ability of different dimensions in a three dimension mirt model.
library("mirt")
setwd('~/My_R')
# data
D27 <- read.csv("all27.csv")
# the confirmatory model:
mod1.mixed27<-mirt.model
('
F1 = 1,6,9,13,14,17,19,21,23
F2 = 2,3,4,7,8,10,15,16,18,20,24,25,27
F3 = 5,11,12,22,26
START = (14,a1,1.0),(2,a2,1.0),(5,a3,1.0)
FIXED = (14,a1),(2,a2),(5,a3)
COV = F1*F2*F3
')
# data types
mixed.types <- rep("gpcm", 27)
mixed.types[3] <- "2PL"
mixed.types[11] <- "2PL"
mixed.types[21] <- "2PL"
mod1.27 <- mirt(D27, mod1.mixed27, itemtype = mixed.types, method = 'MHRM')
coef1.27.sim <- coef(mod1.27, simplify=TRUE)
write.csv(coef1.27.sim$items, file="coef1_27_sim__iterms.csv")
mod1.27.abil <- fscores(mod1.27, QMC = TRUE)
write.csv(mod1.27.abil, file="mod1_27_abil.csv")