Skip to contents

EAP ability scores for new respondents under a fitted 2PL model

Usage

morie_irt_eap(fit, responses, n_quad = 41L)

Arguments

fit

A morie_irt_2pl fit (or any list with discrimination and difficulty).

responses

0/1 matrix of new response patterns.

n_quad

Quadrature points.

Value

A data frame with theta (EAP) and se.

Examples

set.seed(63)
th <- rnorm(1000)
a <- c(0.8, 1.2, 1.6, 1.0, 1.4); b <- c(-1, -0.5, 0, 0.5, 1)
X <- vapply(seq_along(a), function(j) rbinom(1000, 1, plogis(a[j] * (th - b[j]))), numeric(1000))
fit <- morie_irt_2pl(X)
morie_irt_eap(fit, X)$theta[1:5]
#> [1]  1.3633218 -0.8728535 -0.3110472 -0.8728535  0.2163830