如何通过考虑工作日的日期序列对xts进行子集化

时间:2016-08-21 06:48:03

标签: r date subset xts

我使用as.Date序列对xts对象进行了子集化:

library("quantmod")

getSymbols("AAPL",from="2003-01-01")
# indx sequence
indx <- seq(as.Date('2003-03-31'), length.out=200, by='4 weeks')

SELECT <- AAPL[paste(indx)]

如果您查看SELECT中的最后两行,我会看到它从2016-06-06跳到2016-08-01,这不是4周。它缺少2016-07-04,但由于这不是一个工作日,它会跳过它。如果SELECT不是工作日,如何返回indx,这将返回下一个工作日?在此示例中,它应返回2016-07-05 ...

1 个答案:

答案 0 :(得分:1)

在工作日,我猜您的意思是AAPL的交易日,在这种情况下,您的工作日实际上是AAPL安全的时间索引。

使用事实日期的第一原则方法可以增加1:

indx <- seq(as.Date('2003-03-31'), length.out=200, by='4 weeks')
indx <- indx[indx < as.Date(end(AAPL))]
while(!all(indx %in% as.Date(index(AAPL)))) {
    # You ask for the next available business day:
    indx[!indx %in% as.Date(index(AAPL))] <- indx[!indx %in% as.Date(index(AAPL))] + 1
    # Careful that the last indx value does not go past as.Date(end(AAPL))
    if (indx[length(indx)] > as.Date(end(AAPL))) {
        indx[length(indx)] <- as.Date(end(AAPL))
    }
}

SELECT <- AAPL[indx]
tail(SELECT)
# AAPL.Open AAPL.High AAPL.Low AAPL.Close AAPL.Volume AAPL.Adjusted
# 2016-03-14    101.91    102.91   101.78     102.52    25076100     101.35055
# 2016-04-11    108.97    110.61   108.83     109.02    29407500     107.77640
# 2016-05-09     93.00     93.77    92.59      92.79    32936400      92.29005
# 2016-06-06     97.99    101.89    97.55      98.63    23292500      98.09858
# 2016-07-05     95.39     95.40    94.46      94.99    27705200      94.47819
# 2016-08-01    104.41    106.15   104.41     106.05    38167900     105.47860

您还可以通过timeDate包找到解决方案,更常用于商业日期类型子集。例如http://stackoverflow.com/questions/13673895/r-time-series-data-daily-only-working-days