将不规则时间序列转换为规则时间序列

Question

我在将不规则时间序列转换为规则时间序列时遇到问题。下面是一个简化的示例：

require(zoo)
t <- as.character(c(1981,1984,1985))
d <- c(1,3,6)
dt <- data.frame(d,t)
t <- as.Date(t,"%Y")
z <- zoo(d,t)
plot(z)
ts.d <- as.ts(as.zooreg(z,freq=1)) # create a regular ts object
ts.d # regular time series

我想创建一个常规时间序列 ts.d，如下所示 c(1981,NA,NA,1984,1985)。

令人惊奇的是，我第一次运行这个：它有效！但是当我想再次运行它或重复它（as.ts()行）时，它停止工作并且我获得一个很长的时间序列：

ts.d # regular time series
Time Series:
Start = 4299 
End = 5760 
Frequency = 1 
  [1]  1 NA NA NA NA NA NA NA NA NA NA NA NA NA
 [15] NA NA NA NA NA NA NA NA 

等等

。出了什么问题？

Answer 1

正如已经指出的那样， as.Date(as.character(t), "%Y") 是不正确的，因为它没有给出所需的月份和日期。如果我们想将年份转换为 "Date" 类，我们可以使用动物园的 as.yearmon 来执行此 as.Date(as.yearmon(t)) 操作代码>;然而，这样我们就会遇到进一步的问题，不同的年份有不同的天数，因此无法使用日期来表示年份的常规序列。

事实上，我们一开始就不想要约会。我们只想使用年份，在这种情况下，它可以简化为：

> z <- zoo(c(1, 3, 6), c(1981, 1984, 1985))
> 
> as.ts(z)
Time Series:
Start = 1981 
End = 1985 
Frequency = 1 
[1]  1 NA NA  3  6

或者如果我们想要安全，我们可以这样做，这将强制它是每年的，即使输入偶然具有较低的频率：频率(z)＜-1； as.ts(z) 或者只是将原始动物园系列从一开始就定义为频率为 1：

> z <- zoo(c(1, 3, 6), c(1981, 1984, 1985), frequency = 1)
> as.ts(z)
Time Series:
Start = 1981 
End = 1985 
Frequency = 1 
[1]  1 NA NA  3  6

在这个例子中，它没有什么区别，但在这种情况下 z <- Zoo( c(1, 3, 6), c(1981, 1983, 1985), 频率 = 1) 需要显式的 频率 来防止其频率为 0.5 。

Answer 2

这不是一个错误。您的时间序列中有 4 年跨越 1,461 天。当我第一次运行它时，它对我不起作用。 as.Date(t,"%Y") 不知道使用哪个月/日来创建日期，因此它使用今天的月/日。这不利于可重复的分析。试试这个：

t <- c(1981,1984,1985)
d <- c(1,3,6)
z <- zoo(d,t)
z <- merge(z,zoo(,c(1981,1982,1983,1984,1985)))
ts.d <- as.ts(z)

产生：

> ts.d
Time Series:
Start = 1981 
End = 1985 
Frequency = 1 
[1]  1 NA NA  3  6