Placing the grid along date tickmarks

Question

I have the following data:

x=strptime(20010101:20010110)
y=1:10
z=data.frame(x,y)

So my data is this:

            x  y
1  2001-01-01  1
2  2001-01-02  2
3  2001-01-03  3
4  2001-01-04  4
5  2001-01-05  5
6  2001-01-06  6
7  2001-01-07  7
8  2001-01-08  8
9  2001-01-09  9
10 2001-01-10 10

When I create a plot in base using:

plot(x,y)
grid(NULL,NULL)

My vertical grid does not align with the date tick marks. I know this seems like a pretty simple problem, but I have not found a solution to this anywhere. Is there a way to get the vertical grid to align with the date tick marks using base that does not require me to do this:

abline(v=as.numeric(strptime(c(20010102,20010104,20010106,20010108,20010110),'%Y%m%d')))

I have a lot of plots with different dates and I would really like to automate this as much as possible, hopefully using base.

score 8 · Answer 1 · answered Feb 02 '12 at 21:22

8

The function axis draws your axes, tick marks and labels, and returns the tick mark positions as a vector.

Since you have Date data, you need to use axis.Date to do this, and then use abline to plot the grid:

z=data.frame(
  x=seq(as.Date("2001-01-01"), by="+1 month", length.out=10)
  y=1:10
)
plot(y~x, data=z)
abline(v=axis.Date(1, z$x), col="grey80")

enter image description here

answered Feb 02 '12 at 21:22

Andrie

176,377
47
447
496

Thanks for this alternative! I love how there are 80 different ways to do everything in R :) – thequerist Feb 02 '12 at 21:53
Comman missing in data frame assignment ... z=data.frame(x=seq(as.Date("2001-01-01"), by="+1 month", length.out=10), y=1:10) – Markm0705 Nov 25 '20 at 07:19

Josh O'Brien · Answer 2 · 2012-02-02T21:30:54.217

3

As the ?grid help file says, "if more fine tuning is required, [you can] use abline(h = ., v = .) directly". It's a little bit more work, but not much, and would be easy enough to wrap up in a function if you want to use it often.

Here's one way to do it:

plot(x,y)
abline(v = pretty(extendrange(z$x)), 
       h = pretty(extendrange(z$y)),
       col = 'lightgrey', lty = "dotted")
points(x,y, pch=16)

enter image description here

edited Feb 02 '12 at 21:30

answered Feb 02 '12 at 21:24

Josh O'Brien

159,210
26
366
455

Thanks, how did you get information on what extendedrange does? I wish I could accept more than one answer, I gave it to P Lapointe because he also added the function to his. – thequerist Feb 02 '12 at 21:52
No problem. Just check `?extendrange` to learn a bit more. It's a function used by plotting functions to determine the extents of the x- and y- axes. (Compare `plot(1:10)` and `plot(1:10, xaxs="i")` to see what I mean.) – Josh O'Brien Feb 02 '12 at 22:07
Thanks, I kept typing `?extend_ed_range` and was wondering why I was getting nowhere. – thequerist Feb 02 '12 at 22:14

score 3 · Accepted Answer · answered Feb 02 '12 at 21:39

abline can extract the date ticks from your POSIXlt vector (via strptime).

x=strptime(20010101:20010110,format="%Y%m%d")
y=1:10

plot(x,y)
grid(nx=NA, ny=NULL)
abline(v=axis.POSIXct(1, x=pretty(x)),col = "lightgray", lty = "dotted", lwd = par("lwd"))

I would suggest creating your own function which would add both horizontal and vertical grids.

my.grid <-function(){
grid(nx=NA, ny=NULL)
abline(v=axis.POSIXct(1, x=pretty(x)),col = "lightgray", lty = "dotted", lwd =
par("lwd"))
}

plot(x,y)
my.grid()

score 1 · Answer 4 · answered Feb 02 '12 at 20:52

1

I'm not familiar with the intricacies of the base plots. But, ggplot does this effectively.

x <- strptime(20010101:20010110, format='%Y%m%d')
y <- 1:10
z <- data.frame(x, y)
qplot(x,y,data=z,'point')

answered Feb 02 '12 at 20:52

Justin

42,475
9
93
111

Placing the grid along date tickmarks

4 Answers4

Linked

Related