## load data
data("toydata")
## fit 2PL model using mirt
mirt_model <- mirt::mirt(toydata$resp, model = 1)
## test DIF along a contiuous covariate
DIFtest <- mstDIF(mirt_model, DIF_covariate = toydata$group_cont,
method = "analytical")
## print
DIFtest
## summary
summary(DIFtest)
Run the code above in your browser using DataLab