set.seed(123)
library(Qval)
# \donttest{
################################################################
# Example 1: summary a CDM object #
################################################################
Q <- sim.Q(3, 20)
IQ <- list(P0 = runif(20, 0, 0.2), P1 = runif(20, 0.8, 1))
data.obj <- sim.data(Q, N = 500, IQ = IQ,
model = "GDINA", distribute = "horder")
CDM.obj <- CDM(data.obj$dat, Q, model = "GDINA", method = "BM")
summary(CDM.obj)
CDM.updated <- update(CDM.obj, method = "EM", maxitr = 1000)
summary(CDM.updated)
################################################################
# Example 2: summary a validation object #
################################################################
MQ <- sim.MQ(Q, 0.1)
CDM.obj2 <- CDM(data.obj$dat, MQ)
validation.obj <- validation(data.obj$dat, MQ, CDM.obj2,
method = "GDI")
summary(validation.obj)
validation.updated <- update(validation.obj, method = "Hull")
summary(validation.updated)
################################################################
# Example 3: summary a sim.data object #
################################################################
data.obj2 <- sim.data(Q = sim.Q(3, 10), N = 1000)
summary(data.obj2)
data.updated <- update(data.obj2, N = 200)
summary(data.updated)
################################################################
# Example 4: summary a fit object #
################################################################
fit.obj <- fit(data.obj$dat, Q, model = "GDINA")
summary(fit.obj)
data.updated <- update(fit.obj, model = "DINA")
summary(data.updated)
# }
Run the code above in your browser using DataLab