Sequential for loop benchmark V2

args <- commandArgs(trailingOnly=TRUE)
n <- if (length(args) > 0) as.integer(args[1]) else 1000000
m <- if (length(args) > 1) as.integer(args[2]) else ceiling(n / 100)
set.seed(107)
td <- data.frame(val=rnorm(n), id=sample(m, n, replace=TRUE))

start <- proc.time()[3]
res <- rep(0, NROW(td))
for (i in unique(td$id))
  res[td$id == i] <- mean(td$val[td$id != i])
elapsed <- proc.time()[3] - start

cat("Sequential for loop:\n")
cat(sprintf("Rows: %d, Unique IDs: %d\n", length(td$id), length(unique(td$id))))
cat(sprintf("Elapsed time: %f\n", elapsed))

library(digest)
cat(sprintf("MD5 hash: %s\n", digest(res)))