if (FALSE) {
dat <- expand.table(Bock1997)
head(dat)
itemstats(dat, use_ts=FALSE)
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