[tor-commits] [metrics-tasks/master] Add R code for #2686.

[karsten at torproject.org]

commit 3b7c299d3e869d92af69d01f3410dfd566ccc9d8
Author: Karsten Loesing <karsten.loesing at gmx.net>
Date:   Thu Apr 28 10:45:36 2011 +0200

    Add R code for #2686.
---
 task-2686/.gitignore          |    2 +
 task-2686/README              |    3 +
 task-2686/cbt-cutoff-guards.R |   88 +++++++++++++++++++++++++++++++++++++++++
 3 files changed, 93 insertions(+), 0 deletions(-)

diff --git a/task-2686/.gitignore b/task-2686/.gitignore
new file mode 100644
index 0000000..29a7166
--- /dev/null
+++ b/task-2686/.gitignore
@@ -0,0 +1,2 @@
+*.csv
+
diff --git a/task-2686/README b/task-2686/README
new file mode 100644
index 0000000..cf0ccf5
--- /dev/null
+++ b/task-2686/README
@@ -0,0 +1,3 @@
+$ R --slave -f filter.R --args *.mergedata
+$ R --slave -f cbt-cutoff-guards.R
+
diff --git a/task-2686/cbt-cutoff-guards.R b/task-2686/cbt-cutoff-guards.R
new file mode 100644
index 0000000..08a913c
--- /dev/null
+++ b/task-2686/cbt-cutoff-guards.R
@@ -0,0 +1,88 @@
+# Load ggplot library without printing out stupid warnings
+options(warn = -1)
+suppressPackageStartupMessages(library("ggplot2"))
+
+# Read data
+cbt <- read.csv("filtered.csv", stringsAsFactors = FALSE)
+
+# Transform data to put it in an ECDF.
+transform <- function(data, filenamepart, guards, cbtcutoff, filesize) {
+
+  # Remove unfinished or timed out downloads
+  data <- data[data$timeout == FALSE & data$failure == FALSE, ]
+
+  # Remove all other cbt cutoffs and file sizes
+  data <- data[data$guards == filenamepart & data$filesize == filesize, ]
+
+  # Order data frame by completion time
+  data <- data[order(data$completemillis), ]
+
+  # Convert to a data frame that has ordered completion times on the x
+  # axis and CDF value on the y axis, plus source and filesize for
+  # distinguishing these values from the other data frames we're going to
+  # bind together.
+  data <- data.frame(
+          x = data$completemillis,
+          y = (1:length(data$completemillis)) / length(data$completemillis),
+          guards = guards,
+          cbtcutoff = cbtcutoff,
+          filesize = filesize,
+          stringsAsFactors = FALSE)
+
+  # Return the transformed data frame
+  data
+}
+
+data <- rbind(
+  transform(cbt, "slow50cbt", "slow", "50", "50kb"),
+  transform(cbt, "slow60cbt", "slow", "60", "50kb"),
+  transform(cbt, "slow70cbt", "slow", "70", "50kb"),
+  transform(cbt, "slow80cbt", "slow", "80", "50kb"),
+  transform(cbt, "slow99cbt", "slow", "99", "50kb"),
+  transform(cbt, "slow50cbt", "slow", "50", "1mb"),
+  transform(cbt, "slow60cbt", "slow", "60", "1mb"),
+  transform(cbt, "slow70cbt", "slow", "70", "1mb"),
+  transform(cbt, "slow80cbt", "slow", "80", "1mb"),
+  transform(cbt, "slow99cbt", "slow", "99", "1mb"),
+  transform(cbt, "fast50cbt", "fast", "50", "50kb"),
+  transform(cbt, "fast60cbt", "fast", "60", "50kb"),
+  transform(cbt, "fast70cbt", "fast", "70", "50kb"),
+  transform(cbt, "fast80cbt", "fast", "80", "50kb"),
+  transform(cbt, "fast99cbt", "fast", "99", "50kb"),
+  transform(cbt, "fast50cbt", "fast", "50", "1mb"),
+  transform(cbt, "fast60cbt", "fast", "60", "1mb"),
+  transform(cbt, "fast70cbt", "fast", "70", "1mb"),
+  transform(cbt, "fast80cbt", "fast", "80", "1mb"),
+  transform(cbt, "fast99cbt", "fast", "99", "1mb"),
+  transform(cbt, "regular50cbt", "regular", "50", "50kb"),
+  transform(cbt, "regular60cbt", "regular", "60", "50kb"),
+  transform(cbt, "regular70cbt", "regular", "70", "50kb"),
+  transform(cbt, "regular80cbt", "regular", "80", "50kb"),
+  transform(cbt, "regular99cbt", "regular", "99", "50kb"),
+  transform(cbt, "regular50cbt", "regular", "50", "1mb"),
+  transform(cbt, "regular60cbt", "regular", "60", "1mb"),
+  transform(cbt, "regular70cbt", "regular", "70", "1mb"),
+  transform(cbt, "regular80cbt", "regular", "80", "1mb"),
+  transform(cbt, "regular99cbt", "regular", "99", "1mb"))
+
+data <- data[data$y <= 0.9975, ]
+data <- data[(data$x <= 120000 & data$filesize == "50kb") |
+             (data$x <= 600000 & data$filesize == "1mb"), ]
+
+data[data$guards == "slow", "guards"] <- "a) slowest overall"
+data[data$guards == "regular", "guards"] <- "c) default"
+data[data$guards == "fast", "guards"] <- "e) fastest overall"
+data[data$filesize == "50kb", "filesize"] <- "a) 50 KB"
+data[data$filesize == "1mb", "filesize"] <- "b) 1 MB"
+
+ggplot(data, aes(x = x / 1000, y = y, colour = cbtcutoff)) +
+geom_line(size = 0.5) +
+facet_grid(guards ~ filesize, scale = "free_x") +
+scale_x_continuous("\nDownload completion time in seconds") +
+scale_y_continuous("Fraction of completed downloads\n",
+  limits = c(0.85, 1), breaks = seq(0.9, 1, 0.05), formatter = "percent") +
+scale_colour_manual("CBT cutoff", breaks = c("50", "60", "70", "80", "99"),
+  values = alpha("black", c(0.2, 0.35, 0.5, 0.7, 1))) +
+opts(title = "Influence of guard nodes and circuit build timeouts on worst case performance\n")
+ggsave("cbt-cutoff.png", width = 8, height = 6, dpi = 150)
+