web-dev-qa-db-ja.com

Rでミリ秒単位の実行時間を取得します

私はtic(), toc() functionsを使用してこれに対する解決策を読みました

tic <- function(gcFirst = TRUE, type=c("elapsed", "user.self", "sys.self"))
{
   type <- match.arg(type)
   assign(".type", type, envir=baseenv())
   if(gcFirst) gc(FALSE)
   tic <- proc.time()[type]         
   assign(".tic", tic, envir=baseenv())
   invisible(tic)
}

toc <- function()
{
   type <- get(".type", envir=baseenv())
   toc <- proc.time()[type]
   tic <- get(".tic", envir=baseenv())
   print(toc - tic)
   invisible(toc)
}




tic();
-----code----
toc();


elapsed 
   0.15

しかし、ミリ秒単位で多くの精度を取得したいと思いますか？

また、私はこれを使用していました

ptm <- proc.time()
---code
proc.time() - ptm

そしてこれを入手

   user  system elapsed 
   1.55    0.25    1.84

より多くの小数またはより高い精度を取得するにはどうすればよいですか？

rbenchmarkingexecution-timemilliseconds
22
cMinor

1）タイミングはオペレーティングシステムに依存します。 Windowsでは、ミリ秒しか取得できない場合があります。

2）tic()とtoc()を定義する必要はありません。Rにはsystem.time()があります。次に例を示します。

R> system.time(replicate(100, sqrt(seq(1.0, 1.0e6))))
   user  system elapsed 
  2.210   0.650   2.867 
R>

3）優れたアドオンパッケージ rbenchmark および microbenchmark があります。

3.1） rbenchmark はコマンドの比較に特に役立ちますが、直接使用することもできます。

R> library(rbenchmark)
R> x <- seq(1.0, 1.0e6); benchmark(sqrt(x), log(x))
     test replications elapsed relative user.self sys.self user.child sys.child
2  log(x)          100   5.408  2.85835      5.21     0.19          0         0
1 sqrt(x)          100   1.892  1.00000      1.62     0.26          0         0
R>

3.2） マイクロベンチマーク 最高精度の測定に優れています：

R> library(microbenchmark)
R> x <- seq(1.0, 1.0e6); microbenchmark(sqrt(x), log(x))
Unit: nanoseconds
     expr      min       lq   median       uq      max
1  log(x) 50589289 50703132 55283301 55353594 55917216
2 sqrt(x) 15309426 15412135 15452990 20011418 39551819
R>

そして最後の、特にLinuxでは、すでにナノ秒を与えています。結果などをプロットすることもできるので、そのパッケージを詳しく見てみましょう。

34
Dirk Eddelbuettel

これは良いです：

options(digits.secs = 6) # This is set so that milliseconds are displayed

start.time <- Sys.time()

...Relevant code...

end.time <- Sys.time()
time.taken <- end.time - start.time
time.taken

ここ から取得。

6
Denis Rasulev

場所 start_timeコードの前とend_timeコードの後。

つまり.

start_time <- as.numeric(as.numeric(Sys.time())*1000, digits=15) # place at start

-----code----

end_time <- as.numeric(as.numeric(Sys.time())*1000, digits=15) # place at end

end_time - start_time    # run time (in milliseconds)
0
stevec