r 中 xts 上的最后一次观察结转（locf）？ （不像na.approx？）

Question

我希望将“最后一次观察结转 (LOCF)”从一个 xts 时间序列转移到另一个包含我想使用的时间戳的 xts 系列。我尝试过使用na.locf，但xout 似乎不起作用，正如我从na.approx 所知道的那样。

任何建议将不胜感激！

这是一个例子：

我有两个 xts 时间序列数据集。一个带有我的数据，另一个带有我想使用的时间戳。

xts_with_data <- as.xts(read.zoo(text='
2016-07-01 00:00:10,   1.0
2016-07-01 00:00:14,   2.0
2016-07-01 00:00:20,   3.0
2016-07-01 00:00:29,   4.0
2016-07-01 00:00:34,   5.0
2016-07-01 00:00:39,   6.0
', sep=',', index=1, tz='', format="%Y-%m-%d %H:%M:%S"))
names(xts_with_data) <- c('x')

xts_with_timestamps <- as.xts(read.zoo(text='
2016-07-01 00:00:15,   0.0
2016-07-01 00:00:20,   0.0
2016-07-01 00:00:30,   0.0
2016-07-01 00:00:35,   0.0
2016-07-01 00:00:38,   0.0
', sep=',', index=1, tz='', format="%Y-%m-%d %H:%M:%S"))

我想要的是这个：

                    [,1]
2016-07-01 00:00:15    2
2016-07-01 00:00:20    3
2016-07-01 00:00:30    4
2016-07-01 00:00:35    5
2016-07-01 00:00:38    5

（即来自xts_with_timestamps 的时间戳和来自xts_with_data 的相应locf）。

我想我可以像这样使用xout 来做到这一点（这对na.approx 很有效）：

na.locf(xts_with_data, xout = index(xts_with_timestamps))

但这只是返回我原来的xts_with_data。

有什么建议吗？

提前致谢。

Answer 1

在申请na.locf()之前使用merge()加入这两个系列可以是一种方法。

xts_all <- merge(xts_with_data, xts_with_timestamps)[,-2]
na.locf(xts_all)[index(xts_with_timestamps)]

#                     x
# 2016-07-01 00:00:15 2
# 2016-07-01 00:00:20 3
# 2016-07-01 00:00:30 4
# 2016-07-01 00:00:35 5
# 2016-07-01 00:00:38 5

似乎xout 适用于na.approx() 但不适用于na.locf() 的原因似乎是后者有一个xts 方法似乎不能正确传递额外的参数。但是，您可以强制使用默认方法，从而恢复对 xout 的支持

na.locf.default(xts_with_data, xout=index(xts_with_timestamps))
# or
na.locf(zoo(xts_with_data), xout=index(xts_with_timestamps))

#                     x
# 2016-07-01 00:00:15 2
# 2016-07-01 00:00:20 3
# 2016-07-01 00:00:30 4
# 2016-07-01 00:00:35 5
# 2016-07-01 00:00:38 5