2016-09-17 115 views
1
m <- lm(Sepal.Length ~ Sepal.Width + Petal.Width, data = iris) 
beta = dfbetas(m) 
> head(beta) 
    (Intercept) Sepal.Width Petal.Width 
1 -0.0018633253 0.0054762565 -0.0096031648 
2 0.0094916858 -0.0062007468 -0.0137816086 
3 -0.0221770886 0.0069280848 0.0540485812 
4 -0.0408776612 0.0219247324 0.0731671391 
5 0.0071436202 -0.0134636336 0.0150509697 
6 0.0006264958 -0.0007979264 0.0001755277 


apply(data.frame(beta), 2, function(x) which(abs(x) < 0.1632993)) 

我有一个名为beta的矩阵,它包含3列。在将其转换为data.frame类型后,我想使用which观测值的函数查找绝对值为< 0.1632993的索引。基本上我想要每列的索引列表。但我的应用功能似乎没有做我想做的事情。R:找到data.frame中列的索引

+0

它可以是任何形式,只要它告诉我,这观察(S)满足我的腹肌的'标准(价值)<我的data.frame的每列的<0.1632993'。 – Adrian

+0

带'arr.ind = TRUE'的'which'给出行/列索引 – akrun

回答

1

我们可以对整个数据集使用whicharr.ind=TRUE得到row/column指标

which(abs(beta) < 0.1632993, arr.ind=TRUE)