## Not run:
# #############################################################################
# # EXAMPLE 1: Expected counts gdm function
# #############################################################################
#
# data( data.fraction1 )
# dat <- data.fraction1$data
# theta.k <- seq( -6 , 6 , len=11 ) # discretized ability
#
# #--- Model 1: Rasch model
# mod1 <- gdm( dat , irtmodel="1PL" , theta.k=theta.k , skillspace="normal" ,
# centered.latent=TRUE )
# emod1 <- IRT.expectedCounts(mod1)
# str(emod1)
#
# #############################################################################
# # EXAMPLE 2: Expected counts gdina function
# #############################################################################
#
# data(sim.dina)
# data(sim.qmatrix)
#
# #--- Model 1: estimation of the GDINA model
# mod1 <- gdina( data = sim.dina , q.matrix = sim.qmatrix)
# summary(mod1)
# emod1 <- IRT.expectedCounts( mod1 )
# str(emod1)
#
# #--- Model 2: GDINA model with two groups
# mod2 <- gdina( data = sim.dina , q.matrix = sim.qmatrix , group = rep(1:2, each=200) )
# summary(mod2)
# emod2 <- IRT.expectedCounts( mod2 )
# str(emod2)
# ## End(Not run)
Run the code above in your browser using DataLab