dat <- expand.table(LSAT6)
head(dat)
model <- mirt.model('F = 1-5
CONSTRAIN = (1-5, a1)')
(mod <- mirt(dat, model))
coef(mod)
#equivalentely, but with a different parameterization
mod2 <- mirt(dat, 1, itemtype = 'Rasch')
anova(mod, mod2) #equal
coef(mod2)
sqrt(coef(mod2)$GroupPars[2]) #latent SD equal to the slope in mod
Run the code above in your browser using DataLab