在 R 中传递大型矩阵的替代方法答案

【问题标题】：alternative to passing large matrices in R在 R 中传递大型矩阵的替代方法
【发布时间】：2013-05-12 01:35:11
【问题描述】：

Using Globals instead of passing large arrays in Matlab

我的问题与上面的问题完全相同，但我想要关于 R 的答案。我现在正在函数之间传递巨大的矩阵。矩阵中的数据在这些函数中没有改变。我只是使用矩阵。我的代码运行缓慢。我想知道是否有替代方法，例如使用全局变量或面向对象的方式。谢谢

【问题讨论】：

虽然从概念上讲参数被复制到函数中，但实际上底层实现足够智能，不会复制未更改的参数，因此不会降低性能。

【解决方案1】：

R 具有传递引用（有点）。当您将一个对象分配给另一个变量或传递给一个函数时，会创建另一个引用。但是，如果您通过其中一个引用来修改对象，则会生成实际副本。

f <- function(m) {
  .Internal(inspect(m))
}
g <- function(m) {
  m[1] <- 0
  .Internal(inspect(m))
}

m <- matrix(1,1)
.Internal(inspect(m))
## @452e308 14 REALSXP g0c1 [NAM(2),ATT] (len=1, tl=0) 1
## ATTRIB:
##   @42c8ee8 02 LISTSXP g0c0 [] 
##     TAG: @2faaf98 01 SYMSXP g0c0 [MARK,LCK,gp=0x4000] "dim" (has value)
##     @452e2d8 13 INTSXP g0c1 [NAM(2)] (len=2, tl=0) 1,1

# f shows that this is the same object (@452e308):
f(m)
## @452e308 14 REALSXP g0c1 [NAM(2),ATT] (len=1, tl=0) 1
## ATTRIB:
##   @42c8ee8 02 LISTSXP g0c0 [] 
##     TAG: @2faaf98 01 SYMSXP g0c0 [MARK,LCK,gp=0x4000] "dim" (has value)
##     @452e2d8 13 INTSXP g0c1 [NAM(2)] (len=2, tl=0) 1,1

# g shows a newly allocated object (@3941998):
g(m)
## @3941998 14 REALSXP g0c1 [NAM(1),ATT] (len=1, tl=0) 0
## ATTRIB:
##   @3b9fc80 02 LISTSXP g0c0 [] 
##     TAG: @2faaf98 01 SYMSXP g1c0 [MARK,LCK,gp=0x4000] "dim" (has value)
##     @3941ae8 13 INTSXP g0c1 [NAM(2)] (len=2, tl=0) 1,1

【讨论】：

只有在你完成作业之前才是正确的。然后复制开始。
@DWin 我以为我是这么说的。
我猜是这样，而且我的陈述并不是对 pass-by-promise 的特别好的总结。我被代码误导以为你在说 R 使用了传递引用。

【解决方案2】：

R 中没有“全局变量”。如果您想使用传递引用语义，您可以使用可能有用的 'data.table' 包，或者您可以使用环境或参考类。目前这个问题太模糊，无法透露更多细节。

【讨论】：