2012-11-26 66 views
-1

做了方差分析和建築效果表之後,我只需要捕獲那些最小值的表中的5個表格。由於一長串的因素,因爲從x1到x100有大約100個因子,所以我無法看到所有表格。在R中搜索具有最小值的aov效果表

model<-aov(y~., data=data) 
effects<-model.tables(model, "effects") 

任何條款的表的標籤名稱是

names(effects$tables$x1) 
"1"  "2" 

我只需要那些長期的表2捕獲與標籤「2」的最小值。

問題編輯:

effects$tables[1:4] 
$x1 
x1 
     1   2 
-0.01099232 0.01053045 

$x2 
x2 
     1   2 
-0.03292931 0.03321318 

$x3 
x3 
    1   2 
0.2881996 -0.3008399 

$x4 
x4 
     1   2 
-0.010151743 0.009236422 

結果可以是-0.3008399和0.009236422,如果我們只選擇2中給出的表上述最小值出4。

+1

這個問題很清楚,whatdoe你用'最小均標籤「2」的值是什麼值? – mnel

+0

@mnel。我編輯了這個問題。可以看出,每個級別都有效果。我想捕捉2級的效果,但只有那些前5個最小值的術語表或者可以說我們是否將所有這些術語表分類,然後是前5個或最重要的那些。 – Shahzad

回答

0

可能:(在沒有重複的例子,未經測試)

eff2 <- effects$tables$x1[['2']] 
eff2[ order(eff2) ][1:5] 

或許(在conitnued沒有例子的還是未經)

efftabls <- sapply(eff2, '[', 2) 
head(efftabls[order(efftabls)], 5) # 5 used to be the number requested 
+0

但是,此示例僅使用x1項。我有大約100個條款,即從x1到x100或影響$表格[1:100]。我需要用最小值來確定這些術語;可能是5個術語表中的5個最小值。 – Shahzad