# NOT RUN {
# }
# NOT RUN {
dat <- expand.table(Bock1997)
head(dat)
mod <- mirt(dat, 1, 'nominal')
#reproduce table 3 in Bock (1997)
fs <- round(fscores(mod, verbose = FALSE, full.scores = FALSE)[,c('F1','SE_F1')],2)
fttd <- residuals(mod, type = 'exp')
table <- data.frame(fttd[,-ncol(fttd)], fs)
table
mod <- mirt(dat, 1, 'nominal')
coef(mod)
# }
Run the code above in your browser using DataLab