Description Usage Arguments Details Value Note Author(s) References See Also Examples

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.

1 2 |

`th` |
numeric: the ability value. |

`it` |
numeric: a suitable matrix of item parameters. See |

`model` |
either |

`D` |
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 `NA`

values.

Under dichotomous IRT models, a list with four arguments:

`Pi` |
the vector with response probabilities (one value per item). |

`dPi` |
the vector with first derivatives of the response probabilities (one value per item). |

`d2Pi` |
the vector with second derivatives of the response probabilities (one value per item). |

`d3Pi` |
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 `Item1`

, `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.

David Magis

[email protected]

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")
``` |

catR documentation built on July 28, 2018, 5:01 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.