检查矩阵元素的长度

时间:2017-07-17 13:47:08

标签: r matrix

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矩阵对象。每个单元格都包含一个数值列表。我想在i, j对象中找到索引mat,指示哪些单元格包含的向量是&gt; 1但是&lt; 5个元素长。换句话说,我希望看到指数:(1,2),(1,3)和(3,1)。

我已尝试使用which(length(mat) < 5 & length(mat) > 1),但这不起作用。

2 个答案:

答案 0 :(得分:5)

一种可能的解决方案是使用基础R中的lengths。因为在这种情况下有什么输出是什么,这里有我的详细信息,

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

如果您使用的是版本&lt; 3.4.1,然后采用两步法(如@lmo提出的方法),即

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

答案 1 :(得分:2)

这是我的解决方案:

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