p.2pl <- est(Scored, model = "2PL", engine = "ltm")
# plot items 1:5 in different colours, label
npp(Scored, items = 1:5, co = NA, label = TRUE)
# For item 7, compare npp with the 2PL parametric IRF
npp(Scored, items = 7, bands = TRUE)
plot(irf(ip = p.2pl[7, ]), co = 3, add = TRUE)
Run the code above in your browser using DataLab