# NOT RUN {
library(mirt)
# loading 100-item medical admission test data sets
data(dataMedical)
# binary data set
dataBin <- dataMedical[, 1:100]
# fit Rasch model with mirt package
fit <- mirt(dataBin, model = 1, itemtype = "Rasch")
# factor scores
theta <- as.vector(fscores(fit))
# difficulty estimates
b <- coef(fit, simplify = TRUE)$items[, "d"]
ggWrightMap(theta, b)
item.names <- paste("Item", 1:20)
ggWrightMap(theta, b, item.names = item.names)
# }
Run the code above in your browser using DataLab