这个问题与这里提出的问题非常相似:Element-wise mean of several big.matrix objects in R。
但是,现在我的目标是计算 18 个文件支持的 big.matrix 对象(每个 > 4GB)的元素平均值,同时考虑到 na.rm = TRUE。
对于正常大小的矩阵,我会执行以下操作:
#Generate a list of 3 example matrices
my.list <- list(matrix(c(1:9), nrow = 3),
matrix(c(10:18), nrow = 3),
matrix(c(19:27), nrow = 3))
my.list[[1]][1,1] <- NA
my.list[[3]][2,1] <- NA
#Take element-wise mean with na.rm = TRUE
test <- apply(simplify2array(my.list), 1:2, mean, na.rm= TRUE)
#This results in:
> test
[,1] [,2] [,3]
[1,] 14.5 13 16
[2,] 6.5 14 17
[3,] 12.0 15 18
我正在使用具有 16GB RAM 的 64x Windows10 机器。
任何建议都非常受欢迎。谢谢!