p.2pl <- est(Scored, model = "2PL", engine = "ltm")
# plot IRF for all items in red, label with item number
plot(irf(p.2pl), co = "red", label = TRUE)
# plot IRF for items 2, 3, and 7 in different colours
plot(irf(p.2pl[c(2, 3, 7), ]), co = NA)
Run the code above in your browser using DataLab