7

我的代码如下所示（与原始版本相比，它有点简化，但仍然反映了问题）。

require(VGAM)

Median.sum  = vector(mode="numeric", length=75) 
AA.sum      = vector(mode="numeric", length=75)                                                    
BB.sum      = vector(mode="numeric", length=75)                   
Median      = array(0, dim=c(75 ,3)) 
AA          = array(0, dim=c(75 ,3))                                                    
BB          = array(0, dim=c(75 ,3))                              

y.sum     = vector(mode="numeric", length=100000)
y         = array(0, dim=c(100000,3))
b.size    = vector(mode="numeric", length=3) 
c.size    = vector(mode="numeric", length=3) 


for (h in 1:40)
{
  for (j in 1:75)
  {  
    for (i in 1:100000)
    {
      y.sum[i] = 0

      for (f in 1:3)
      {
        b.size[f] = rbinom(1, 30, 0.9)
        c.size[f] = 30 - rbinom(1, 30, 0.9) + 1
        y[i, f] = sum( rlnorm(b.size[f], 8.5, 1.9) ) + 
          sum( rgpd(c.size[f], 120000, 1870000, 0.158) )
        y.sum[i] = y.sum[i] + y[i, f]
      }
    }

    Median.sum[j] = median(y.sum)
    AA.sum[j] = mean(y.sum)
    BB.sum[j] = quantile(y.sum, probs=0.85)

    for (f in 1:3)
    {
      Median[j,f] = median(y[,f])
      AA[j,f] = mean(y[,f])
      BB[j,f] = quantile(y[,f], probs=0.85)
    }
  }
  #gc()
}

它在执行过程中中断 (h=7, j=1, i=93065) 并出现错误：

Error: cannot allocate vector of size 526.2 Mb

收到此消息后，我已经阅读了 this，this和this，但这还不够。问题是，垃圾收集器 (gc()) 和清除工作区中的所有对象都没有帮助。我的意思是我已经尝试将我的代码都放入我的代码中：垃圾收集器和操作删除所有变量并在循环中再次声明它们（看看 #gc() 所在的位置 - 但后者不包括在内在我发布的代码中）。

对我来说这似乎很奇怪，因为所有过程在循环的每个步骤中都使用相同的对象（=> 并且应该在循环的每个步骤中消耗相同的内存量）。为什么内存消耗会随着时间的推移而增加？

更糟糕的是，如果我想在 R 的同一会话中工作甚至执行：

rm(list=ls())
gc()

即使我想声明一些次要的内容，我仍然会收到相同的错误消息：

abc = array(0, dim=c(10,3))

只有关闭 R 并开始新的会话有帮助。为什么？也许有一些方法可以重新编码我的循环？

R：2.15.1（32 位），操作系统：Windows XP（32 位）

我在这里很新，所以每一个小费都很感激！提前致谢。

编辑：（来自阿伦）。我发现仅通过一个简单的示例就可以更轻松地重现此行为。启动一个新的 R 会话并复制并粘贴此代码，然后观察系统监视器中的内存增长情况。

mm <- rep(0, 1e4) # initialise a vector
for (i in 1:1e3) {
    for (j in 1:1e3) {
        for (k in 1:1e4) {
            mm[k] <- k # already pre-allocated
         }
    }
}
4

2 回答 2

4

gc()在for (i in 1:100000)循环内添加调用。

在 Arun 代码的紧密循环中添加对 的调用可以gc()消除其内存增长。

这显示了内存增长：

mm <- rep(0, 1e4) # initialise a vector
for (i in 1:1e3) {
    for (j in 1:1e3) {
        for (k in 1:1e4) {
            mm[k] <- k # already pre-allocated
         }
     }
 }

这不会：

mm <- rep(0, 1e4) # initialise a vector
for (i in 1:1e3) {
    for (j in 1:1e3) {
        for (k in 1:1e4) {
            mm[k] <- k # already pre-allocated
            gc()
         }
     }
 }

这里的自动垃圾收集有些问题。如前所述，收集器在第一种情况下被调用gcinfo(TRUE)。但是记忆增长得非常快。

于 2013-03-24T17:36:35.427 回答
2

这似乎可行（将最内层循环放入函数中）。我直到最后才运行它，因为它很慢，但我没有注意到你的代码中的内存膨胀。

require(VGAM)

Median.sum  = vector(mode="numeric", length=75) 
AA.sum      = vector(mode="numeric", length=75)                                                    
BB.sum      = vector(mode="numeric", length=75)                   
Median      = array(0, dim=c(75 ,3)) 
AA          = array(0, dim=c(75 ,3))                                                    
BB          = array(0, dim=c(75 ,3))                              


inner.fun <- function() {
  y.sum     = vector(mode="numeric", length=100000)
  y         = array(0, dim=c(100000,3))
  b.size    = vector(mode="numeric", length=3) 
  c.size    = vector(mode="numeric", length=3) 
  for (i in 1:100000)
    {
      y.sum[i] = 0

      for (f in 1:3)
      {
        b.size[f] = rbinom(1, 30, 0.9)
        c.size[f] = 30 - rbinom(1, 30, 0.9) + 1
        y[i, f] = sum( rlnorm(b.size[f], 8.5, 1.9) ) + 
          sum( rgpd(c.size[f], 120000, 1870000, 0.158) )
        y.sum[i] = y.sum[i] + y[i, f]
      }
    }
    list(y.sum, y)
}

for (h in 1:40)
{
  cat("\nh =", h,"; j = ")
  for (j in 1:75)
  {  
    cat(j," ")
    result = inner.fun()
    y.sum = result[[1]]
    y = result[[2]]
    Median.sum[j] = median(y.sum)
    AA.sum[j] = mean(y.sum)
    BB.sum[j] = quantile(y.sum, probs=0.85)

    for (f in 1:3)
    {
      Median[j,f] = median(y[,f])
      AA[j,f] = mean(y[,f])
      BB[j,f] = quantile(y[,f], probs=0.85)
    }
  }
}
于 2013-03-25T15:34:41.467 回答