計算eRm（Rating Scale Model）中每個評級類別的服裝均方差

我試圖在R中複製用WINSTEPS執行的已發佈Rasch分析。我正在使用eRm軟件包中的RSM功能。計算eRm（Rating Scale Model）中每個評級類別的服裝均方差

發表的論文的作者使用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函數生成了項目圖，我可以在這個圖中看到每個類別的infit t統計量。但我找不到一種方法來計算每個評級類別的服裝均方。

是否有計算這些函數的函數？我怎樣才能得到他們作爲數字輸出？

非常感謝，亞歷

來源

2016-08-12 Alexandra Dima

讓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]

而且這將修剪你的表。

來源

2016-10-05 20:28:30

感謝Guillaume，我實際上是在爲評級比例模型（不適合）的每個評級類別尋找服裝均方。同時，我檢查了軟件包的作者，這些在eRm中尚不可用。我購買了Winsteps並在同一軟件中複製了分析。如果任何人在eRm中寫出/知道某個函數（該函數獲得了圖表或表格中每個評級類別的服裝均方塊），將會非常適合它（您也可以聯繫作者以將其包含在eRm中）以及該功能將被確認）。 –

計算eRm（Rating Scale Model）中每個評級類別的服裝均方差

回答

相關問題