##from rdm, item parameters denote 4 item 1PL test, cut score at theta=.5
##only return marginal indices
params<-matrix(c(1,1,1,1,-2,1,0,1,0,0,0,0),4,3)
rdm<-sim(params, rnorm(100))
class.Rud(2, params, rdm = rdm)$Marginal
##or from 40 quadrature points and weights
quad <- normal.qu(40)
class.Rud(2, params, quadrature = quad)$Marginal
Run the code above in your browser using DataLab