N <- 1000000 Rprof("profile_out.txt") # one-at-a-time version draws <- matrix(NA,N,1) for (t in 0:N) { draws[t,1] <- rnorm(1) } # all at once version #draws <- rnorm(N) Rprof(NULL) summaryRprof("profile_out.txt")$by.self