compute theta EAP estimation

compute theta EAP estimation

Usage

compute_eap_dichotomous(
  bank,
  x,
  D = 1,
  priorPar = c(0, 1),
  lower = -4,
  upper = 4,
  nqp = 33
)

Arguments

bank

matrix of item parameters

x

vector of item responses

D

metric constant can be 1 or 1.702

priorPar

mean and sd for normal distribution default is mean 0 and sd 1

lower

lower bound for numerical integration

upper

upper bound for numerical integration

nqp

number of quadrature points

Note

this function should return the same result as the catR::eapEst function

Examples

bank<-matrix(c(1.75,1.80,1.85,1.90,1.95,
              -3,-2,0,2,3,
               0.1,0.1,0.1,0.1,0.1,
               0.9,0.9,0.9,0.9,0.9),
             nrow=5,
             dimnames=list(c(1,2,3,4,5),
                           c("a","b","c","d")))
response<-c(1,0,0,0,0)
catR::eapEst(bank,response)
#> [1] -0.8954934
compute_eap_dichotomous(bank,response)
#> [1] -0.8954934