data(data.ecpe)
# Model 1: LCDM model
mod1 <- gdina( data.ecpe$data[,-1], q.matrix= data.ecpe$q.matrix , link="logit")
summary(mod1)
# Model 2: DINA model
mod2 <- gdina( data.ecpe$data[,-1], q.matrix= data.ecpe$q.matrix , rule="DINA")
summary(mod2)
# Model comparison using likelihood ratio test
anova(mod1,mod2)
## Model loglike Deviance Npars AIC BIC Chisq df p
## 2 Model 2 -42841.61 85683.23 63 85809.23 86185.97 206.0359 18 0
## 1 Model 1 -42738.60 85477.19 81 85639.19 86123.57 NA NA NA
Run the code above in your browser using DataLab