data(Scored)
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