有效地将稀疏矩阵的元素设置为0

洪井井

假设我有一个稀疏矩阵：

library(Matrix)
m <- Matrix(rpois(100, 1), 10, 10, sparse=TRUE)
m
#10 x 10 sparse Matrix of class "dgCMatrix"

# [1,] . . 2 3 1 1 1 . 3 1
# [2,] 1 . 1 . 1 . 2 . . .
# [3,] 1 4 1 . . . . 5 1 1
# [4,] . 1 1 . . . 4 1 1 1
# [5,] 1 . 4 . 1 . . 1 . 2
# [6,] 1 1 2 1 . 1 . 1 2 1
# [7,] 4 2 1 1 . . . 1 1 3
# [8,] 2 2 1 1 1 1 . 1 1 2
# [9,] 2 1 1 . 3 . 2 . . 1
#[10,] . 2 1 1 . . . 1 3 3

我想将所有小于或等于1的元素设置为零。这是一种安全的方法吗？请注意，m2从技术上讲，它不再是压缩的稀疏矩阵，因为某些“非零”元素现在实际上为零。

m2@x[m2@x < 2] <- 0
m2
#10 x 10 sparse Matrix of class "dgCMatrix"

# [1,] . . 2 3 0 0 0 . 3 0
# [2,] 0 . 0 . 0 . 2 . . .
# [3,] 0 4 0 . . . . 5 0 0
# [4,] . 0 0 . . . 4 0 0 0
# [5,] 0 . 4 . 0 . . 0 . 2
# [6,] 0 0 2 0 . 0 . 0 2 0
# [7,] 4 2 0 0 . . . 0 0 3
# [8,] 2 2 0 0 0 0 . 0 0 2
# [9,] 2 0 0 . 3 . 2 . . 0
#[10,] . 2 0 0 . . . 0 3 3

具体地讲，我不希望创建大规模指示器矩阵，使该元件以改变的位置（如在m[m < 2] <- 0手动通过元素），或循环。

罗兰

我认为您的方法没有任何不安全之处，甚至可能会建议您这样做。如果您想要一个完全稀疏的矩阵，可以执行以下操作：

drop0(m2)
#10 x 10 sparse Matrix of class "dgCMatrix"
#                         
# [1,] 2 . . . . . . . . .
# [2,] . . . 2 . . . . . .
# [3,] . . 2 . . . . . 2 .
# [4,] 2 2 . . 4 . . 2 2 3
# [5,] . 4 . 3 . 3 . . . .
# [6,] 2 . . . . 3 2 . . .
# [7,] . . . . . . 2 . . 2
# [8,] . . . . . . 2 . . 4
# [9,] 2 . . 2 . . . . . 2
#[10,] . . . . . . 2 . . .

本文收集自互联网，转载请注明来源。

如有侵权，请联系 [email protected] 删除。