É muito frequente ouvir (ou ler) que os loops não são eficientes no R
e devem ser evitados (neste link ou noutro link ou mesmo neste).
E provar esta afirmação é simples:
numeros <- rnorm(10000)
com_loop <- function(vetor) {
res <- 0
for (i in seq_along(vetor)) {
res <- res + vetor[i]
}
res
}
microbenchmark::microbenchmark(
loop = com_loop(numeros),
vetorizado = sum(numeros)
)
Unit: microseconds
expr min lq mean median uq max neval
loop 494.709 512.670 562.71062 514.723 551.9285 3074.480 100
vetorizado 9.750 10.263 10.77702 10.264 10.2640 28.226 100
As perguntas que faço são:
- Por que os loops são lentos no
R
r? - Quais alternativas exitemexistem? (pacotes, estratégias, etc)