Sometimes it is useful to simulate tests and items in Item Response Theory, this vignette specifies how to simulate tests and use or interpret the output, using the IRTpp package.
To simulate a test use the function simulateTest as follows
library(IRTpp)
test <- simulateTest(model="2PL",items=10,individuals=1000)This runs a simulation with a 2 parameter logistic model, using 10 items and 1000 individuals. A list is returned with the item parameters, the test (Correct responses of the individuals are marked as a 1, in dichotomic models). For the test we just simulated, item parameters are :
test$itempars## $a
## [1] 0.4969640 3.2754086 0.3422372 2.2515297 1.2174463 3.2933905 2.8389663
## [8] 3.5995106 0.3512622 0.6957358
##
## $b
## [1] 0.96848929 1.96536781 0.88632253 0.03054026 0.94955742
## [6] -0.57673001 0.72152335 0.61909890 0.02100582 0.27485018
##
## $c
## [1] 0 0 0 0 0 0 0 0 0 0
##
## $model
## [1] "2PL"Where test$itempars$a indicates the parameters belonging to the discrimination parameter, and so on. Notice that the c parameter is a vector of 0, since we are simulating a 2 parameter logistic model.
The test element of the list returned by simulateTest contain a list of tests generated, notice that if you simulate only one test, the response data will be in :
responses <-test$test
responses[1:10,]## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
## [1,] 1 0 0 0 0 0 0 0 0 0
## [2,] 1 0 1 1 0 1 0 0 1 0
## [3,] 1 0 0 1 1 1 0 0 1 0
## [4,] 1 0 1 1 0 1 1 1 1 0
## [5,] 0 0 0 1 1 1 1 1 1 1
## [6,] 1 0 0 0 0 0 0 0 1 0
## [7,] 1 0 0 0 1 1 1 0 1 1
## [8,] 1 0 0 1 0 1 1 0 1 1
## [9,] 0 0 1 1 0 1 0 1 0 1
## [10,] 1 0 1 0 0 1 0 0 0 0responses[1:10,] are the first 10 rows of the simulated test.
IRTpp has built-in capabilities for simulating multiple tests at once, use the option reps to simulate multiple tests at once. All the tests are simulated with the same individual and item parameters.
For example simulating 3 tests of the one parameter logistic model :
t = simulateTest(model="1PL",items=4,individuals=5)
length(t$test)## [1] 20Sometimes you want the item parameters to fall in a range, use the boundaries option to achieve this. The boundaries option receives a list of the parameter bounds you want to implement. Let’s look at the parameter boundaries in an ordinary 3PL simulated test.
t3 = simulateTest(model="3PL",items=500,individuals=10);
summary(t3$itempars$c)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000825 0.066092 0.124335 0.125722 0.192200 0.249775As you can see, the item parameters will most likely fall inside 0 and 0.35 which are the package defaults, however we can impose a lower boundary for the c parameter as follows:
bd = list(c_lower=0.2)
t3 = simulateTest(model="3PL",items=500,individuals=10,boundaries=bd);
summary(t3$itempars$c)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.2002 0.2132 0.2249 0.2251 0.2384 0.2500Notice how the lower boundary has been imposed in the c parameter.
When simulating tests, sometimes, perfect responses by an individual can affect procedures made with the tests themselves, or it is required for the simulation that the individuals do not answer less (or more) than a percent of the items, to impose restrictions in the individual classic scores use the threshold parameter.
t3 = simulateTest(model="3PL",items=10,individuals=100,threshold=0.2);This threshold ensures that the individuals do not answer less than 20% of the answers or more than 80% of the answers.
response <- t3$test
summary(rowSums(response))## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.00 3.00 4.00 4.34 6.00 10.00