set.seed(123)
library(Qval)
# \donttest{
K <- 4
I <- 20
IQ <- list(
P0 = runif(I, 0.2, 0.4),
P1 = runif(I, 0.6, 0.8)
)
################################################################
# Example 1: sim.data object #
################################################################
Q <- sim.Q(K, I)
data.obj <- sim.data(Q = Q, N = 500, IQ = IQ,
model = "GDINA", distribute = "horder")
plot(data.obj)
################################################################
# Example 2: CDM object #
################################################################
CDM.obj <- CDM(data.obj$dat, Q, model = "GDINA",
method = "EM", maxitr = 2000, verbose = 1)
plot(CDM.obj)
################################################################
# Example 3: validation object (Hull plot) #
################################################################
MQ <- sim.MQ(Q, 0.1)
CDM.obj <- CDM(data.obj$dat, MQ)
############### ESA ###############
Hull.obj <- validation(data.obj$dat, MQ, CDM.obj,
method = "Hull", search.method = "ESA")
## plot Hull curve for item 20
plot(Hull.obj, 20)
############### PAA ###############
Hull.obj <- validation(data.obj$dat, MQ, CDM.obj,
method = "Hull", search.method = "PAA")
## plot Hull curve for item 20
plot(Hull.obj, 20)
# }
Run the code above in your browser using DataLab