r - timeBasedSeq 函数会重复它创建的某些日子！

Question

我正在使用 xts 包中的 timeBasedSeq 函数用作时间序列/动物园对象中的索引，但它会重复它创建的某些日子！这会导致动物园出现问题，因为“'order.by'中的索引条目不是唯一的”。

例如：

timeBasedSeq("19860601/19861231")

将创建

...“1986-10-23”“1986-10-24”“1986-10-25”“1986-10-26”“1986-10-26”“1986-10-27”“1986-10- 28" "1986-10-29""...

所以由于某种原因它重复了 26 日，我该如何避免呢？

Answer 1

该错误在 xts 中。该函数使用seq.POSIXct，并且可以通过以下方式产生相同的行为：

seq(as.POSIXct("1986-10-01"), as.POSIXct("1986-11-01"), by="day")

更令我惊讶的是seq.POSIXlt

seq(as.POSIXlt("1986-10-01"), as.POSIXlt("1986-11-01"), by="day")

但是这种行为在 seq.POSIXt 中有很好的记录，并且有一个规定，by="DSTday"xts 作者可能应该在天是隐含间隔的情况下使用它。临时解决方法是：

timeBasedSeq("19860601/19861231")[ !duplicated(timeBasedSeq("19860601/19861231") ]

或更简洁：

unique(timeBasedSeq("19860601/19861231"))

Answer 2

扩展我的评论：这个错误是由夏令时的变化来解释的。此时添加 24 小时不会将您移动到第二天，因为已经取消了一个小时。

> as.POSIXct("1986-10-26")
[1] "1986-10-26 BST"
> as.POSIXct("1986-10-26")+24*60*60
[1] "1986-10-26 23:00:00 GMT"

一种解决方案是将日期设置在中午而不是午夜，并用于strftime去除无关信息并as.POSIXlt转换回时间类：

as.POSIXlt(strftime(timeBasedSeq("19860601 1200/19861231/d"),"%Y-%m-%d"))

...
[146] "1986-10-24" "1986-10-25" "1986-10-26" "1986-10-27" "1986-10-28"
...