How to do a stacked area plot chart with a custom X-axis labeling

561 views Asked by At

I would like to do a stacked area plot chart with a customized x-axis. Specifically, I would like to label the tick marks on x-axis with cex.axis=0.8 (make it smaller to fit all the values) and I would like to use red color to highlight some of the x values (for example, weekend dates as contrasted with weekdays).

stackpoly in plotrix package does a nice job of creating the plot. But I am not able to control cex.axis or to suppress the plotting of x-axis.

Here are my two attempts to get this done:

## dates in winter of 2014 (these are the labels for my X-axis)
date.labels=paste( c(rep(1,2), rep(2,21)), c(30, 31, 1:21), sep = "/")
> date.labels
[1] "1/30" "1/31" "2/1"  "2/2"  "2/3"  "2/4"  "2/5"  "2/6"  "2/7"  "2/8"        "2/9"  "2/10" "2/11"
[14] "2/12" "2/13" "2/14" "2/15" "2/16" "2/17" "2/18" "2/19" "2/20" "2/21"

weekend=c(FALSE, FALSE, TRUE, TRUE, rep(FALSE, 5), TRUE, TRUE, rep(FALSE, 5), TRUE, TRUE, rep(FALSE, 5))

set.seed(10)
#This is the data I want to plot
fake.data=data.frame(abs(rnorm(23)), abs(rnorm(23)), abs(rnorm(23)))

library("plotrix")

#First attempt: use cex.axis in stackpoly to control font size. 
#overplot the weekend dates with red via axis()
stackpoly(fake.data, stack=TRUE, xat=1:23, xaxlab=date.labels, cex.axis=0.8,
      xlab="dates in 2014", sub="weekend dates are marked red")
#stackpoly ignores cex.axis=0.8

axis(1, (1:23)[weekend], labels= date.labels[ weekend], col.axis="red",     cex.axis=0.8)
#but axis does not ignore cex.axis=0.8

## Another attempt: suppress x axis altogether in stackpoly

stackpoly(fake.data, stack=TRUE, xlim=c(1,23), xaxt="n",
      xlab="dates in 2014", sub="weekend dates are marked red")
#stackpoly ignores xaxt="n"

axis(1, (1:23)[weekend], labels= date.labels[ weekend], col.axis="red",    cex.axis=0.8)
axis(1, (1:23)[!weekend], labels= date.labels [ !weekend], col.axis="black", cex.axis=0.8)
#but axis does not ignore cex.axis=0.8
1

There are 1 answers

0
jlhoward On

Try setting par(...) before the call to stackpoly(...)

par(mfrow=c(1,2))
# original axis text size
stackpoly(fake.data, stack=TRUE, xat=1:23, xaxlab=date.labels, 
          xlab="dates in 2014", sub="weekend dates are marked red")
# smaller axis text size
par(cex.axis=0.8)
stackpoly(fake.data, stack=TRUE, xat=1:23, xaxlab=date.labels, 
          xlab="dates in 2014", sub="weekend dates are marked red")