3

考虑这个简单的稀疏矩阵

> (X <- sparseMatrix(c(1, 2, 1), c(1, 1, 2), x = 0:2))
2 x 2 sparse Matrix of class "dgCMatrix"

[1,] 0 2
[2,] 1 .

如何将其转换为指示相应元素是否为非空的矩阵?这是我现在正在做的事情,但为 0 并不等于“空”,这种方法并没有区分它们。

> (Y <- X != 0)
2 x 2 sparse Matrix of class "lgCMatrix"

[1,] : |
[2,] | .

为了澄清,所需的输出可能只包含TRUEorFALSE但不包含NA. 它可以是 amatrix或 a sparseMatrix。甚至更优选地,它可以是list,在这种情况下,每个槽对应于 的列X。例如,答案X应该是

     [,1]  [,2]
[1,] TRUE  TRUE
[2,] TRUE FALSE

或者

$`1`
[1] TRUE TRUE

$`2`
[1]  TRUE FALSE
4

2 回答 2

3
Y <- as(X, "lgCMatrix") #should be more efficient than X != 0
Y@x[] <- TRUE #set all values to TRUE
as.matrix(Y)
#     [,1]  [,2]
#[1,] TRUE  TRUE
#[2,] TRUE FALSE
于 2019-08-02T06:33:44.847 回答
2
pex <- function(p) {
    dp <- diff(p)
    rep(seq_along(dp), dp)
}

m = matrix(FALSE, nrow = nrow(X), ncol = ncol(X))
m[cbind(X@i + 1, pex(X@p))] = TRUE
m
#     [,1]  [,2]
#[1,] TRUE  TRUE
#[2,] TRUE FALSE

pex是从这里

于 2019-08-02T03:31:12.580 回答