-1

在进行方差分析和构建效果表之后,我只需要用最小值捕获这些术语表中的 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。

4

1 回答 1

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
于 2012-11-26T02:42:33.637 回答