[R] Aggregate records to 10min

Fankhauser GEP Data Consulting rolf.fankhauser at gepdata.ch
Sun Jan 17 16:17:08 CET 2016


Hi Jim,
Thanks a lot! It works now. I didn't remember how to access the 
datetimes in w10min. names(...) is the solution!

Rolf

Jim Lemon wrote:
> Hi Rolf,
> If I get the above, perhaps if you change the names of w10min after 
> applying the calculation:
>
> raindata<-data.frame(value=round(runif(60,0,4),1),
> ptime=paste("2016-01-17 ","15:",0:59,sep=""))
> t10min <- 600*floor(as.integer(as.POSIXct(raindata$ptime))/600)
> w10min <- tapply(raindata$value,t10min,sum)
> names(w10min)<-format(as.POSIXct(as.numeric(names(w10min)),
>  tz="AEST",origin="1970-01-01"),"%m/%d/%Y %H:%M")
>
> Jim
>
>
> On Sun, Jan 17, 2016 at 5:45 AM, Rolf Fankhauser 
> <rolf.fankhauser at gepdata.ch <mailto:rolf.fankhauser at gepdata.ch>> wrote:
>
>     Hi
>
>     I would like to aggregate a rainfall series with 1min records
>     (timestamp and value of 0.1mm from a tipping bucket raingauge) to
>     10min values by summing up the values.
>
>     # ptime is a POSIXlt datetime value with tz="GMT"
>
>     t10min <- 600*floor(as.integer(as.POSIXct(data$ptime))/600)
>     w10min <- tapply(data$value, format(as.POSIXct(t10min, tz="GMT",
>     origin = "1970-01-01"), "%Y-%m-%d %H:%M"), sum)
>     write.table(as.matrix(w10min),"data 10min.txt", row.names=TRUE,
>     col.names=FALSE, quote=FALSE)
>
>     This code works but I would like to have the result in datetime
>     format of %m/%d/%Y %H:%M. When I output this format the records
>     are not chronologically sorted but text-sorted because dimnames of
>     w10min is of type character (because of the format function).
>     Is there an easier way summing up the records to 10min records?
>
>     Thanks,
>     Rolf
>
>     ______________________________________________
>     R-help at r-project.org <mailto:R-help at r-project.org> mailing list --
>     To UNSUBSCRIBE and more, see
>     https://stat.ethz.ch/mailman/listinfo/r-help
>     PLEASE do read the posting guide
>     http://www.R-project.org/posting-guide.html
>     and provide commented, minimal, self-contained, reproducible code.
>
>


-- 
______________________________

Fankhauser GEP Data Consulting
Hegenheimerstrasse 129
4055 Basel

Tel:    ++41-(0)61-321-4525
Mobile: ++41-(0)79-440-7706
rolf.fankhauser at gepdata.ch
www.gepdata.ch



More information about the R-help mailing list