estimate_mixed function

Estimation of the Mixed Format Model

Estimate the mixed format model

model_mixed_eap(u, items, D = 1.702, priors = c(0, 1),
  bounds_t = c(-4, 4))

model_mixed_map(u, items, D = 1.702, priors = c(0, 1),
  bounds_t = c(-4, 4), iter = 30, conv = 0.001)

Arguments

• u: the response data, 2d marix
• items: a list of 3pl, gpcm, grm items
• D: the scaling constant
• priors: the prior distribution
• bounds_t: the lower- and upper-bound of the parameter
• iter: the maximum number of newton-raphson iterations
• conv: the convergence criterion

Returns

model_mixed_eap returns a list of point estimates and standard error of the ability parameters

model_mixed_map returns a list of point estimates of the ability parameters

Examples

x <- model_mixed_gendata(200, 30, 5, 5, 3)
y <- model_mixed_eap(x$u, x$items)
c('corr'=cor(x$t, y$t), 'rmse'=rmse(x$t, y$t))
x <- model_mixed_gendata(200, 30, 5, 5, 3)
y <- model_mixed_map(x$u, x$items)
c('corr'=cor(x$t, y$t), 'rmse'=rmse(x$t, y$t))