#############################################################################
# EXAMPLE 1: data.read. The 12 items are arranged in 4 testlets
#############################################################################
data(data.read)
# estimate the Rasch model
mod <- rasch.mml2( data.read)
mod$item
# estmate WLEs
mod.wle <- wle.rasch( dat = data.read , b = mod$item$b )
# Yen's Q3 statistic
mod.q3 <- yen.q3( dat = data.read , theta = mod.wle$theta , b = mod$item$b )
# Yen's Q3 statistic with testlets
items <- colnames(data.read)
testlet.matrix <- cbind( substring( items,1,1) , items )
mod.testletq3 <- testlet.yen.q3( q3.res=mod.q3,testlet.matrix=testlet.matrix)
mod.testletq3
Run the code above in your browser using DataLab