# NOT RUN {
# # Not run
# #
# # data(dichotomousItemParameters)
# # threePlParameters <- dichotomousItemParameters
# # isNot3Pl <- ((dichotomousItemParameters[['focal']][, 3] == 0) |
# # (dichotomousItemParameters[['reference']][, 3] == 0))
# #
# # threePlParameters[['focal']] <- threePlParameters[['focal']][!isNot3Pl, ]
# # threePlParameters[['reference']] <- threePlParameters[['reference']][!isNot3Pl, ]
# # threePlParameters[['focal']][, 3] <- threePlParameters[['focal']][, 3] + 0.1
# # threePlParameters[['reference']][, 3] <- threePlParameters[['reference']][, 3] + 0.1
# # threePlParameters[['focal']][, 2] <- threePlParameters[['focal']][, 2] + 1.5
# # threePlParameters[['reference']][, 2] <- threePlParameters[['reference']][, 2] + 1.5
# # threePlParameters[['focal']] <- threePlParameters[['focal']][-c(12, 16, 28), ]
# # threePlParameters[['reference']] <- threePlParameters[['reference']][-c(12, 16, 28), ]
# #
# # threePlAse <- list()
# # threePlAse[["focal"]] <- AseIrt(itemParameters = threePlParameters[["focal"]],
# # logistic = TRUE,
# # sampleSize = 10000,
# # irtModel = "3pl")
# # threePlAse[["reference"]] <- AseIrt(itemParameters = threePlParameters[["reference"]],
# # logistic = TRUE,
# # sampleSize = 15000,
# # irtModel = "3pl")
# #
# # set.seed(41568)
# # threePlIpr <- Ipr(itemParameters = threePlParameters, itemCovariances = threePlAse,
# # nReplicates = 100)
# # threePlIpr <- Bound3PlIpr(threePlIpr)
# }
Run the code above in your browser using DataLab