Monthly anomalies
This function computes monthly anomalies.
intFun.anom.mly(mly, clim)
mly |
An R object with a monthly time series |
clim |
An R object with climatological monthly mean data |
A raster object with monthly anomalies
# make some data
month <- seq(1,12,1)
month <- rep(month,10)
mod <- runif(length(month), 0,10)
mod <- data.frame(month, mod)
# make an index
index <- list(mod$month)
clim.mly <- apply(mod, 2, function(x) {tapply(x, index, mean, na.rm=TRUE)})
mod.anom <- intFun.anom.mly(mod, clim.mly)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.