Sangyoon Lee

Background

While thinking about ways to represent incoming and outgoing flows in a business process, I thought about using export-import charts like the one shown here in the Learning R blog. However, as the author acknowledges, it is difficult to compare individual values using these charts. Regardless, I still wanted to have this graph for an at-a-glance view before breaking it into facets and comparing individual values.

My Solution

In the Learning R article, the author chooses to show multiple categories of import and export using stacked bars. Instead of representing multiple categories, I decided to use the color intensity on the bars’ fill as visual reinforcement of information the graph already contains. Import and export are represented by red and blue, respectively, and the transparency facilitates the visual comparison the reader must make between bars that are not side by side.

In the below example, I use the same subset of data as in the motivating post. Please refer to the linked article for the data used in this example. Make sure to click on the link that says “Access the subset used in this post in here.” rather than going to the Eurostat website. Save the file as “trade.csv” in the working directory. These are monthly trade data for the 27 European Union countries by broad economic categories (BEC) in millions of Euros.

First, load the necessary packages.

library(ggplot2)
library(plyr)
library(reshape)
library(scales)

For convenient and powerful data manipulation, plyr and reshape provide functions like ddply and melt. A relatively new package, scales is required for scale functions to format the numbers to specific scales within ggplot2.

Next, import the data, calculate the trade balance (export – import), and melt the data for ggplot2.

trade <- read.csv("trade.csv", header = TRUE, 
   stringsAsFactors = FALSE)
balance <- ddply(trade, .(Time), summarise, 
   balance = sum(EXP - IMP))
trade.m <- melt(trade, id.vars = c("BEC","Time"))

After the melt step, add another line to aggregate over BEC. This will further simplify the structure.

trade.a <- ddply(trade.m, c("Time", "variable"), summarise,
                 value = sum(value))

At this point, the data will look like this:

> head(trade.a)
     Time variable    value
1 2008M05      EXP 273153.2
2 2008M05      IMP 260789.1
3 2008M06      EXP 284994.7
4 2008M06      IMP 273033.0
5 2008M07      EXP 284681.6
6 2008M07      IMP 271122.2

We step through one layer at a time.

Layer 1: Start with export bars. We will add import data on the bottom of this graph.

ggplot(trade.a, aes(x=Time)) + 
   geom_bar(data = subset(trade.a, variable == "EXP"),
   aes(y=value, fill = value), stat = "identity")

Layer 2: Add the import data and attach it back-to-back to the export data. Label the x-axis and the y-axis accordingly.

last_plot() + geom_bar(data = subset(trade.a, variable == "IMP"), 
  aes(y=-value, fill = -value), stat = 'identity') +
  scale_y_continuous(labels = comma) + xlab("") + 
  ylab("Export  -  Import") + 
  scale_fill_gradient2(low = muted("red"),
  mid = "white", high = muted("blue"), midpoint = 0,space = "rgb")

Layer 3. Now add the balance trend line, remove the meaningless legend, and format the y-axis with commas.

last_plot() + geom_line(data = balance, aes(Time, balance, group = 1), size = 1) + geom_hline(yintercept = 0,colour = "grey90") + opts(legend.position = "none") 

Layer 4: Finally, change the x-axis to make it easy for viewers to read. The following result is my final product.