我有一个气候时间序列,缺少日期(没有值)。例如:
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”向量,我如何“填充”(使用插值法)这些缺失点?
谢谢!
最佳答案
与具有所有所需日期的“空”对象合并,然后使用na.approx
(或na.spline
等)填充缺少的值。
x <- merge(serie, zoo(,seq(start(serie),end(serie),by="day")), all=TRUE)
x <- na.approx(x)