# NOT RUN {
dat <- expand.table(deAyala)
mod <- mirt(dat, 1)
# marginal estimate
marginal_rxx(mod)
# }
# NOT RUN {
# empirical estimate (assuming the same prior)
fscores(mod, returnER = TRUE)
# empirical rxx the alternative way, given theta scores and SEs
fs <- fscores(mod, full.scores.SE=TRUE)
head(fs)
empirical_rxx(fs)
# }
Run the code above in your browser using DataLab