plot {ltm} | R Documentation |
Produces the Item Characteristic or Item Information Curves for fitted IRT models.
## S3 method for class 'grm': plot(x, type = c("ICC", "IIC"), items = NULL, category = NULL, legend = FALSE, cx = "top", cy = NULL, ncol = NULL, col = palette(), lty = 1, ...) ## S3 method for class 'ltm': plot(x, type = c("ICC", "IIC"), items = NULL, legend = FALSE, cx = "topleft", cy = NULL, ncol = 1, col = palette(), lty = 1, ...) ## S3 method for class 'rasch': plot(x, type = c("ICC", "IIC"), items = NULL, legend = FALSE, cx = "topleft", cy = NULL, ncol = 1, col = palette(), lty = 1, ...)
x |
an object inheriting from either class grm , class ltm or class rasch . |
type |
the type of plot; "ICC" refers to Item Response Category Characteristic Curves whereas "IIC" to Item Information Curves. |
items |
a numeric vector denoting which items to plot; if NULL all items are plotted; if 0 and
type = "IIC" the Test Information Curve is plotted. |
category |
a scalar indicating the response category for which the curves should be plotted; if NULL
all categories are considered. |
legend |
logical; if TRUE a legend is printed. |
cx |
the x-coordinate to be passed to legend() . |
cy |
the y-coordinate to be passed to legend() . |
ncol |
the ncol argument of legend() . |
col |
the colors to use; recycling is used if necessary. |
lty |
the types of lines to use; recycling is used if necessary. |
... |
extra graphical parameters to be passed to plot() , lines() , legend() and
text() . |
Item response category characteristic curves show how the probability of responding in the kth category, in each item, changes with the values of the latent variable (ability).
The item information curves indicate the relative ability of an item to discriminate among contiguous trait scores at various locations along the trait continuum. The test information curve, which is the sum of item information curves, provides a visual depiction of where along the trait continuum a test is most discriminating (Reise and Waller, 2002).
If x
is an ltm
object representing a linear two-factor model, then the plot of the standardized
loadings (see coef.ltm
) is also produced.
Since the ...
argument is passed to both plot()
and legend()
an error will occur
if you pass an argument that is not shared by the two functions (e.g., plot(object, legend = TRUE,
cex.lab = 1.1)
will fail). In such cases you could specify the required option(s) using par()
; see
the Examples for plot.rasch()
.
Dimitris Rizopoulos dimitris.rizopoulos@med.kuleuven.be
Reise, S. and Waller, N. (2002) Item response theory for dichotomous assessment data. In Drasgow, F. and Schmitt, N., editors, Measuring and Analyzing Behavior in Organizations. San Francisco: Jossey-Bass.
information
,
grm
,
ltm
,
rasch
# Examples for plot.grm() m <- grm(Science[c(1,3,4,7)]) ## Item Response Category Characteristic Curves for ## the Science data op <- par(mfrow = c(2, 2)) plot(m, lwd = 2, legend = TRUE, ncol = 2) # re-set par() par(op) ## Item Characteristic Curves for the 2nd category, ## and items 1 and 3 plot(m, category = 2, items = c(1, 3), lwd = 2, legend = TRUE, cx = "right") ## Item Information Curves for the Science data; plot(m, type = "IIC", legend = TRUE, cx = "topright", lwd = 2, cex = 1.4) ## Test Information Function for the Science data; plot(m, type = "IIC", items = 0, lwd = 2) ################################################### # Examples for plot.ltm() ## Item Characteristic Curves for the two-parameter logistic ## model; plot only items 1, 2, 4 and 6: m <- ltm(Wirs ~ z1) plot(m, items = c(1, 2, 4, 6), lwd = 3, cex = 1.4) ## Test Information Function under the two-parameter logistic ## model for the Lsat data m <- ltm(Lsat ~ z1) plot(m, type = "IIC", items = 0, lwd = 2, cex.lab = 1.2, cex.main = 1.3) info <- information(m, c(-3, 0)) text(x = 2, y = 0.5, labels = paste("Total Information:", round(info$Total, 3), "\n\nInformation in (-3, 0):", round(info$Info, 3), paste("(", round(100 * info$Proportion, 2), "%)", sep = "")), cex = 1.2) ## Item Characteristic Surfaces for the interaction model: m <- ltm(Wirs ~ z1 * z2) plot(m) ################################################### # Examples for plot.rasch() ## Item Characteristic Curves for the Wirs data; ## plot only items 1, 3 and 5: m <- rasch(Wirs) plot(m, items = c(1, 3, 5), lwd = 3, cex = 1.4) m <- rasch(Lsat) ## Item Characteristic Curves for the Lsat data; ## plot all items plus a legend; increase 'cex.lab' ## and 'cex.main' using par() op <- par(cex.lab = 1.6, cex.main = 2) plot(m, legend = TRUE, cx = "right", lwd = 3, cex = 1.4) par(op) ## Item Information Curves, for the first 3 items; include a legend plot(m, type = "IIC", items = 1:3, legend = TRUE, lwd = 2, cx = "topright") ## Test Information Function plot(m, type = "IIC", items = 0, lwd = 2, cex.lab = 1.1) info <- information(m, c(-2, 0)) text(x = 2, y = 0.5, labels = paste("Total Information:", round(info$Total, 3), "\n\nInformation in (-2, 0):", round(info$Info, 3), paste("(", round(100 * info$Prop, 2), "%)", sep = "")), cex = 1.2)