myAbilities <- rnorm(1000, 0, 1)
myItems <- runif (10, -2, 3)
myExProbs <- matrix(rnorm(100, 0, 1), 10, 10)
myExPerturb <- list(data.frame(item = 1, type = "flat", pivot = 0.5, factor = 25))
myProbs <- perturbProbs(myProbs = myExProbs, perturbP = myExPerturb)
myPerturbDelta <- perturbDelta(myAbilities = myAbilities, myItems = myItems, myProbs = myProbs)
Run the code above in your browser using DataLab