This command returns the item response probabilities for a given ability value and a given matrix of item parameters under either the 4PL model or any suitable polytomous model. Numerical values of the first, second and third derivatives of the response probabilities are also returned.
numeric: the ability value.
numeric: a suitable matrix of item parameters. See Details.
numeric: the metric constant. Default is
Whatever the IRT model, the response probabilities and first, second, and third derivatives are computed algebraically. These derivatives are necessary for both the estimation of ability and the computation of related standard errors.
Dichotomous IRT models are considered whenever
model is set to
NULL (default value). In this case,
it must be a matrix with one row per item and four columns, with the values of the discrimination, the difficulty, the pseudo-guessing and the inattention parameters (in this order). These are the parameters of the four-parameter logistic (4PL) model
(Barton and Lord, 1981).
Polytomous IRT models are specified by their respective acronym:
"GRM" for Graded Response Model (Samejima, 1969),
"MGRM" for Modified Graded Response Model (Muraki, 1990),
"PCM" for Partical Credit Model (Masters, 1982),
"GPCM" for Generalized Partial Credit Model (Muraki, 1992),
"RSM" for Rating Scale Model (Andrich, 1978) and
"NRM" for Nominal Response Model (Bock, 1972). The
it still holds one row per item, end the number of columns and their content depends on the model. See
genPolyMatrix for further information and illustrative examples of suitable polytomous item banks.
The output list contains the response probabilities and the respective derivatives. In case of dichotomous models, only a vector of such values is returned, with one value per item. In case of polytomous models, matrices are returned instead, with one row per item and one column per response category. In case of unequal numbers of response categories (which may happen under GRM, PCM, GPCM and NRM), values for empty response categories are returned as
Under dichotomous IRT models, a list with four arguments:
the vector with response probabilities (one value per item).
the vector with first derivatives of the response probabilities (one value per item).
the vector with second derivatives of the response probabilities (one value per item).
the vector with third derivatives of the response probabilities (one value per item).
Under polytomous IRT models, the aforementioned vectors are replaced by matrices with one row per item (labeled as
Item2 etc.) and one row per response category.
For dichotomous IRT models, response probabilites exactly equal to zero are returned as
1e-10 values, as well as probabilities exactly equal to one which are returned as
1-1e-10 values. This is to permit the computation of ability estimates (with the
thetaEst function) in such extreme cases.
Many thanks to Pan Tong (University of Texas MD Anderson Cancer Center, USA) who noticed this problem.
Andrich, D. (1978). A rating formulation for ordered response categories. Psychometrika, 43, 561-573. doi: 10.1007/BF02293814
Barton, M.A., and Lord, F.M. (1981). An upper asymptote for the three-parameter logistic item-response model. Research Bulletin 81-20. Princeton, NJ: Educational Testing Service.
Bock, R. D. (1972). Estimating item parameters and latent ability when responses are scored in two or more nominal categories. Psychometrika, 37, 29-51. doi: 10.1007/BF02291411
Haley, D.C. (1952). Estimation of the dosage mortality relationship when the dose is subject to error. Technical report no 15. Palo Alto, CA: Applied Mathematics and Statistics Laboratory, Stanford University.
Magis, D. and Barrada, J. R. (2017). Computerized Adaptive Testing with R: Recent Updates of the Package catR. Journal of Statistical Software, Code Snippets, 76(1), 1-18. doi: 10.18637/jss.v076.c01
Magis, D., and Raiche, G. (2012). Random Generation of Response Patterns under Computerized Adaptive Testing with the R Package catR. Journal of Statistical Software, 48 (8), 1-31. doi: 10.18637/jss.v048.i08
Masters, G. N. (1982). A Rasch model for partial credit scoring. Psychometrika, 47, 149-174. doi: 10.1007/BF02296272
Muraki, E. (1990). Fitting a polytomous item response model to Likert-type data. Applied Psychological Measurement, 14, 59-71. doi: 10.1177/014662169001400106
Muraki, E. (1992). A generalized partial credit model: Application of an EM algorithm. Applied Psychological Measurement, 16, 19-176. doi: 10.1177/014662169201600206
Samejima, F. (1969). Estimation of latent ability using a response pattern of graded scores. Psychometrika Monograph (vol. 17).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
## Dichotomous models ## # Loading the 'tcals' parameters data(tcals) # Selecting item parameters only tcals <- as.matrix(tcals[,1:4]) # Response probabilities and derivatives (various th and D values) Pi(th = 0, tcals) Pi(th = 0, tcals, D = 1.702) Pi(th = 1, tcals) ## Polytomous models ## # Generation of an item bank under GRM with 100 items and at most 4 categories m.GRM <- genPolyMatrix(100, 4, "GRM") m.GRM <- as.matrix(m.GRM) # Computation of probabilities and derivatives for ability level 0 Pi(0, m.GRM, model = "GRM") # Loading the cat_pav data data(cat_pav) cat_pav <- as.matrix(cat_pav) # Computation of probabilities and derivatives for ability level 1 Pi(1, cat_pav, model = "GPCM")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.