if (FALSE) {
dat <- expand.table(deAyala)
itemstats(dat)
mod <- mirt(dat)
theta_se <- fscores(mod, full.scores.SE = TRUE)
empirical_rxx(theta_se)
theta_se <- fscores(mod, full.scores.SE = TRUE, method = 'ML')
empirical_rxx(theta_se)
empirical_rxx(theta_se, T_as_X = TRUE)
}
Run the code above in your browser using DataLab