781 lines
28 KiB
Plaintext
781 lines
28 KiB
Plaintext
# Dates and times {#sec-dates-and-times}
|
|
|
|
```{r}
|
|
#| results: "asis"
|
|
#| echo: false
|
|
source("_common.R")
|
|
status("polishing")
|
|
|
|
# https://github.com/tidyverse/lubridate/issues/1058
|
|
options(warnPartialMatchArgs = FALSE)
|
|
```
|
|
|
|
## Introduction
|
|
|
|
This chapter will show you how to work with dates and times in R.
|
|
At first glance, dates and times seem simple.
|
|
You use them all the time in your regular life, and they don't seem to cause much confusion.
|
|
However, the more you learn about dates and times, the more complicated they seem to get!
|
|
|
|
To warm up think about how many days there are in a year, and how many hours there are in a day.
|
|
You probably remembered that most years have 365 days, but leap years have 366.
|
|
Do you know the full rule for determining if a year is a leap year[^datetimes-1]?
|
|
The number of hours in a day is a little less obvious: most days have 24 hours, but in places that use daylight saving time (DST), one day each year has 23 hours and another has 25.
|
|
|
|
[^datetimes-1]: A year is a leap year if it's divisible by 4, unless it's also divisible by 100, except if it's also divisible by 400.
|
|
In other words, in every set of 400 years, there's 97 leap years.
|
|
|
|
Dates and times are hard because they have to reconcile two physical phenomena (the rotation of the Earth and its orbit around the sun) with a whole raft of geopolitical phenomena including months, time zones, and DST.
|
|
This chapter won't teach you every last detail about dates and times, but it will give you a solid grounding of practical skills that will help you with common data analysis challenges.
|
|
|
|
We'll begin by showing you how to create date-times from various inputs, and then once you've got a date-time, how you can extract components like year, month, and day.
|
|
We'll then dive into the tricky topic of working with time spans, which come in a variety of flavors depending on what you're trying to do.
|
|
We'll conclude with a brief discussion of the additional challenges posed by time zones.
|
|
|
|
### Prerequisites
|
|
|
|
This chapter will focus on the **lubridate** package, which makes it easier to work with dates and times in R.
|
|
lubridate is not part of core tidyverse because you only need it when you're working with dates/times.
|
|
We will also need nycflights13 for practice data.
|
|
|
|
```{r}
|
|
#| message: false
|
|
library(tidyverse)
|
|
|
|
library(lubridate)
|
|
library(nycflights13)
|
|
```
|
|
|
|
## Creating date/times {#sec-creating-datetimes}
|
|
|
|
There are three types of date/time data that refer to an instant in time:
|
|
|
|
- A **date**.
|
|
Tibbles print this as `<date>`.
|
|
|
|
- A **time** within a day.
|
|
Tibbles print this as `<time>`.
|
|
|
|
- A **date-time** is a date plus a time: it uniquely identifies an instant in time (typically to the nearest second).
|
|
Tibbles print this as `<dttm>`.
|
|
Base R calls these POSIXct, but doesn't exactly trip off the tongue.
|
|
|
|
In this chapter we are going to focus on dates and date-times as R doesn't have a native class for storing times.
|
|
If you need one, you can use the **hms** package.
|
|
|
|
You should always use the simplest possible data type that works for your needs.
|
|
That means if you can use a date instead of a date-time, you should.
|
|
Date-times are substantially more complicated because of the need to handle time zones, which we'll come back to at the end of the chapter.
|
|
|
|
To get the current date or date-time you can use `today()` or `now()`:
|
|
|
|
```{r}
|
|
today()
|
|
now()
|
|
```
|
|
|
|
Otherwise, the following sections describe the four ways you're likely to create a date/time:
|
|
|
|
- While reading a file with readr.
|
|
- From a string.
|
|
- From individual date-time components.
|
|
- From an existing date/time object.
|
|
|
|
### During import
|
|
|
|
If your CSV contains an ISO8601 date or date-time, you don't need to do anything; readr will automatically recognize it:
|
|
|
|
```{r}
|
|
#| message: false
|
|
csv <- "
|
|
date,datetime
|
|
2022-01-02,2022-01-02 05:12
|
|
"
|
|
read_csv(csv)
|
|
```
|
|
|
|
If you haven't heard of **ISO8601** before, it's an international standard[^datetimes-2] for writing dates where the components of a date are organised from biggest to smallest separated by `-`. For example, in ISO8601 March 5 2022 is `2022-05-03`. ISO8601 dates can also include times, where hour, minute, and second are separated by `:`, and the date and time components are separated by either a `T` or a space.
|
|
For example, you could write 4:26pm on March 5 2022 as either `2022-05-03 16:26` or `2022-05-03T16:26`.
|
|
|
|
[^datetimes-2]: <https://xkcd.com/1179/>
|
|
|
|
For other date-time formats, you'll need to use `col_types` plus `col_date()` or `col_datetime()` along with a date-time format.
|
|
The date-time format used by readr is a standard used across many programming languages, describing a date component with a `%` followed by a single character.
|
|
For example, `%Y-%m-%d` specifies a date that's a year, `-`, month (as number) `-`, day.
|
|
Table @tbl-date-formats lists all the options.
|
|
|
|
| Type | Code | Meaning | Example |
|
|
|-------|-------|--------------------------------|-----------------|
|
|
| Year | `%Y` | 4 digit year | 2021 |
|
|
| | `%y` | 2 digit year | 21 |
|
|
| Month | `%m` | Number | 2 |
|
|
| | `%b` | Abbreviated name | Feb |
|
|
| | `%B` | Full name | Februrary |
|
|
| Day | `%d` | Two digits | 02 |
|
|
| | `%e` | One or two digits | 2 |
|
|
| Time | `%H` | 24-hour hour | 13 |
|
|
| | `%I` | 12-hour hour | 1 |
|
|
| | `%p` | AM/PM | pm |
|
|
| | `%M` | Minutes | 35 |
|
|
| | `%S` | Seconds | 45 |
|
|
| | `%OS` | Seconds with decimal component | 45.35 |
|
|
| | `%Z` | Time zone name | America/Chicago |
|
|
| | `%z` | Offset from UTC | +0800 |
|
|
| Other | `%.` | Skip one non-digit | : |
|
|
| | `%*` | Skip any number of non-digits | |
|
|
|
|
: All date formats understood by readr {#tbl-date-formats}
|
|
|
|
And this code shows some a few options applied to a very ambiguous date:
|
|
|
|
```{r}
|
|
#| messages: false
|
|
|
|
csv <- "
|
|
date
|
|
01/02/15
|
|
"
|
|
|
|
read_csv(csv, col_types = cols(date = col_date("%m/%d/%y")))
|
|
|
|
read_csv(csv, col_types = cols(date = col_date("%d/%m/%y")))
|
|
|
|
read_csv(csv, col_types = cols(date = col_date("%y/%m/%d")))
|
|
```
|
|
|
|
Note that no matter how you specify the date format, it's always displayed the same way once you get it into R.
|
|
|
|
If you're using `%b` or `%B` and working with non-English dates, you'll also need to provide a `locale()`.
|
|
See the list of built-in languages in `date_names_langs()`, or create your own with `date_names()`,
|
|
|
|
### From strings
|
|
|
|
The date-time specification language is powerful, but requires careful analysis of the date format.
|
|
An alternative approach is to use lubridate's helpers which attempt to automatically determine the format once you specify the order of the component.
|
|
To use them, identify the order in which year, month, and day appear in your dates, then arrange "y", "m", and "d" in the same order.
|
|
That gives you the name of the lubridate function that will parse your date.
|
|
For example:
|
|
|
|
```{r}
|
|
ymd("2017-01-31")
|
|
mdy("January 31st, 2017")
|
|
dmy("31-Jan-2017")
|
|
```
|
|
|
|
`ymd()` and friends create dates.
|
|
To create a date-time, add an underscore and one or more of "h", "m", and "s" to the name of the parsing function:
|
|
|
|
```{r}
|
|
ymd_hms("2017-01-31 20:11:59")
|
|
mdy_hm("01/31/2017 08:01")
|
|
```
|
|
|
|
You can also force the creation of a date-time from a date by supplying a timezone:
|
|
|
|
```{r}
|
|
ymd("2017-01-31", tz = "UTC")
|
|
```
|
|
|
|
### From individual components
|
|
|
|
Instead of a single string, sometimes you'll have the individual components of the date-time spread across multiple columns.
|
|
This is what we have in the `flights` data:
|
|
|
|
```{r}
|
|
flights |>
|
|
select(year, month, day, hour, minute)
|
|
```
|
|
|
|
To create a date/time from this sort of input, use `make_date()` for dates, or `make_datetime()` for date-times:
|
|
|
|
```{r}
|
|
flights |>
|
|
select(year, month, day, hour, minute) |>
|
|
mutate(departure = make_datetime(year, month, day, hour, minute))
|
|
```
|
|
|
|
Let's do the same thing for each of the four time columns in `flights`.
|
|
The times are represented in a slightly odd format, so we use modulus arithmetic to pull out the hour and minute components.
|
|
Once we've created the date-time variables, we focus in on the variables we'll explore in the rest of the chapter.
|
|
|
|
```{r}
|
|
make_datetime_100 <- function(year, month, day, time) {
|
|
make_datetime(year, month, day, time %/% 100, time %% 100)
|
|
}
|
|
|
|
flights_dt <- flights |>
|
|
filter(!is.na(dep_time), !is.na(arr_time)) |>
|
|
mutate(
|
|
dep_time = make_datetime_100(year, month, day, dep_time),
|
|
arr_time = make_datetime_100(year, month, day, arr_time),
|
|
sched_dep_time = make_datetime_100(year, month, day, sched_dep_time),
|
|
sched_arr_time = make_datetime_100(year, month, day, sched_arr_time)
|
|
) |>
|
|
select(origin, dest, ends_with("delay"), ends_with("time"))
|
|
|
|
flights_dt
|
|
```
|
|
|
|
With this data, we can visualize the distribution of departure times across the year:
|
|
|
|
```{r}
|
|
#| fig.alt: >
|
|
#| A frequency polyon with departure time (Jan-Dec 2013) on the x-axis
|
|
#| and number of flights on the y-axis (0-1000). The frequency polygon
|
|
#| is binned by day so you see a time series of flights by day. The
|
|
#| pattern is dominated by a weekly pattern; there are fewer flights
|
|
#| on weekends. The are few days that stand out as having a surprisingly
|
|
#| few flights in early Februrary, early July, late November, and late
|
|
#| December.
|
|
flights_dt |>
|
|
ggplot(aes(dep_time)) +
|
|
geom_freqpoly(binwidth = 86400) # 86400 seconds = 1 day
|
|
```
|
|
|
|
Or within a single day:
|
|
|
|
```{r}
|
|
#| fig.alt: >
|
|
#| A frequency polygon with departure time (6am - midnight Jan 1) on the
|
|
#| x-axis, number of flights on the y-axis (0-17), binned into 10 minute
|
|
#| increments. It's hard to see much pattern because of high variability,
|
|
#| but most bins have 8-12 flights, and there are markedly fewer flights
|
|
#| before 6am and after 8pm.
|
|
flights_dt |>
|
|
filter(dep_time < ymd(20130102)) |>
|
|
ggplot(aes(dep_time)) +
|
|
geom_freqpoly(binwidth = 600) # 600 s = 10 minutes
|
|
```
|
|
|
|
Note that when you use date-times in a numeric context (like in a histogram), 1 means 1 second, so a binwidth of 86400 means one day.
|
|
For dates, 1 means 1 day.
|
|
|
|
### From other types
|
|
|
|
You may want to switch between a date-time and a date.
|
|
That's the job of `as_datetime()` and `as_date()`:
|
|
|
|
```{r}
|
|
as_datetime(today())
|
|
as_date(now())
|
|
```
|
|
|
|
Sometimes you'll get date/times as numeric offsets from the "Unix Epoch", 1970-01-01.
|
|
If the offset is in seconds, use `as_datetime()`; if it's in days, use `as_date()`.
|
|
|
|
```{r}
|
|
as_datetime(60 * 60 * 10)
|
|
as_date(365 * 10 + 2)
|
|
```
|
|
|
|
### Exercises
|
|
|
|
1. What happens if you parse a string that contains invalid dates?
|
|
|
|
```{r}
|
|
#| eval: false
|
|
|
|
ymd(c("2010-10-10", "bananas"))
|
|
```
|
|
|
|
2. What does the `tzone` argument to `today()` do?
|
|
Why is it important?
|
|
|
|
3. For each of the following date-times show how you'd parse it using a readr column-specification and a lubridate function.
|
|
|
|
```{r}
|
|
d1 <- "January 1, 2010"
|
|
d2 <- "2015-Mar-07"
|
|
d3 <- "06-Jun-2017"
|
|
d4 <- c("August 19 (2015)", "July 1 (2015)")
|
|
d5 <- "12/30/14" # Dec 30, 2014
|
|
t1 <- "1705"
|
|
t2 <- "11:15:10.12 PM"
|
|
```
|
|
|
|
## Date-time components
|
|
|
|
Now that you know how to get date-time data into R's date-time data structures, let's explore what you can do with them.
|
|
This section will focus on the accessor functions that let you get and set individual components.
|
|
The next section will look at how arithmetic works with date-times.
|
|
|
|
### Getting components
|
|
|
|
You can pull out individual parts of the date with the accessor functions `year()`, `month()`, `mday()` (day of the month), `yday()` (day of the year), `wday()` (day of the week), `hour()`, `minute()`, and `second()`.
|
|
|
|
```{r}
|
|
datetime <- ymd_hms("2026-07-08 12:34:56")
|
|
|
|
year(datetime)
|
|
month(datetime)
|
|
mday(datetime)
|
|
|
|
yday(datetime)
|
|
wday(datetime)
|
|
```
|
|
|
|
For `month()` and `wday()` you can set `label = TRUE` to return the abbreviated name of the month or day of the week.
|
|
Set `abbr = FALSE` to return the full name.
|
|
|
|
```{r}
|
|
month(datetime, label = TRUE)
|
|
wday(datetime, label = TRUE, abbr = FALSE)
|
|
```
|
|
|
|
We can use `wday()` to see that more flights depart during the week than on the weekend:
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A bar chart with days of the week on the x-axis and number of
|
|
#| flights on the y-axis. Monday-Friday have roughly the same number of
|
|
#| flights, ~48,0000, decreasingly slightly over the course of the week.
|
|
#| Sunday is a little lower (~45,000), and Saturday is much lower
|
|
#| (~38,000).
|
|
flights_dt |>
|
|
mutate(wday = wday(dep_time, label = TRUE)) |>
|
|
ggplot(aes(x = wday)) +
|
|
geom_bar()
|
|
```
|
|
|
|
There's an interesting pattern if we look at the average departure delay by minute within the hour.
|
|
It looks like flights leaving in minutes 20-30 and 50-60 have much lower delays than the rest of the hour!
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A line chart with minute of actual departure (0-60) on the x-axis and
|
|
#| average delay (4-20) on the y-axis. Average delay starts at (0, 12),
|
|
#| steadily increases to (18, 20), then sharply drops, hitting at minimum
|
|
#| at ~23 minute past the hour and 9 minutes of delay. It then increases
|
|
#| again to (17, 35), and sharply decreases to (55, 4). It finishes off
|
|
#| with an increase to (60, 9).
|
|
flights_dt |>
|
|
mutate(minute = minute(dep_time)) |>
|
|
group_by(minute) |>
|
|
summarize(
|
|
avg_delay = mean(dep_delay, na.rm = TRUE),
|
|
n = n()) |>
|
|
ggplot(aes(minute, avg_delay)) +
|
|
geom_line()
|
|
```
|
|
|
|
Interestingly, if we look at the *scheduled* departure time we don't see such a strong pattern:
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A line chart with minute of scheduled departure (0-60) on the x-axis
|
|
#| and average delay (4-16). There is relatively little pattern, just a
|
|
#| small suggestion that the average delay decreases from maybe 10 minutes
|
|
#| to 8 minutes over the course of the hour.
|
|
sched_dep <- flights_dt |>
|
|
mutate(minute = minute(sched_dep_time)) |>
|
|
group_by(minute) |>
|
|
summarize(
|
|
avg_delay = mean(arr_delay, na.rm = TRUE),
|
|
n = n())
|
|
|
|
ggplot(sched_dep, aes(minute, avg_delay)) +
|
|
geom_line()
|
|
```
|
|
|
|
So why do we see that pattern with the actual departure times?
|
|
Well, like much data collected by humans, there's a strong bias towards flights leaving at "nice" departure times.
|
|
Always be alert for this sort of pattern whenever you work with data that involves human judgement!
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A line plot with departure minute (0-60) on the x-axis and number of
|
|
#| flights (0-60000) on the y-axis. Most flights are scheduled to depart
|
|
#| on either the hour (~60,000) or the half hour (~35,000). Otherwise,
|
|
#| all most all flights are scheduled to depart on multiples of five,
|
|
#| with a few extra at 15, 45, and 55 minutes.
|
|
ggplot(sched_dep, aes(minute, n)) +
|
|
geom_line()
|
|
```
|
|
|
|
### Rounding
|
|
|
|
An alternative approach to plotting individual components is to round the date to a nearby unit of time, with `floor_date()`, `round_date()`, and `ceiling_date()`.
|
|
Each function takes a vector of dates to adjust and then the name of the unit round down (floor), round up (ceiling), or round to.
|
|
This, for example, allows us to plot the number of flights per week:
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A line plot with week (Jan-Dec 2013) on the x-axis and number of
|
|
#| flights (2,000-7,000) on the y-axis. The pattern is fairly flat from
|
|
#| February to November with around 7,000 flights per week. There are
|
|
#| far fewer flights on the first (approximately 4,500 flights) and last
|
|
#| weeks of the year (approximately 2,500 flights).
|
|
flights_dt |>
|
|
count(week = floor_date(dep_time, "week")) |>
|
|
ggplot(aes(week, n)) +
|
|
geom_line() +
|
|
geom_point()
|
|
```
|
|
|
|
You can use rounding to show the distribution of flights across the course of a day by computing the difference between `dep_time` and the earliest instant of that day:
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A line plot with depature time on the x-axis. This is units of seconds
|
|
#| since midnight so it's hard to interpret.
|
|
flights_dt |>
|
|
mutate(dep_hour = dep_time - floor_date(dep_time, "day")) |>
|
|
ggplot(aes(dep_hour)) +
|
|
geom_freqpoly(binwidth = 60 * 30)
|
|
```
|
|
|
|
Computing the difference between a pair of date-times yields a difftime (more on that in @sec-intervals).
|
|
We can convert that to an `hms` object to get a more useful x-axis:
|
|
|
|
```{r}
|
|
#| fig-alt: >
|
|
#| A line plot with depature time (midnight to midnight) on the x-axis
|
|
#| and number of flights on the y-axis (0 to 15,000). There are very few
|
|
#| (<100) flights before 5am. The number of flights then rises rapidly
|
|
#| to 12,000 / hour, peaking at 15,000 at 9am, before falling to around
|
|
#| 8,000 / hour for 10am to 2pm. Number of flights then increases to
|
|
#| around 12,000 per hour until 8pm, when they rapidly drop again.
|
|
flights_dt |>
|
|
mutate(dep_hour = hms::as_hms(dep_time - floor_date(dep_time, "day"))) |>
|
|
ggplot(aes(dep_hour)) +
|
|
geom_freqpoly(binwidth = 60 * 30)
|
|
```
|
|
|
|
### Modifying components
|
|
|
|
You can also use each accessor function to modify the components of a date/time:
|
|
|
|
```{r}
|
|
(datetime <- ymd_hms("2026-07-08 12:34:56"))
|
|
|
|
year(datetime) <- 2030
|
|
datetime
|
|
month(datetime) <- 01
|
|
datetime
|
|
hour(datetime) <- hour(datetime) + 1
|
|
datetime
|
|
```
|
|
|
|
Alternatively, rather than modifying an existing variabke, you can create a new date-time with `update()`.
|
|
This also allows you to set multiple values in one step:
|
|
|
|
```{r}
|
|
update(datetime, year = 2030, month = 2, mday = 2, hour = 2)
|
|
```
|
|
|
|
If values are too big, they will roll-over:
|
|
|
|
```{r}
|
|
update(ymd("2023-02-01"), mday = 30)
|
|
update(ymd("2023-02-01"), hour = 400)
|
|
```
|
|
|
|
### Exercises
|
|
|
|
1. How does the distribution of flight times within a day change over the course of the year?
|
|
|
|
2. Compare `dep_time`, `sched_dep_time` and `dep_delay`.
|
|
Are they consistent?
|
|
Explain your findings.
|
|
|
|
3. Compare `air_time` with the duration between the departure and arrival.
|
|
Explain your findings.
|
|
(Hint: consider the location of the airport.)
|
|
|
|
4. How does the average delay time change over the course of a day?
|
|
Should you use `dep_time` or `sched_dep_time`?
|
|
Why?
|
|
|
|
5. On what day of the week should you leave if you want to minimise the chance of a delay?
|
|
|
|
6. What makes the distribution of `diamonds$carat` and `flights$sched_dep_time` similar?
|
|
|
|
7. Confirm my hypothesis that the early departures of flights in minutes 20-30 and 50-60 are caused by scheduled flights that leave early.
|
|
Hint: create a binary variable that tells you whether or not a flight was delayed.
|
|
|
|
## Time spans
|
|
|
|
Next you'll learn about how arithmetic with dates works, including subtraction, addition, and division.
|
|
Along the way, you'll learn about three important classes that represent time spans:
|
|
|
|
- **Durations**, which represent an exact number of seconds.
|
|
- **Periods**, which represent human units like weeks and months.
|
|
- **Intervals**, which represent a starting and ending point.
|
|
|
|
How do you pick between duration, periods, and intervals?
|
|
As always, pick the simplest data structure that solves your problem.
|
|
If you only care about physical time, use a duration; if you need to add human times, use a period; if you need to figure out how long a span is in human units, use an interval.
|
|
|
|
### Durations
|
|
|
|
In R, when you subtract two dates, you get a difftime object:
|
|
|
|
```{r}
|
|
# How old is Hadley?
|
|
h_age <- today() - ymd("1979-10-14")
|
|
h_age
|
|
```
|
|
|
|
A difftime class object records a time span of seconds, minutes, hours, days, or weeks.
|
|
This ambiguity can make difftimes a little painful to work with, so lubridate provides an alternative which always uses seconds: the **duration**.
|
|
|
|
```{r}
|
|
as.duration(h_age)
|
|
```
|
|
|
|
Durations come with a bunch of convenient constructors:
|
|
|
|
```{r}
|
|
dseconds(15)
|
|
dminutes(10)
|
|
dhours(c(12, 24))
|
|
ddays(0:5)
|
|
dweeks(3)
|
|
dyears(1)
|
|
```
|
|
|
|
Durations always record the time span in seconds.
|
|
Larger units are created by converting minutes, hours, days, weeks, and years to seconds: 60 seconds in a minute, 60 minutes in an hour, 24 hours in a day, and 7 days in a week.
|
|
Larger time units are more problematic.
|
|
A year is uses the "average" number of days in a year, i.e. 365.25.
|
|
There's no way to convert a month to a duration, because there's just too much variation.
|
|
|
|
You can add and multiply durations:
|
|
|
|
```{r}
|
|
2 * dyears(1)
|
|
dyears(1) + dweeks(12) + dhours(15)
|
|
```
|
|
|
|
You can add and subtract durations to and from days:
|
|
|
|
```{r}
|
|
tomorrow <- today() + ddays(1)
|
|
last_year <- today() - dyears(1)
|
|
```
|
|
|
|
However, because durations represent an exact number of seconds, sometimes you might get an unexpected result:
|
|
|
|
```{r}
|
|
one_pm <- ymd_hms("2026-03-12 13:00:00", tz = "America/New_York")
|
|
|
|
one_pm
|
|
one_pm + ddays(1)
|
|
```
|
|
|
|
Why is one day after 1pm March 12, 2pm March 13?
|
|
If you look carefully at the date you might also notice that the time zones have changed.
|
|
March 12 only has 23 hours because it's when DST starts, so if we add a full days worth of seconds we end up with a different time.
|
|
|
|
### Periods
|
|
|
|
To solve this problem, lubridate provides **periods**.
|
|
Periods are time spans but don't have a fixed length in seconds, instead they work with "human" times, like days and months.
|
|
That allows them to work in a more intuitive way:
|
|
|
|
```{r}
|
|
one_pm
|
|
one_pm + days(1)
|
|
```
|
|
|
|
Like durations, periods can be created with a number of friendly constructor functions.
|
|
|
|
```{r}
|
|
hours(c(12, 24))
|
|
days(7)
|
|
months(1:6)
|
|
```
|
|
|
|
You can add and multiply periods:
|
|
|
|
```{r}
|
|
10 * (months(6) + days(1))
|
|
days(50) + hours(25) + minutes(2)
|
|
```
|
|
|
|
And of course, add them to dates.
|
|
Compared to durations, periods are more likely to do what you expect:
|
|
|
|
```{r}
|
|
# A leap year
|
|
ymd("2024-01-01") + dyears(1)
|
|
ymd("2024-01-01") + years(1)
|
|
|
|
# Daylight Savings Time
|
|
one_pm + ddays(1)
|
|
one_pm + days(1)
|
|
```
|
|
|
|
Let's use periods to fix an oddity related to our flight dates.
|
|
Some planes appear to have arrived at their destination *before* they departed from New York City.
|
|
|
|
```{r}
|
|
flights_dt |>
|
|
filter(arr_time < dep_time)
|
|
```
|
|
|
|
These are overnight flights.
|
|
We used the same date information for both the departure and the arrival times, but these flights arrived on the following day.
|
|
We can fix this by adding `days(1)` to the arrival time of each overnight flight.
|
|
|
|
```{r}
|
|
flights_dt <- flights_dt |>
|
|
mutate(
|
|
overnight = arr_time < dep_time,
|
|
arr_time = arr_time + days(if_else(overnight, 0, 1)),
|
|
sched_arr_time = sched_arr_time + days(overnight * 1)
|
|
)
|
|
```
|
|
|
|
Now all of our flights obey the laws of physics.
|
|
|
|
```{r}
|
|
flights_dt |>
|
|
filter(overnight, arr_time < dep_time)
|
|
```
|
|
|
|
### Intervals {#sec-intervals}
|
|
|
|
It's obvious what `dyears(1) / ddays(365)` should return: one, because durations are always represented by a number of seconds, and a duration of a year is defined as 365 days worth of seconds.
|
|
|
|
What should `years(1) / days(1)` return?
|
|
Well, if the year was 2015 it should return 365, but if it was 2016, it should return 366!
|
|
There's not quite enough information for lubridate to give a single clear answer.
|
|
What it does instead is give an estimate:
|
|
|
|
```{r}
|
|
years(1) / days(1)
|
|
```
|
|
|
|
If you want a more accurate measurement, you'll have to use an **interval**.
|
|
An interval is a pair of starting and ending date times, or you can think of it as a duration with a starting point.
|
|
|
|
You can create an interval by writing `start %--% end`:
|
|
|
|
```{r}
|
|
y2023 <- ymd("2023-01-01") %--% ymd("2024-01-01")
|
|
y2024 <- ymd("2024-01-01") %--% ymd("2025-01-01")
|
|
|
|
y2023
|
|
y2024
|
|
```
|
|
|
|
You could then divide it by `days()` to find out how many days fit in the year:
|
|
|
|
```{r}
|
|
y2023 / days(1)
|
|
y2024 / days(1)
|
|
```
|
|
|
|
### Exercises
|
|
|
|
1. Explain `days(overnight * 1)` to someone who has just started learning R.
|
|
How does it work?
|
|
|
|
2. Create a vector of dates giving the first day of every month in 2015.
|
|
Create a vector of dates giving the first day of every month in the *current* year.
|
|
|
|
3. Write a function that given your birthday (as a date), returns how old you are in years.
|
|
|
|
4. Why can't `(today() %--% (today() + years(1))) / months(1)` work?
|
|
|
|
## Time zones
|
|
|
|
Time zones are an enormously complicated topic because of their interaction with geopolitical entities.
|
|
Fortunately we don't need to dig into all the details as they're not all important for data analysis, but there are a few challenges we'll need to tackle head on.
|
|
|
|
<!--# https://www.ietf.org/timezones/tzdb-2018a/theory.html -->
|
|
|
|
The first challenge is that everyday names of time zones tend to be ambiguous.
|
|
For example, if you're American you're probably familiar with EST, or Eastern Standard Time.
|
|
However, both Australia and Canada also have EST!
|
|
To avoid confusion, R uses the international standard IANA time zones.
|
|
These use a consistent naming scheme `{area}/{location}`, typically in the form `{continent}/{city}` or `{ocean}/{city}`.
|
|
Examples include "America/New_York", "Europe/Paris", and "Pacific/Auckland".
|
|
|
|
You might wonder why the time zone uses a city, when typically you think of time zones as associated with a country or region within a country.
|
|
This is because the IANA database has to record decades worth of time zone rules.
|
|
Over the course of decades, countries change names (or break apart) fairly frequently, but city names tend to stay the same.
|
|
Another problem is that the name needs to reflect not only the current behavior, but also the complete history.
|
|
For example, there are time zones for both "America/New_York" and "America/Detroit".
|
|
These cities both currently use Eastern Standard Time but in 1969-1972 Michigan (the state in which Detroit is located), did not follow DST, so it needs a different name.
|
|
It's worth reading the raw time zone database (available at <https://www.iana.org/time-zones>) just to read some of these stories!
|
|
|
|
You can find out what R thinks your current time zone is with `Sys.timezone()`:
|
|
|
|
```{r}
|
|
Sys.timezone()
|
|
```
|
|
|
|
(If R doesn't know, you'll get an `NA`.)
|
|
|
|
And see the complete list of all time zone names with `OlsonNames()`:
|
|
|
|
```{r}
|
|
length(OlsonNames())
|
|
head(OlsonNames())
|
|
```
|
|
|
|
In R, the time zone is an attribute of the date-time that only controls printing.
|
|
For example, these three objects represent the same instant in time:
|
|
|
|
```{r}
|
|
x1 <- ymd_hms("2024-06-01 12:00:00", tz = "America/New_York")
|
|
x1
|
|
|
|
x2 <- ymd_hms("2024-06-01 18:00:00", tz = "Europe/Copenhagen")
|
|
x2
|
|
|
|
x3 <- ymd_hms("2024-06-02 04:00:00", tz = "Pacific/Auckland")
|
|
x3
|
|
```
|
|
|
|
You can verify that they're the same time using subtraction:
|
|
|
|
```{r}
|
|
x1 - x2
|
|
x1 - x3
|
|
```
|
|
|
|
Unless otherwise specified, lubridate always uses UTC.
|
|
UTC (Coordinated Universal Time) is the standard time zone used by the scientific community and is roughly equivalent to GMT (Greenwich Mean Time).
|
|
It does not have DST, which makes a convenient representation for computation.
|
|
Operations that combine date-times, like `c()`, will often drop the time zone.
|
|
In that case, the date-times will display in your local time zone:
|
|
|
|
```{r}
|
|
x4 <- c(x1, x2, x3)
|
|
x4
|
|
```
|
|
|
|
You can change the time zone in two ways:
|
|
|
|
- Keep the instant in time the same, and change how it's displayed.
|
|
Use this when the instant is correct, but you want a more natural display.
|
|
|
|
```{r}
|
|
x4a <- with_tz(x4, tzone = "Australia/Lord_Howe")
|
|
x4a
|
|
x4a - x4
|
|
```
|
|
|
|
(This also illustrates another challenge of times zones: they're not all integer hour offsets!)
|
|
|
|
- Change the underlying instant in time.
|
|
Use this when you have an instant that has been labelled with the incorrect time zone, and you need to fix it.
|
|
|
|
```{r}
|
|
x4b <- force_tz(x4, tzone = "Australia/Lord_Howe")
|
|
x4b
|
|
x4b - x4
|
|
```
|
|
|
|
## Summary
|
|
|
|
This chapter has introduced you to the tools that lubridate provides to help you work with date-time data.
|
|
Working with dates and times can seem harder than necessary, but hopefully this chapter has helped you see why --- date-times are more complex than they seem at first glance, and handling every possible situation adds complexity.
|
|
Even if your data never crosses a day light savings boundary or involves a leap year, the functions need to be able to handle it.
|
|
|
|
The next chapter gives a round up of missing values.
|
|
You've seen them in a few places and have no doubt encounter in your own analysis, and it's how time to provide a grab bag of useful techniques for dealing with them.
|