Создание регулярных 15-минутных временных рядов из нерегулярных временных рядов
У меня есть неправильный временной ряд (с DateTime и RainfallValue) в CSV-файлеC:\SampleData.csv
:
<code> DateTime,RainInches 1/6/2000 11:59,0 1/6/2000 23:59,0.01 1/7/2000 11:59,0 1/13/2000 23:59,0 1/14/2000 0:00,0 1/14/2000 23:59,0 4/14/2000 3:07,0.01 4/14/2000 3:12,0.03 4/14/2000 3:19,0.01 12/31/2001 22:44,0 12/31/2001 22:59,0.07 12/31/2001 23:14,0 12/31/2001 23:29,0 12/31/2001 23:44,0.01 12/31/2001 23:59,0.01 </code>
Note: Нерегулярные временные шаги могут составлять 1 минуту, 15 минут, 1 час и т. Д. Кроме того, может быть несколько наблюдений в желаемом 15-минутном интервале.
Я пытаюсь создать регулярные 15-минутные временные ряды с 2000-01-01 по 2001-12-31, которые должны выглядеть следующим образом:
<code> 2000-01-01 00:15:00 0.00 2000-01-01 00:30:00 0.00 2000-01-01 00:45:00 0.00 ... 2001-12-31 23:30:00 0.01 2001-12-31 23:45:00 0.01 </code>
Note: Временные ряды регулярны с 15-минутными интервалами, заполняя недостающие данные 0. Если в 15-минутных интервалах имеется более одной точки данных, они суммируются.
Вот мой код:
<code> library(zoo) library(xts) filename = "C:\\SampleData.csv" ReadData <- read.zoo(filename, format = "%m/%d/%Y %H:%M", sep=",", tz="UTC", header=TRUE) # read .csv as a ZOO object RawData <- aggregate(ReadData, index(ReadData), sum) # Merge duplicate time stamps and SUM the corresponding data (CAUTION) RawDataSeries <- as.xts(RawData,order.by =index(RawData)) #convert to an XTS object RegularTimes <- seq(as.POSIXct("2000-01-01 00:00:00", tz = "UTC"), as.POSIXct("2001-12-31 23:45:00", tz = "UTC"), by = 60*15) BlankTimeSeries <- xts((rep(0,length(RegularTimes))),order.by = RegularTimes) MergedTimeSeries <- merge(RawDataSeries,BlankTimeSeries) TS_sum15min <- period.apply(MergedTimeSeries,endpoints(MergedTimeSeries, "minutes", 15), sum, na.rm = TRUE ) TS_align15min <- align.time( TS_sum15min [endpoints(TS_sum15min , "minutes", 15)], n=60*15) </code>
Problem: Выходной временной рядTS_align15min
:
(a) has repeating blocks of time-stamps
(b) starts (mysteriously) from 1999, as:
<code>1999-12-31 19:15:00 0 1999-12-31 19:30:00 0 1999-12-31 19:45:00 0 1999-12-31 20:00:00 0 1999-12-31 20:15:00 0 1999-12-31 20:30:00 0 </code>
Что я делаю неправильно?
Спасибо за любое направление!