2017-07-17 119 views
2
mat <- structure(list(c(1, 2, 3, 4, 5), 2, c(3, 2, 1), numeric(0), numeric(0), 
    numeric(0), c(1, 2, 3, 6), c(1, 2, 3, 4, 5), 1, numeric(0), 
    numeric(0), numeric(0), c(3, 4, 2), 3, c(1, 2, 3, 4, 5), 
    numeric(0), numeric(0), numeric(0), numeric(0), numeric(0), 
    numeric(0), 1.358, numeric(0), numeric(0), numeric(0), numeric(0), 
    numeric(0), numeric(0), 0.0223257970827299, numeric(0), numeric(0), 
    numeric(0), numeric(0), numeric(0), numeric(0), 1.493), .Dim = c(6L, 
6L)) 

> mat 
    [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  
[1,] Numeric,5 Numeric,4 Numeric,3 Numeric,0 Numeric,0 Numeric,0 
[2,] 2   Numeric,5 3   Numeric,0 Numeric,0 Numeric,0 
[3,] Numeric,3 1   Numeric,5 Numeric,0 Numeric,0 Numeric,0 
[4,] Numeric,0 Numeric,0 Numeric,0 1.358  Numeric,0 Numeric,0 
[5,] Numeric,0 Numeric,0 Numeric,0 Numeric,0 0.0223258 Numeric,0 
[6,] Numeric,0 Numeric,0 Numeric,0 Numeric,0 Numeric,0 1.493 

我有一个6x6矩阵对象。每个单元格都包含一个数值列表。我想找到mat对象中的指数i, j,指示哪些单元格包含大于1但长度为5个元素的矢量。换句话说,我希望看到指数:(1,2),(1,3)和(3,1)。检查矩阵元素的长度

我试过使用which(length(mat) < 5 & length(mat) > 1)但没有奏效。

+2

FWIW'长度(垫)'给你的长度 –

+0

所以我想只是'这(长度(垫的一个简单的载体)<5的长度(垫)> 1, arr.ind = TRUE)' – Sotos

+2

@Sotos'长度'只会返回一个向量。这两行将起作用:'lMat < - 矩阵(长度(mat),6)'然后'哪个(lMat> 1&lMat <5,arr.ind = TRUE)'。 – lmo

回答

5

一个可能的解决方案是使用lengths从基地R.由于有什么在这种情况下,这里有我的详细信息输出什么各种评论,

sessionInfo() 
R version 3.4.1 (2017-06-30) 
Platform: x86_64-w64-mingw32/x64 (64-bit) 
Running under: Windows >= 8 x64 (build 9200) 

mat 
    [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  
[1,] Numeric,5 Numeric,4 Numeric,3 Numeric,0 Numeric,0 Numeric,0 
[2,] 2   Numeric,5 3   Numeric,0 Numeric,0 Numeric,0 
[3,] Numeric,3 1   Numeric,5 Numeric,0 Numeric,0 Numeric,0 
[4,] Numeric,0 Numeric,0 Numeric,0 1.358  Numeric,0 Numeric,0 
[5,] Numeric,0 Numeric,0 Numeric,0 Numeric,0 0.0223258 Numeric,0 
[6,] Numeric,0 Numeric,0 Numeric,0 Numeric,0 Numeric,0 1.493  

str(mat) 
List of 36 
$ : num [1:5] 1 2 3 4 5 
$ : num 2 
$ : num [1:3] 3 2 1 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num [1:4] 1 2 3 6 
$ : num [1:5] 1 2 3 4 5 
$ : num 1 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num [1:3] 3 4 2 
$ : num 3 
$ : num [1:5] 1 2 3 4 5 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num 1.36 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num 0.0223 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num(0) 
$ : num 1.49 
- attr(*, "dim")= int [1:2] 6 6 

然后,

m2 <- lengths(mat) 
which(m2 < 5 & m2 > 1, arr.ind = TRUE) 
#  row col 
#[1,] 3 1 
#[2,] 1 2 
#[3,] 1 3 

在如果您使用的版本是< 3.4.1,那么两步法(如由@lmo提出的方法)将起作用,即

m1 <- matrix(lengths(mat), 6) 
which(m1 > 1 & m1 < 5, arr.ind=TRUE) 
+1

美丽,谢谢! – Adrian

2

这里是我的解决方案:

lMat <- apply(mat, 2, function(x) sapply(x, length)) 
which(lMat > 1 & lMat < 5, arr.ind = TRUE)