This function computes the inter-annual variability. All data must start in January. All months after the last Dec will be dropped if data does not end in December.
intFun.iav(anom)
R object with monthly anomalies
R object with inter-annual variability
# NOT RUN {
# make some data
month <- seq(1,12,1)
month <- rep(month,10)
mod <- runif(length(month), 0,100)
mod <- data.frame(month, mod)
# compute climatological mean monthly values
index <- list(mod$month)
mod.clim.mly <- apply(mod, 2, function(x) {tapply(x, index, mean, na.rm=TRUE)})
mod.anom <- intFun.anom.mly(mod, mod.clim.mly)
mod.iav <- intFun.iav(mod.anom)
# }
Run the code above in your browser using DataLab