¡Hola! Quede remetido con tu problema, aunque no soy un experto en R, hice las siguientes pruebas:
trials <- 10000
ptime <- system.time({
r <- foreach(icount(trials), .combine=cbind) %do% {
ind <- sample(100, 100, replace=TRUE)
}
})[3]
ptime
Me arroja 2.25s (do)
trials <- 10000
ptime <- system.time({
r <- foreach(icount(trials), .combine=cbind) %dopar% {
ind <- sample(100, 100, replace=TRUE)
}
})[3]
ptime
Me arroja 2.62s (dopar)
x <- iris[which(iris[,5] != "setosa"), c(1,5)]
trials <- 10000
ptime <- system.time({
r <- foreach(icount(trials), .combine=cbind) %do% {
ind <- sample(100, 100, replace=TRUE)
result1 <- glm(x[ind,2]~x[ind,1], family=binomial(logit))
coefficients(result1)
}
})[3]
ptime
Me arrojo 21.53s (do)
x <- iris[which(iris[,5] != "setosa"), c(1,5)]
trials <- 10000
ptime <- system.time({
r <- foreach(icount(trials), .combine=cbind) %dopar% {
ind <- sample(100, 100, replace=TRUE)
result1 <- glm(x[ind,2]~x[ind,1], family=binomial(logit))
coefficients(result1)
}
})[3]
ptime
Me arroja 11.03s (dopar)
No habiendo un uso excesivo de RAM ni CPU en ningún caso, lo único que se me ocurre es que al usar sample solo, se terminan gastando más recursos y tiempo en la asignación en múltiples CPUs, que si lo hiciera directamente un solo CPU. Probablemente sample se va ejecutando muy rápido y la asignación termina tomando más tiempo que cada ejecución, quizás incluso llegando a saturar al sistema en tu caso.
¡Saludos y suerte!