# loading packages
library(VGAM)
# loading data
data(Science, package = "mirt")
# total score calculation
score <- rowSums(Science)
Science[, 1] <- factor(Science[, 1], levels = sort(unique(Science[, 1])), ordered = TRUE)
# adjacent category logit model for item 1
fit <- vglm(Science[, 1] ~ score, family = acat(reverse = FALSE, parallel = TRUE))
# coefficients for item 1
coef(fit)
plotAdjacent(fit, matching.name = "Total score")
Run the code above in your browser using DataLab