if (FALSE) {
mod1 <- GDINA(realdata_Tatsuoka1990$dat,realdata_Tatsuoka1990$Q,model="DINA")
mod1
summary(mod1)
# Higher order model
mod2 <- GDINA(realdata_Tatsuoka1990$dat,
realdata_Tatsuoka1990$Q,model="DINA",
att.dist="higher.order")
mod2
anova(mod1,mod2)
}
Run the code above in your browser using DataLab