s2dverification (version 2.5.0)

Ano_CrossValid: Computes Anomalies In Cross-Validation Mode

Description

This function computes anomalies from experimental and observational matrices output from load() by subtracting the climatologies computed in a cross-validation mode and with a per-pair method.

Usage

Ano_CrossValid(var_exp, var_obs, memb = TRUE)

Arguments

var_exp

Model data: c(nmod/nexp, nmemb/nparam, nsdates, nltime) up to c(nmod/nexp, nmemb/nparam, nsdates, nltime, nlevel, nlat, nlon)

var_obs

Observational data: c(nobs, nmemb, nsdates, nltime) up to c(nobs, nmemb, nsdates, nltime, nlevel, nlat, nlon)

memb

memb: TRUE/FALSE (1 climatology for each member/1 climatology averaging all the members). Default = TRUE.

Value

$ano_exp

Matrix with same dimensions as var_exp

$ano_obs

Matrix with same dimensions as var_obs

Examples

# Load sample data as in Load() example:
example(Load)
anomalies <- Ano_CrossValid(sampleData$mod, sampleData$obs)
PlotAno(anomalies$ano_exp, anomalies$ano_obs, startDates, 
        toptitle = paste('anomalies'), ytitle = c('K', 'K', 'K'), 
        legends = 'ERSST', biglab = FALSE, fileout = 'tos_ano_crossvalid.eps')