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