在R中搜索具有最小值的aov效果表

时间:2012-11-26 02:08:45

标签: r anova

在进行方差分析和构建效果表之后,我需要使用最小值捕获这些术语表中的5个。由于一系列因素;因为从x1到x100有大约100个因素,所以我无法看到所有表格。

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 

如果我们只从上面给出的4个表中选择2个最小值,结果可能是-0.3008399和0.009236422。

1 个答案:

答案 0 :(得分:0)

也许:(在没有可重现的例子的情况下未经测试)

eff2 <- effects$tables$x1[['2']]
eff2[ order(eff2) ][1:5]
或者或许(在没有例子的情况下仍未经测试)

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