sum_thread_iops.r

#!/usr/local/bin/Rscript

filenames <- list.files("./",pattern="output*",full.names=TRUE)
files <- lapply(filenames, function(filename){ read.csv(file=filename,sep=';',stringsAsFactors=FALSE,skip=4) })
post_ramp_subsets <- lapply(files, function(file) { subset(file,file[50]>=60000) })
results <- lapply(post_ramp_subsets, function(thread_subset) {
  threads <- unique(thread_subset[,3])
  thread_subsets <- lapply(threads, function(thread){ subset(thread_subset,thread_subset[,3]==thread) })
  means <- lapply(thread_subsets, function(thread_subset){ mean(thread_subset[,49]) })
  thread_sum <- Reduce("+",means)
})
Reduce("+",Filter(Negate(is.null),results))