mod <- mirt(Science, 1)
extr.2 <- extract.item(mod, 2)
Theta <- matrix(seq(-6,6, length.out=200))
expected <- expected.item(extr.2, Theta, min(Science[,1])) #min() of first item
head(data.frame(expected, Theta=Theta))
expected.item(extr.2, Theta, min(Science[,1]), include.var=TRUE)
Run the code above in your browser using DataLab