dat <- expand.table(deAyala)
model <- mirt.model('F = 1-5
CONSTRAIN = (1-5, a1)')
mod <- mirt(dat, model)
Theta <- matrix(seq(-6,6,.01))
tscore <- expected.test(mod, Theta)
tail(cbind(Theta, tscore))
Run the code above in your browser using DataLab