我有一个时间数据的数据库。我想将数据插值为特定的时间步长。
Id Time humid humtemp prtemp press t
1 2012-01-21 18:41:50 47.7 14.12 13.870 1005.70 -0.05277778
1 2012-01-21 18:46:43 44.5 15.37 15.100 1005.20 0.02861111
1 2012-01-21 18:51:35 43.2 15.88 15.576 1005.10 0.10972222
1 2012-01-21 18:56:28 42.5 16.17 15.833 1004.90 0.19111111
1 2012-01-21 19:01:21 42.2 16.31 15.986 1004.80 0.27250000
1 2012-01-21 19:06:14 41.8 16.47 16.118 1004.60 0.35388889
1 2012-01-21 19:11:07 41.6 16.51 16.177 1004.60 0.43527778
我希望得到以下时间步长的数据并进行插值。
Id Time humid humtemp prtemp press t
1 2012-01-21 18:45:00 .... ... ..... .... ....
1 2012-01-21 18:50:00 ....
1 2012-01-21 18:55:00 ....
1 2012-01-21 19:00:00 ....
1 2012-01-21 19:05:00 ....
1 2012-01-21 19:10:00 ....
我尝试了不同的方法,但没有找到解决方案。例如,我创建了一个动物园对象。
z <- zoo(MTS01m,order.by=MTS01m$Time)
tstart2<-asP("2012-01-21 18:45:00")
Ts<-1*60
y <- merge(z, zoo(order.by=seq(tstart2, end(z), by=Ts)))
xa <- na.approx(y)
xs <- na.spline(y)
但是出现了错误:
Errore in approx(x[!na], y[!na], xout, ...) :
need at least two non-NA values to interpolate
Inoltre: Warning message:
In xy.coords(x, y) : si è prodotto un NA per coercizione
我创建了一个名为t的辅助索引,该索引从我想要的数据开始,但我不知道如何使用这个索引。你有什么建议吗?
asP
是哪个包中的? - GSee用于将字符串或自1970-01-01以来的秒数轻松转换为POSIXct
switch(class(timeVal [1])[1], character = as.POSIXct(timeVal,tz = tz,...), POSIXct = timeVal, POSIXlt = timeVal, ISOdate(1970,1,1,0)+ timeVal ) } - Marco Giuliani