commit ba24732e4d23671da527b45d0fb8e35c54470f95
Author: Karsten Loesing <karsten.loesing(a)gmx.net>
Date: Wed Jul 25 12:31:38 2012 +0200
Update #4255 graphing code to ggplot2 0.9.0.
---
task-4255/stability.R | 35 +++++++++++++++++++----------------
1 files changed, 19 insertions(+), 16 deletions(-)
diff --git a/task-4255/stability.R b/task-4255/stability.R
index c46266e..b4f162e 100644
--- a/task-4255/stability.R
+++ b/task-4255/stability.R
@@ -1,4 +1,7 @@
library(ggplot2)
+library(scales)
+library(grid)
+library(reshape)
stability <- read.csv("stability.csv", stringsAsFactors = FALSE)
d <- stability[stability$time > '2010-07' & stability$time < '2011-07', ]
@@ -18,8 +21,8 @@ d <- rbind(d,
y = NA, ymin = NA, ymax = NA))
ggplot(d, aes(x = as.Date(x), y = y, ymin = ymin, ymax = ymax)) +
geom_line() +
-scale_x_date(name = "", major = "3 months", minor = "1 month",
- format = "%b %Y") +
+scale_x_date("", breaks = "3 months", minor_breaks = "1 month",
+ labels = date_format("%b %Y")) +
scale_y_continuous(name = "Running \nbridges ",
limits = c(0, max(d_mean$running, na.rm = TRUE))) +
opts(axis.title.x = theme_text(size = 12 * 0.8, face = "bold",
@@ -35,8 +38,8 @@ d <- stability[stability$time > '2010-07-10' &
stability$time < '2010-07-31', ]
a <- ggplot(d, aes(x = as.POSIXct(time), y = running)) +
geom_point(size = 0.75) +
-scale_x_datetime(name = "", major = "1 week", minor = "1 day",
- format = "%b %d, %Y") +
+scale_x_datetime(name = "", breaks = "1 week", minor_breaks = "1 day",
+ labels = date_format("%b %d, %Y")) +
scale_y_continuous(name = "Running \nbridges ",
limits = c(0, max(d$running, na.rm = TRUE))) +
opts(axis.title.x = theme_text(size = 12 * 0.8, face = "bold",
@@ -53,8 +56,8 @@ d <- rbind(
data.frame(time = e$time, running = 687, colour = "grey"))
b <- ggplot(d, aes(x = as.POSIXct(time), y = running, colour = colour)) +
geom_point(size = 0.75) +
-scale_x_datetime(name = "", major = "1 week", minor = "1 day",
- format = "%b %d, %Y") +
+scale_x_datetime(name = "", breaks = "1 week", minor_breaks = "1 day",
+ labels = date_format("%b %d, %Y")) +
scale_y_continuous(name = "Running \nbridges ",
limits = c(0, max(d$running, na.rm = TRUE))) +
scale_colour_manual(values = c("black", "grey60")) +
@@ -108,8 +111,8 @@ geom_line() +#colour = "grey30") +
geom_hline(data = e, aes(yintercept = yintercept), colour = "gray40",
linetype = 2) +
facet_grid(var ~ ., scales = "free_y") +
-scale_x_date(name = "", major = "3 months", minor = "1 month",
- format = "%b %Y") +
+scale_x_date("", breaks = "3 months", minor_breaks = "1 month",
+ labels = date_format("%b %Y")) +
scale_y_continuous(name = "") +
opts(axis.title.x = theme_text(size = 12 * 0.8, face = "bold",
hjust = 0.5),
@@ -132,9 +135,9 @@ d <- melt(d, id = "date")
ggplot(d, aes(x = date, y = value / 10000, linetype = variable)) +
geom_line() +
scale_y_continuous(name = paste("10th perc. \nWFU in \n",
- "the future ", sep = ""), formatter = "percent", limits = c(0, 1)) +
-scale_x_date(name = "", major = "3 months", minor = "1 month",
- format = "%b %Y") +
+ "the future ", sep = ""), labels = percent, limits = c(0, 1)) +
+scale_x_date("", breaks = "3 months", minor_breaks = "1 month",
+ labels = date_format("%b %Y")) +
scale_linetype_manual(name = paste("Requirements for\nconsidering",
"a\nbridge as stable\n"), breaks = c("perc10wfu50wfu50wmtbac",
"perc10wfu50wfu0wmtbac", "perc10wfu0wfu50wmtbac",
@@ -166,8 +169,8 @@ scale_y_continuous(name = paste("10th perc. \ntime on \nthe same \n",
breaks = seq(0, max(d$value / 86400, na.rm = TRUE), 7),
minor = seq(0, max(d$value / 86400, na.rm = TRUE), 1),
limits = c(0, max(d$value / 86400, na.rm = TRUE))) +
-scale_x_date(name = "", major = "3 months", minor = "1 month",
- format = "%b %Y") +
+scale_x_date("", breaks = "3 months", minor_breaks = "1 month",
+ labels = date_format("%b %Y")) +
scale_linetype_manual(name = paste("Requirements for\nconsidering",
"a\nbridge as stable\n"), breaks = c("perc10tosa50wfu50wmtbac",
"perc10tosa0wfu50wmtbac", "perc10tosa50wfu0wmtbac",
@@ -207,9 +210,9 @@ d <- rbind(d,
ggplot(d, aes(x = x, y = y, linetype = variable)) +
geom_line() +
scale_y_continuous(name = "Fraction of \nRunning \nbridges ",
- formatter = "percent", limits = c(0, max(d$y, na.rm = TRUE))) +
-scale_x_date(name = "", major = "3 months", minor = "1 month",
- format = "%b %Y") +
+ labels = percent, limits = c(0, max(d$y, na.rm = TRUE))) +
+scale_x_date("", breaks = "3 months", minor_breaks = "1 month",
+ labels = date_format("%b %Y")) +
scale_linetype_manual(name = paste("\nRequirements for\nconsidering",
"a\nbridge as stable\n"), values = c(3, 2, 4)) +
opts(axis.title.x = theme_text(size = 12 * 0.8, face = "bold",