r - Rと恥ずかしいほど並列のデータを分割しますか？

Question

R を使用して並行して作業したい大きな RDS ファイルがあります。ファイルをロードすると、7.3 GB の RAM が必要になります。

多くのコアを使用しようとすると、R がメモリ不足でクラッシュします。オブジェクトのコピーを作成する代わりに共有メモリを使用するように mclapply に指示する方法はありますか?

これは私が持っているコードです:

results <- readRDS('ResultsICC.RDS')

rand <- 0
Icc <- c(.5, 1, 1.5)
n <- c(.1, .5, 1)
phi <- c(0, .5, 1)
parameterSettings <- expand.grid(rand=rand, Icc=Icc, n=n, phi=phi)

rr <- list()
Ns <- results[[1]][[1]][[2]][,c('Country', 'n')]
EstimatedBestPFiveArmRaw <- matrix(NA, 26, 1000)
EstimatedBestP <- matrix(NA, 26, 1000)


outterloop <- function(dataIN){
  for(k in 1:1000){ #1000
    best <- dataIN[[k]][[2]]
    EstimatedBestPFiveArmRaw[,k] <- rep(weighted.mean(best$estimatedBestPFiveArmRaw, best$n), 26) 
    pHat <- dataIN[[k]][[3]]  
    best <- Ns
    best$estimatedBest <- best$estimateBestP <- NA
    for(j in 1:26){ #26
      best$estimatedBest <- sapply(split(pHat[,paste0('cohort', j+1, 'pHat')], pHat$Country), 
                                   which.max)
      for(i in 1:nrow(best)) #nrow(best)
        best$estimatedBestP[i] <- pHat$p[pHat$Country==best$Country[i] &
                                           pHat$treatNum==best$estimatedBest[i]]  
      EstimatedBestP[j, k] <- weighted.mean(best$estimatedBestP, best$n)
    }
    rr <- (EstimatedBestP/EstimatedBestPFiveArmRaw-1)*100
  }
  return(rr)
}
library(parallel)
rr <- mclapply(X = results, FUN = outterloop, mc.cores = 27, mc.preschedule = T)

これを 32 コアと 64GB の RAM を搭載した Linux ボックスで実行しています。

ありがとう！