s2dverification (version 2.4.0)

PlotAno: Plot Raw Or Smoothed Anomalies

Description

Plots timeseries of raw or smoothed anomalies of any index output from Load() or Ano() or or Ano_CrossValid() or Smoothing() and organized in matrices with dimensions: c(nmod/nexp, nmemb/nparam, nsdates, nltime) for the model data c(nobs, nmemb, nsdates, nltime) for the observational data

Usage

PlotAno(exp_ano, obs_ano = NULL, sdates, toptitle = c("", "", "", "", "", "", 
        "", "", "", "", "", "", "", "", ""), ytitle = c("", "", "", "", "", 
        "", "", "", "", "", "", "", "", "", ""), limits = NULL, 
        legends = NULL, freq = 12, biglab = FALSE, fill = TRUE, memb = TRUE, 
        ensmean = TRUE, linezero = FALSE, points = FALSE, vlines = NULL, 
        fileout = c("output1_plotano.eps", "output2_plotano.eps", 
        "output3_plotano.eps", "output4_plotano.eps", "output5_plotano.eps"), 
        sizetit = 1)

Examples

# Load sample data as in Load() example:
example(Load)
clim <- Clim(sampleData$mod, sampleData$obs)
ano_exp <- Ano(sampleData$mod, clim$clim_exp)
ano_obs <- Ano(sampleData$obs, clim$clim_obs)
runmean_nb_months <- 12
dim_to_smooth <- 4  # Smooth along lead-times
smooth_ano_exp <- Smoothing(ano_exp, runmean_nb_months, dim_to_smooth)
smooth_ano_obs <- Smoothing(ano_obs, runmean_nb_months, dim_to_smooth)
PlotAno(smooth_ano_exp, smooth_ano_obs, startDates, 
        toptitle = paste('smoothed anomalies'), ytitle = c('K', 'K', 'K'), 
        legends = 'ERSST', biglab = FALSE, fileout = 'tos_ano.eps')