Now let us fit the DINA model to the data and estimate attribute profiles using three different approaches
Code
library(GDINA)est <-GDINA(Y, Q15, model ="DINA", verbose =0)MLE <-personparm(est, what ="MLE")MAP <-personparm(est, what ="MAP")EAP <-personparm(est, what ="EAP")