EstCRMperson | R Documentation |
Estimate person parameters for the Continuous Response Model(Samejima,1973)
EstCRMperson(data,ipar,min.item,max.item)
data |
a data frame with N rows and m columns, with N denoting the number of subjects and m denoting the number of items. |
ipar |
a matrix with m rows and three columns, with m denoting the number of items. The first column is the a parameters, the second column is the b parameters, and the third column is the alpha parameters |
min.item |
a vector of length m indicating the minimum possible score for each item. |
max.item |
a vector of length m indicating the maximum possible score for each item. |
Samejima(1973) derived the closed formula for the maximum likelihood estimate (MLE) of the person parameters (ability levels) under Continuous Response Model. Given that the item parameters are known, the parameter for person i (ability level) can be estimated from the following equation:
Please see the manual for the equation!
As suggested in Ferrando (2002), the ability level estimates are re-scaled once they are obtained from the equation above. So, the mean and the variance of the sample estimates equal those of the latent distribution on which the estimation is based.
Samejima (1973) also derived that the item information is equal to:
Please see the manual for the equation!
So, the item information is constant along the theta continuum in the Continuous Response Model.
The asymptotic standard error of the MLE ability level estimate is the square root of the reciprocal of the total test information:
Please see the manual for the equation!
In the Continuous Response Model, the standard error of the ability estimate is same for all examinees unless an examinee has missing data. If the examinee has missing data, the ability level estimate and its standard error is obtained from the available item scores.
EstCRMperson()
returns an object of class "CRMtheta
". An object of class "CRMtheta
" contains the following component:
thetas |
a three column matrix, the first column is the examinee ID, the second column is the re-scaled theta estimate, the third column is the standard error of the theta estimate. |
See the examples of simCRM
for an R code that runs a simulation to examine the recovery of person parameter estimates by EstCRMperson
.
Cengiz Zopluoglu
Ferrando, P.J.(2002). Theoretical and Empirical Comparison between Two Models for Continuous Item Responses. Multivariate Behavioral Research, 37(4), 521-542.
Samejima, F.(1973). Homogeneous Case of the Continuous Response Model. Psychometrika, 38(2), 203-219.
Wang, T. & Zeng, L.(1998). Item Parameter Estimation for a Continuous Response Model Using an EM Algorithm. Applied Psychological Measurement, 22(4), 333-343.
EstCRMitem
for estimating item parameters,
fitCRM
for computing item-fit residual statistics and drawing empirical 3D item category response curves,
plotCRM
for drawing theoretical 3D item category response curves,
simCRM
for generating data under CRM.
## Not run: #Load the dataset EPIA data(EPIA) ##Define the vectors "max.item" and "min.item". The maximum possible ##score was 112 and the minimum possible score was 0 for all items max.item <- c(112,112,112,112,112) min.item <- c(0,0,0,0,0) ##Estimate item parameeters CRM <- EstCRMitem(EPIA, max.item, min.item, max.EMCycle = 500, converge = 0.01) par <- CRM$param ##Estimate the person parameters CRMthetas <- EstCRMperson(EPIA,par,min.item,max.item) theta.par <- CRMthetas$thetas theta.par #Load the dataset SelfEff data(SelfEff) ##Define the vectors "max.item" and "min.item". The maximum possible ##score was 11 and the minimum possible score was 0 for all items max.item <- c(11,11,11,11,11,11,11,11,11,11) min.item <- c(0,0,0,0,0,0,0,0,0,0) ##Estimate the item parameters CRM2 <- EstCRMitem(SelfEff, max.item, min.item, max.EMCycle=200, converge=.01) par2 <- CRM2$param ##Estimate the person parameters CRMthetas2 <- EstCRMperson(SelfEff,par2,min.item,max.item) theta.par2 <- CRMthetas2$thetas theta.par2 ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.