One implementation: each row goes to morie_psy_map_theta().
This is the MODE where morie_psy_theta_eap() is the MEAN, and
they differ whenever the posterior is skewed – short tests and
extreme patterns, which is most of the time.
Usage
morie_psy_theta_map(X, items, prior = c(0, 1))
Arguments
- X
response matrix.
- items
item-parameter matrix.
- prior
c(mean, sd).
Value
list: theta, se, n_examinees, n_items, mode_not_mean,
alias_of, method.
References
Mislevy (1986), Psychometrika 51:177-195.
Examples
X <- matrix(c(1, 0, 1, 1), 1)
morie_psy_theta_map(X, cbind(rep(1, 4), c(-1, -0.5, 0.5, 1)))$theta
#> [1] 0.5403839
