## Not run:
#
# # item parameters
# ipar2 <- ipar1 <- seq(-3, 3, length.out = 20)
# # model differential item function (DIF)
# ipar2[10] <- ipar1[11]
# ipar2[11] <- ipar1[10]
# # simulation for b = 100, 200, 300, 400, 500
# simres <- pwr.rasch(seq(100, 500, by = 100), ipar = list(ipar1, ipar2))
# plot(simres)
# ## End(Not run)
Run the code above in your browser using DataLab