# NOT RUN {
## Load Verbal Aggression data
data("VerbalAggression", package = "psychotools")
## Fit Rasch, rating scale and partial credit
## model to VerbalAggression data
rmmod <- raschmodel(VerbalAggression$resp2)
rsmod <- rsmodel(VerbalAggression$resp)
pcmod <- pcmodel(VerbalAggression$resp)
## Profile plots of the item parameters of
## the three fitted IRT models
plot(rmmod, type = "profile", what = "items", col = 4)
plot(rsmod, type = "profile", what = "items", col = 2, add = TRUE)
plot(pcmod, type = "profile", what = "items", col = 3, add = TRUE)
legend(x = "topleft", legend = c("RM", "RSM", "PCM"),
col = 1, bg = c(4, 2, 3), pch = 21, bty = "n")
## Profile plots of the threshold parameters of type 'mode'
plot(rmmod, type = "profile", what = "thresholds",
parg = list(type = "mode"))
plot(rsmod, type = "profile", what = "thresholds",
parg = list(type = "mode"))
plot(pcmod, type = "profile", what = "thresholds",
parg = list(type = "mode"))
## Profile plot of the discrimination parameters
## of a dichotomous Rasch model (all zero by definition)
plot(rmmod, type = "profile", what = "discrimination")
# }
Run the code above in your browser using DataLab