##with simple 4 item 1PL simulated test:
## find True classification accuracy of theta.hat and total score
## use a cut score of theta = 0 (assume true score for total score cut score)
params<-matrix(c(1,1,1,1,-2,1,0,1,0,0,0,0),4,3)
theta <-rnorm(100)
rdm<-sim(params, theta)
theta.hat <- MLE(rdm,params)
total.score <- rowSums(rdm)
t.cut<-0
x.cut <- sum(irf(params, t.cut)$f)
TOtable.F(theta, theta.hat, t.cut, t.cut)$P
TOtable.F(theta, total.score, t.cut, x.cut)$P
# for 1PL model total score is sufficient for theta.hat, repeat with 2PL or 3PL items...
Run the code above in your browser using DataLab