This function calculates the proportion of person variance that is not due to error. The concept of person separation reliability is very similar to reliability indices such as Cronbach's \(\alpha\).
SepRel(pobject)# S3 method for eRm_SepRel
print(x, …)
# S3 method for eRm_SepRel
summary(object, …)
Object of class ppar
(see person.parameter
).
Object of class eRm_SepRel
.
Object of class eRm_SepRel
.
Further arguments.
SepRel
returns a list object of class eRm_SepRel
containing:
the person separation reliability,
the squared standard deviation (i.e., total person variability),
the mean square measurement error (i.e., model error variance).
Returns the person separation reliability \(\frac{\mathrm{SSD}-\mathrm{MSE}}{\mathrm{SSD}}\) where SSD is the squared standard deviation and MSE the mean squared error.
eRm uses a conditional maximum likelihood (CML) framework and handles missing values as separate groups during the estimation of item parameters. Person parameters are computed in a second step using unconditional or joint maximum likelihood (UML or JML) estimation with item parameters assumed to be known from the first step. Other programs might do JML to estimate item and person parameters at the same time, or employ marginal maximum likelihood MML to estimate item parameters, assuming a certain distribution for person parameters. In the latter case person parameters might be obtained by various methods like EAP, MAP, …. Even CML-based programs yield different values, for example, if they use Warm's weighted maximum likelihood estimation WLE to compute person parameters in the second step.
The bottom line is that, since there is not “definite” solution for this problem, you will end up with different values under different circumstances. This is another reason to take results and implications with a grain of salt.
Wright, B.D., and Stone, M.H. (1999). Measurement essentials. Wide Range Inc., Wilmington. (http://www.rasch.org/measess/me-all.pdf 28Mb).
# NOT RUN {
# Compute Separation Reliability for a Rasch Model:
pers <- person.parameter(RM(raschdat1))
res <- SepRel(pers)
res
summary(res)
# }
Run the code above in your browser using DataLab