Saya membandingkan sample
fungsi dalam R dan membandingkannya dengan igraph:sample_seq
dan berlari ke hasil yang aneh.
Ketika saya menjalankan sesuatu seperti:
library(microbenchmark)
library(igraph)
set.seed(1234)
N <- 55^4
M <- 500
(mbm <- microbenchmark(v1 = {sample(N,M)},
v2 = {igraph::sample_seq(1,N,M)}, times=50))
Saya mendapatkan hasil seperti ini:
Unit: microseconds
expr min lq mean median uq max neval
v1 21551.475 22655.996 26966.22166 23748.2555 28340.974 47566.237 50
v2 32.873 37.952 82.85238 81.7675 96.141 358.277 50
Tetapi ketika saya berlari, misalnya,
set.seed(1234)
N <- 100^4
M <- 500
(mbm <- microbenchmark(v1 = {sample(N,M)},
v2 = {igraph::sample_seq(1,N,M)}, times=50))
Saya mendapatkan hasil yang jauh lebih cepat untuk sample
:
Unit: microseconds
expr min lq mean median uq max neval
v1 52.165 55.636 64.70412 58.2395 78.636 88.120 50
v2 39.174 43.504 62.09600 53.5715 73.253 176.419 50
Tampaknya kapan N
kekuatan 10 (atau nomor khusus lainnya?), sample
Jauh lebih cepat daripada yang lebih kecil lainnya N
yang bukan kekuatan 10. Apakah ini perilaku yang diharapkan atau apakah saya kehilangan sesuatu?