set.seed(123)
library(Qval)
## generate Q-matrix and data
K <- 4
I <- 20
IQ <- list(
P0 = runif(I, 0.2, 0.4),
P1 = runif(I, 0.6, 0.8)
)
# \donttest{
Q <- sim.Q(K, I)
data <- sim.data(Q = Q, N = 500, IQ = IQ, model = "GDINA", distribute = "horder")
MQ <- sim.MQ(Q, 0.1)
CDM.obj <- CDM(data$dat, MQ)
############### ESA ###############
Hull.obj <- validation(data$dat, MQ, CDM.obj, method = "Hull", search.method = "ESA")
## plot Hull curve for item 20
plot(Hull.obj, 20)
############### PAA ###############
Hull.obj <- validation(data$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