He leído una solución para esto usandotic(), 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 

¿Pero me gustaría obtener mucha precisión en milisegundos?

También estaba usando este

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

y obtener esto

   user  system elapsed 
   1.55    0.25    1.84 

¿Cómo obtener más decimales o más precisión?