我有一个气候时间序列,缺少日期(不缺少值)。例如:
n = 15
full.dates = seq(Sys.Date(), by = 'day', length = n)
serie.dates = full.dates[c(1:10, 12, 15)] # missing 11, 13, 14
y = rnorm(n)
require(zoo)
serie = zoo(y, serie.dates)
如果给出'full.dates'向量,我如何'填充'(使用插值)这些缺失点? 谢谢!
答案 0 :(得分:18)
合并"空"具有所需日期的对象,然后使用na.approx
(或na.spline
等)填写缺失的值。
x <- merge(serie, zoo(,seq(start(serie),end(serie),by="day")), all=TRUE)
x <- na.approx(x)