[R] adding date/time stamp

Gabor Grothendieck ggrothendieck at gmail.com
Wed Dec 31 00:51:21 CET 2008


Try this:

> aggregate(z$mph, trunc(time(z), "hour"), mean)
(09/01/05 00:00:00) (09/01/05 01:00:00) (09/01/05 02:00:00)
            9.27500            10.08333             9.20000

On Tue, Dec 30, 2008 at 6:30 PM, Sherri Heck <sheck at ucar.edu> wrote:
> Dear All-
>
> I have a dataset that is comprised of the following (LST = yymmddhhMM):
>
>
>
>    LST   in     mph    Deg   DegF  DegF2    %    volts   Deg    mph2   w/m2
> 0509010000   0.00    7.8  216.9   45.1   -999   24.4   -999   -999   10.6
>  0.2
> 0509010005   0.00    8.6  206.6   45.1   -999   25.2   -999   -999   11.7
>  0.2
> 0509010010   0.00    7.8  199.2   44.9   -999   25.4   -999   -999   12.8
>  0.2
> 0509010015   0.00    7.7  197.4   44.8   -999   25.4   -999   -999   10.4
>  0.2
> 0509010020   0.00    7.6  203.9   44.8   -999   25.3   -999   -999   10.0
>  0.2
> 0509010025   0.00    9.3  200.9   44.9   -999   25.3   -999   -999   11.8
>  0.2
> 0509010030   0.00    9.4  200.3   44.7   -999   25.5   -999   -999   12.2
>  0.2
> 0509010035   0.00   10.0  199.2   44.6   -999   25.9   -999   -999   13.0
>  0.2
> 0509010040   0.00    9.5  201.5   44.5   -999   25.9   -999   -999   13.3
>  0.2
> 0509010045   0.00   10.8  200.4   44.5   -999   26.1   -999   -999   13.0
>  0.2
> 0509010050   0.00   11.8  198.4   44.5   -999   26.1   -999   -999   13.3
>  0.2
> 0509010055   0.00   11.0  197.4   44.5   -999   25.5   -999   -999   13.3
>  0.2
> 0509010100   0.00    9.7  202.0   44.6   -999   25.1   -999   -999   13.0
>  0.2
> 0509010105   0.00    9.0  215.1   44.7   -999   24.9   -999   -999   12.2
>  0.2
> 0509010110   0.00   10.1  223.1   44.6   -999   25.1   -999   -999   13.2
>  0.2
> 0509010115   0.00   10.4  231.2   44.5   -999   25.5   -999   -999   12.0
>  0.2
> 0509010120   0.00   11.0  237.4   44.2   -999   25.9   -999   -999   11.7
>  0.2
> 0509010125   0.00   10.6  241.0   44.2   -999   26.0   -999   -999   11.8
>  0.2
> 0509010130   0.00   11.1  242.2   44.1   -999   26.2   -999   -999   12.2
>  0.2
> 0509010135   0.00   10.6  240.0   44.0   -999   26.5   -999   -999   11.5
>  0.2
> 0509010140   0.00   10.1  241.0   44.0   -999   26.4   -999   -999   11.5
>  0.2
> 0509010145   0.00    9.8  243.2   44.0   -999   26.6   -999   -999   10.7
>  0.2
> 0509010150   0.00    9.3  240.3   43.9   -999   27.0   -999   -999   10.0
>  0.2
> 0509010155   0.00    9.3  239.2   43.8   -999   26.8   -999   -999   10.0
>  0.2
> 0509010200   0.00    9.2  240.1   43.8   -999   26.6   -999   -999    9.8
>  0.2
> 0509010205   0.00    9.0  240.0   43.8   -999   26.6   -999   -999    9.4
>  0.2
> 0509010210   0.00    9.2  245.0   43.9   -999   26.3   -999   -999    9.8
>  0.2
> 0509010215   0.00    9.4  253.2   44.1   -999   26.4   -999   -999   10.6
>  0.2
>
> The data are recorded in 5 minute intervals and I would like to condense it
> into hourly means for "mph". For example,  I would like the hourly avg of
> mph so that the output would be as follows:
>
> Year Month Day Hour mph
> 2005 1 1 0 12
> 2005 1 1 1 7
> 2005 1 1 2 11, etc.
>
>
> It seems I am able to get the averages but not output the corresponding
> date/time stamp.  From looking at previous help questions, I think I need to
> us "ts" and "aggregate".   Gabor taught me how to convert the date/time
> stamp to an easier to manage format (his help is shown below).  This is what
> I have so far.
> library(zoo)
> library(chron)
>
> z <- read.zoo("SPL 2005 2008 met data 5 min wout full hdr.txt", header =
> TRUE, na.strings = -999,
> format = "%y%m%d%H%M", FUN = as.chron,
> colClasses = c("character", rep("numeric", 10)))
> z.ts <- ts(z, frequency=12)    #avging 5 min intervals to get hourly avg.
>   ww <- matrix(aggregate(z.ts[,2], FUN=mean))
>
> any thoughts as to how to add the time stamp is greatly welcomed!
>
> sherri heck
>
> ______________________________________________
> R-help at r-project.org mailing list
> 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.
>



More information about the R-help mailing list