#############################################################################
# EXAMPLE 1: Extracting factor scores in the DINA model
#############################################################################
data(sim.dina)
data(sim.qmatrix)
# estimate model with only few iterations
mod1 <- din(sim.dina, sim.qmatrix, maxit=10)
summary(mod1)
# MLE
fsc1a <- IRT.factor.scores(mod1 )
# MAP
fsc1b <- IRT.factor.scores(mod1 , type="MAP")
# EAP
fsc1c <- IRT.factor.scores(mod1 , type="EAP")
# compare classification for skill 1
stats::xtabs( ~ fsc1a[,1] + fsc1b[,1] )
graphics::boxplot( fsc1c[,1] ~ fsc1a[,1] )
Run the code above in your browser using DataLab