我试图在R中复制用WINSTEPS执行的已发布的Rasch分析。我在eRm包中使用RSM功能。
已发表论文的作者使用5点李克特响应格式检查项目集的项目评级结构。 他们考虑的一个标准是每个评级类别的装备均方差<2.0（来自http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.424.2811&rep=rep1&type=pdf）。
我设法按照这里的示例计算项目装备和infit：http://media.metrik.de/uploads/incoming/pub/Literatur/R/eRmvig.pdf。我通过plotPWmap函数生成了Item Map，我可以看到该图中每个类别的infit t统计量。但我找不到计算每个评级类别的装备均方的方法。
是否有计算这些功能的功能？如何将它们作为数字输出？
非常感谢， 亚历
答案 0 :(得分：0)
让data_scored为项目和响应的二进制矩阵。假设您要删除适合人群的服装。
library(data.table)
library(eRm)
res <- RM(data_scored)
pp <- person.parameter(res)
pfit.pp <- personfit(pp2)
table.pp <- data.table(pearson.fit=pfit.pp$p.fit,
outfit=pfit.pp$p.outfitZ)
table.pp $ outfit的值将是装备均值平方。所以现在你可以这样做：
table.pp.clean <- table.pp[outfit < 2 & outfit > 2]
这会削减你的桌子。