我有一个大的稀疏矩阵。在用一些数学运算填充矩阵后,我意识到由于除以零错误,我有一些无穷大的值。如何检查此矩阵的非有限值?
这是一个玩具矩阵。
A <- Matrix(nrow = 150000, ncol = 150000, data = 0, sparse = TRUE)
A[1, 1] = Inf
A[1, 3] = NA
A[2, 1] = -Inf
试图找到它的非有限值给我一个错误:
test <- A[!is.finite(A)]
#Error: cannot allocate vector of size 83.8 Gb
我也试着逐行扫描这个矩阵,但它需要很长时间。
library(magrittr)
for(i in 1:nrow(A)){
if((
A[i, ] %>% .[!is.finite(.)] %>% length
) > 0) print(i)
}
然后我尝试并行运行它,但我认为这太过分了。更何况,还需要很长时间。
library(parallel)
library(magrittr)
numCores <- detectCores() - 1
cl <- makeCluster(numCores)
clusterExport(cl, c("A"))
clusterEvalQ(cl, library(magrittr))
out <- A %>% nrow %>% seq %>% parLapply(cl, X = ., function(i) A[i, ] %>% .[!is.finite(.)])
我该如何继续?
最佳答案
如果我们想知道一个稀疏矩阵 A
有任何 Inf
, -Inf
, NaN
或 NA
, 我们可以做到
any(!is.finite(A@x))
#[1] TRUE
如果我们也想知道他们的位置,我们可以这样做
subset(summary(A), !is.finite(x))
i j x
1 1 1 Inf
2 2 1 -Inf
3 1 3 NA
备注:
参见 R: element-wise matrix division is.infinite
之间的区别, !is.finite
, is.na
和 is.nan
.
关于r - 检查非有限值的大型稀疏矩阵的有效方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/73253735/