Ниже вы можете найти фрагмент кода на R, который я хотел бы преобразовать для запуска в виде параллельного процесса с использованием нескольких процессоров. Я пытался использовать пакет foreach, но далеко не продвинулся. Я не смог найти хороший пример того, как заставить его работать, учитывая тот факт, что у меня есть 3-уровневый вложенный цикл. Помощь будет чрезвычайно признательна. Пример кода ниже - я сделал простую функцию, чтобы она могла служить примером:
celnum <- c(10,20,30)
t2 <- c(1,2,3)
allrepeat <- 10
samplefunction <- function(celnum,t2){
x <- rnorm(100,celnum,t2)
y = sample(x, 1)
z = sample(x,1)
result = y+z
result
}
Получение результатов традиционным способом:
z_grid <- matrix(, nrow = length(celnum), ncol = length(t2))
repetitions <- matrix(, nrow = allrepeat, ncol = 1)
set.seed=20
for(i in 1:length(celnum)){
for (j in 1:length(t2)){
for (k in 1:allrepeat) {
results <- samplefunction(celnum[i],t2[j])
repetitions[k] <- results
z_grid[i,j] <- mean(repetitions,na.rm=TRUE)
}
}
}
z_grid
Теперь пытаемся сделать то же самое с помощью foreach:
set.seed=20
library(foreach)
library(doSNOW)
cl <- makeCluster(3, type = "SOCK")
registerDoSNOW(cl)
set.seed=20
output <- foreach(i=1:length(celnum),.combine='cbind' ) %:%
foreach (j=1:length(t2), .combine='c') %:%
foreach (k = 1:allrepeat) %do% {
mean(samplefunction(celnum[i],t2[j]) )
}
output
Это не работает так, как хотелось бы, так как возвращает матрицу размером 30x2 вместо 3x3. Я намерен смоделировать сценарий для комбинаций i и j k раз и хотел бы получить среднее значение этих k симуляций для каждой комбинации i и j.
