[R] adding date/time stamp
Sherri Heck
sheck at ucar.edu
Wed Dec 31 00:30:52 CET 2008
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
More information about the R-help
mailing list