我正在使用bigmemory和biganalytics软件包,并尝试计算
没有使用
big.matrix
对象的平均值。 biganalytics的文档(例如?biganalytics
)建议mean()
可用于big.matrix
对象,但失败了。x <- big.matrix(5, 2, type="integer", init=0,
+ dimnames=list(NULL, c("alpha", "beta")))
x
# An object of class "big.matrix"
# Slot "address":
# <pointer: 0x00000000069a5200>
x[,1] <- 1:5
x[,]
# alpha beta
# [1,] 1 0
# [2,] 2 0
# [3,] 3 0
# [4,] 4 0
# [5,] 5 0
mean(x)
# [1] NA
# Warning message:
# In mean.default(x) : argument is not numeric or logical: returning NA
虽然有些东西可以正常工作:
colmean(x)
# alpha beta
# 3 0
sum(x)
# [1] 15
mean(x[])
# [1] 1.5
mean(colmean(x))
# [1] 1.5
没有使用
mean()
函数,使用mean(colmean(x))
似乎是下一个最佳选择:# try it on something bigger
x = big.matrix(nrow=10000, ncol=10000, type="integer")
x[] <- c(1:(10000*10000))
mean(colmean(x))
# [1] 5e+07
mean(x[])
# [1] 5e+07
system.time(mean(colmean(x)))
# user system elapsed
# 0.19 0.00 0.19
system.time(mean(x[]))
# user system elapsed
# 0.28 0.11 0.39
可以想象,对于具有大量列的矩形矩阵,mean()
可能仍然更快。
不知道为什么 mean()
对我不起作用,有什么想法吗?
sessionInfo()
吗?我刚刚在 R 3.2.2 上安装了bigmemory
和biganalytics
,没有出现任何问题。 - cdetermanbiganalytics
解决了问题。不确定出了什么问题。 - chriss