我在R中有一个很大的稀疏矩阵(dgCMatrix
),我想知道这个矩阵每个元素的指数.更准确地说,我想对矩阵的每个元素执行1-exp(-x).不幸的是,当我在R中这样做时,有一个sparse->dense coercion
,它需要大量的时间和内存(见下面的例子).
library(Matrix)
i <- sample(20000, 20000); j <- sample(20000, 20000); x <- 7 * (1:20000)
A <- sparseMatrix(i, j, x = x)
1 - exp(-A)
有没有办法在R中避免这种胁迫?当1-exp(0)为0时,或许可以保持稀疏性.