2013-03-27 96 views
2

我想根據他們的cronbach-alpha-if-excluded值得到項目的名稱。 ltm包給了我需要的數據,但我無法理解如何以我需要的方式訪問其結果集。 (我知道,其他包有信度的功能,但我只能找到這樣一個處理我隨機丟失的數據。)訪問R中的結構

library(ltm) 
descript(Abortion)['alpha'] 

結果

$alpha 
        value 
All Items  0.8706719 
Excluding Item 1 0.8572997 
Excluding Item 2 0.8223116 
Excluding Item 3 0.8147668 
Excluding Item 4 0.8430359 

我想他們根據訂購返回值並將兩個最大的項目名稱放在列表中。在這種情況下,我需要一種方法來返回badItems < -c('Item 1','Item 4')。

我該怎麼做?

+0

這是什麼意思?名稱(sort(描述(墮胎)['alpha'] $ alpha [-1,1],遞減= TRUE)[c(1,2)]) – user1981275 2013-03-27 12:12:28

+0

@ user1981275你爲什麼不把它作爲解決方案? – 2013-03-27 12:14:49

回答

1

這可能會實現:

names(sort(descript(Abortion)['alpha']$alpha[-1,1], decreasing=TRUE)[c(1,2)]) 
+0

就是這樣。 只需要提取項目名稱。 pat < - 「不包括(。*)$」; bad_item_list = names(sort(descript(Abortion)['alpha'] $ alpha [-1,1],decrease = TRUE)[c(1,2)]) sub(pat,「\\ 1」,bad_item_list [grepl(pat,bad_item_list)]) – SunWuKung 2013-03-27 12:29:25