if (FALSE) {
#############################################################################
# EXAMPLE 1: data.read | 1PL model in TAM
#############################################################################
data(data.read, package="sirt")
dat <- data.read
#*** Model 1: 1PL model
mod1 <- TAM::tam.mml( resp=dat )
summary(mod1)
# item fit statistics
imod1 <- CDM::IRT.RMSD(mod1)
summary(imod1)
#############################################################################
# EXAMPLE 2: data.math| RMSD and MD statistic for assessing DIF
#############################################################################
data(data.math, package="sirt")
dat <- data.math$data
items <- grep("M[A-Z]", colnames(dat), value=TRUE )
#-- fit multiple group Rasch model
mod <- TAM::tam.mml( dat[,items], group=dat$female )
summary(mod)
#-- fit statistics
rmod <- CDM::IRT.RMSD(mod)
summary(rmod)
#############################################################################
# EXAMPLE 3: RMSD statistic DINA model
#############################################################################
data(sim.dina)
data(sim.qmatrix)
dat <- sim.dina
Q <- sim.qmatrix
#-- fit DINA model
mod1 <- CDM::gdina( dat, q.matrix=Q, rule="DINA" )
summary(mod1)
#-- compute RMSD fit statistic
rmod1 <- CDM::IRT.RMSD(mod1)
summary(rmod1)
}
Run the code above in your browser using DataLab