set.seed(1)
data <- data.frame(date = seq(from = ISOdatetime(2012, 01, 01, 00, 00, 00),
length.out = 4320, by=5),
nox = sample(1:20, 4320, replace=TRUE))
hr.means <- aggregate(data["nox"],
list(hour = cut(data$date, breaks="hour")),
mean, na.rm = TRUE)
hr.means
c2V0LnNlZWQoMSkKZGF0YSA8LSBkYXRhLmZyYW1lKGRhdGUgPSBzZXEoZnJvbSA9IElTT2RhdGV0aW1lKDIwMTIsIDAxLCAwMSwgMDAsIDAwLCAwMCksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxlbmd0aC5vdXQgPSA0MzIwLCBieT01KSwKICAgICAgICAgICAgICAgICAgIG5veCA9IHNhbXBsZSgxOjIwLCA0MzIwLCByZXBsYWNlPVRSVUUpKQoKaHIubWVhbnMgPC0gYWdncmVnYXRlKGRhdGFbIm5veCJdLCAKICAgICAgICAgICAgICAgICAgICAgIGxpc3QoaG91ciA9IGN1dChkYXRhJGRhdGUsIGJyZWFrcz0iaG91ciIpKSwgCiAgICAgICAgICAgICAgICAgICAgICBtZWFuLCBuYS5ybSA9IFRSVUUpCmhyLm1lYW5z