function to calculate Ozone 8-hour moving average for a data.frame
Examples
model_file <- paste(system.file("extdata", package = "eva3dm"),
"/model_o3_ugm3_36km.Rds", sep="")
model <- readRDS(model_file)
model_8h <- ma8h(model)
#> processing hourly data ...
#> processing 8-hour moving avarage ...
plot(model$date,model$Campinas, pch = 19,
main = expression(O[3]~~'['*mu*g*m^-3*']'))
points(model_8h$date,model_8h$Campinas, col = 'blue', pch = 19)
legend('topleft',bty = 'n',
pch = 19,
legend = c('hourly','8h-mov average'),
col = c('black','blue'))