【问题标题】:How can i obtain a tsibble from this tibble without using the as.Date function?我如何在不使用 as.Date 函数的情况下从此 tibble 获取 tsibble?
【发布时间】:2021-12-12 00:08:27
【问题描述】:

我需要将几个“tibble”转换成“tsibble”。

这里是一个简单的例子:

require(tidyverse)
require(lubridate)

time_1 <- c(ymd_hms('20210101 000000'),
            ymd_hms('20210101 080000'),
            ymd_hms('20210101 160000'),
            # ymd_hms('20210102 000000'),
            ymd_hms('20210102 080000'),
            ymd_hms('20210102 160000'))

df_1 <- tibble(time_1, y=rnorm(5))

df_1 %>% 
  as_tsibble(index=time_1)

这段代码按预期工作。 但是,如果日期都是午夜,则此代码会引发错误:

time_2 <- c(ymd_hms('20210101 000000'),
            ymd_hms('20210102 000000'),
            ymd_hms('20210103 000000'),
            # ymd_hms('20210104 000000'),
            ymd_hms('20210105 000000'),
            ymd_hms('20210106 000000'))

df_2 <- tibble(time_2, y=rnorm(5))

df_2 %>% 
  as_tsibble(index=time_2)

我不想用这种方式解决这个问题,因为as.Date 函数改变了列类型。

df_2 %>% 
  mutate(time_2=as.Date(time_2)) %>% 
  as_tsibble(index=time_2)

我也不想以这种方式解决问题,因为在将 tibble 转换为 tsibble 后,我需要应用 fill_gaps 函数,该函数不会在第二种情况下创建 ymd_hms('20210104 000000')

df_2 %>% 
  as_tsibble(index=time_2, regular=FALSE)

这是一个错误吗?

谢谢。

【问题讨论】:

标签: r tidyverse lubridate tsibble


【解决方案1】:

tsibble 的FAQ 中解释了这种行为。

每天午夜测量的基本上是次日数据 (ymd_hms()) 不一定有 1 天(24 小时)的间隔。考虑到由于您所在时区的夏令时,有些日子会发生变化,因此午夜和第二天午夜之间的小时数可能是 23 或 25 小时。

如果您使用的是每天测量的数据,则应使用具有ymd() 精度的日期。如果愿意,您可以使用 as_datetime() 将其转换回日期时间。

我个人认为这不应该产生错误,但是如果它确实会简单得多。也许这里的适当间隔是 1 小时或 30 分钟(或任何适合指定时区中时区转换的时间)。

【讨论】:

  • 感谢您的宝贵意见。当您说“基本上每天午夜测量的次日数据 (ymd_hms()) 不一定有 1 天(24 小时)的间隔。”我完全同意你的看法。但是这个问题也会影响每 8 小时测量一次的时间序列,不是吗?
  • 是的,没错。以 8 小时为间隔进行夏令时更改的次日数据并不总是有 8 小时的间隔,因此精确的间隔更像是 1 小时。如果您希望有 8 小时的间隔并忽略夏令时更改,您可以使用没有夏令时的时区。
  • 在我报告的示例中,我正确地得到了 8 小时的间隔 (df_1),这就是我所期望的。让我产生疑问的问题是,在午夜的时间间隔内,代码会出错。对于间隔为 8 小时的数据集,它失败的原因也应该是失败的。
  • 是的。您在午夜间隔中得到的错误不一定是必需的,但为了防止人们在不考虑时区时犯下的常见错误而提出。也许最好在这里给出警告。您可以在此处获得所需的 8 小时间隔,因为您使用的数据范围不包含由于时区而导致的时间变化。
猜你喜欢
  • 2014-11-06
  • 2019-06-15
  • 2017-05-10
  • 1970-01-01
  • 2019-10-23
  • 2021-09-17
  • 2019-08-31
  • 1970-01-01
  • 1970-01-01
相关资源
最近更新 更多