Create a regular sequence of date-times (POSIXct) using seq()

Question

My goal is to create a vector of POSIXct time stamps given a start, an end and a delta (15min, 1hour, 1day). I hoped I could use seq for this, but I have a problem converting between the numeric and POSIXct representation:

now <- Sys.time()
now
# [1] "2012-01-19 10:30:39 CET"
as.POSIXct(as.double(now), origin="1970-01-01", tz="CET")
# [1] "2012-01-19 09:30:39 CET"
as.POSIXct(as.double(now), origin=as.POSIXct("1970-01-01", tz="CET"), tz="CET")
# [1] "2012-01-19 09:30:39 CET"

One hour gets lost during this conversion. What am I doing wrong?

Gavin Simpson · Accepted Answer · 2016-02-11T16:48:39.100

27

There is a seq() method for objects of class "POSIXt" which is the super class of the "POSIXlt" and "POSIXct" classes. As such you don't need to do any conversion.

> now <- Sys.time()
> tseq <- seq(from = now, length.out = 100, by = "mins")
> length(tseq)
[1] 100
> head(tseq)
[1] "2012-01-19 10:52:38 GMT" "2012-01-19 10:53:38 GMT"
[3] "2012-01-19 10:54:38 GMT" "2012-01-19 10:55:38 GMT"
[5] "2012-01-19 10:56:38 GMT" "2012-01-19 10:57:38 GMT"

edited Feb 11 '16 at 16:48

answered Jan 19 '12 at 10:53

Gavin Simpson

170,508
25
396
453

3

And this is pointed out in the documentation of seq(), where I should have looked in the first place. Thank you! – Karsten W. Jan 19 '12 at 11:08
see @RichieCotton's answer below if you need it by a regular interval like "10 min" – d8aninja Jun 09 '17 at 16:55

score 10 · Answer 2 · answered Jan 19 '12 at 12:11

You have to be aware that when converting from POSIXct to numeric, R takes the timezone into account but always starts counting from a GMT origin :

> xgmt <- as.POSIXct('2011-01-01 14:00:00',tz='GMT')
> xest <- as.POSIXct('2011-01-01 14:00:00',tz='EST')
> (as.numeric(xgmt) - as.numeric(xest)) / 3600
[1] -5

As you see, the time in EST is conceived to be five hours earlier than the time in GMT, which is the time difference between both timezones. It's that value that is saved internally.

The as.POSIXCT() function just adds an attribute containing the timezone. It doesn't alter the value, so you get the time presented in GMT time, but with an attribute telling it is EST. This also means that once you go from POSIXct to numeric, you should treat your data as if it's GMT time. (It's a whole lot more complex than that, but it's the general idea). So you have to calculate the offset as follows:

> nest <- as.numeric(xest)
> origin <- as.POSIXct('1970-01-01 00:00:00',tz='EST')
> offset <- as.numeric(origin)
> as.POSIXct(nest-offset,origin=origin)
[1] "2011-01-01 14:00:00 EST"

This works whatever the timezone is in your locale (in my case, that's actually CET). Also note that behaviour of timezone data can differ between systems.

Richie Cotton · Answer 3 · 2012-01-19T10:57:52.593

9

These time zone issues are always fiddly, but I think the problem is that your origin is being calculated in the wrong time zone (since the string only specifies the date).

Try using origin <- now - as.numeric(now).

Alternatively, use lubridate::origin, which is the string "1970-01-01 UTC".

A full solution, again using lubridate.

start <- now()
seq(start, start + days(3), by = "15 min")

edited Jan 19 '12 at 10:57

answered Jan 19 '12 at 10:45

Richie Cotton

118,240
47
247
360

You have to be careful though, as this only works when converting to your own current locale. A more general solution is given in my answer. The idea is the same, but is independent of your own locale. – Joris Meys Jan 19 '12 at 13:01
1

@Richie Using the `lubridate` example, I am getting `invalid string for "by"` – Brian Feb 22 '17 at 22:50
1

@BrianVanover Double check your typing, and if you still have a problem then post a new question. – Richie Cotton Feb 23 '17 at 13:43
@RichieCotton bet you he wrote 'minutes' and not 'min'. I say that because that's exactly what I did. – d8aninja Jun 09 '17 at 16:54

score 4 · Answer 4 · answered Jan 19 '12 at 10:18

I do not have an answer to your problem, but I do have an alternative way of creating vectors of POSIXct objects. If, for example, you want to create a vector of 1000 timestamps from now with a delta_t of 15 minutes:

now = Sys.time()
dt = 15 * 60 # in seconds
timestamps = now + seq(0, 1000) * dt
> head(timestamps)
[1] "2012-01-19 11:17:46 CET" "2012-01-19 11:32:46 CET"
[3] "2012-01-19 11:47:46 CET" "2012-01-19 12:02:46 CET"
[5] "2012-01-19 12:17:46 CET" "2012-01-19 12:32:46 CET"

The trick is you can add a vector of seconds to a POSIXct object.

score 4 · Answer 5 · answered Jan 19 '12 at 14:10

An alternative to using seq.POSIXt is xts::timeBasedSeq, which allows you to specify the sequence as a string:

library(xts)
now <- Sys.time()
timeBasedSeq(paste("2012-01-01/",format(now),"/H",sep=""))  # Hourly steps
timeBasedSeq(paste("2012-01-01/",format(now),"/d",sep=""))  # Daily steps

score 1 · Answer 6 · answered Nov 09 '15 at 15:13

1

You need to use seq(from=start,to=end, by=step). Note that in step you can either use "days" or an integer defining how many seconds elapse from item to item.

answered Nov 09 '15 at 15:13

Marcelo Bielsa

142
11

You need to read the other answers before posting a new one 3 years later. This was already provided in at least two earlier answers. – David Arenburg Nov 09 '15 at 15:18

Create a regular sequence of date-times (POSIXct) using seq()

6 Answers6

Linked

Related