data(sim.rasch)
# 1PL estimation
mod1 <- tam.mml(resp=sim.rasch)
# 2PL estimation
mod2 <- tam.mml.2pl(resp=sim.rasch , irtmodel="2PL")
# Model comparison
anova( mod1 , mod2 )
## > anova( mod1 , mod2 )
## Model loglike Deviance Npars AIC BIC Chisq df p
## 1 Model 1 -42077.88 84155.77 41 84237.77 84467.40 54.05078 39 0.05508
## 2 Model 2 -42050.86 84101.72 80 84261.72 84709.79 NA NA NA
Run the code above in your browser using DataLab