Calculate correlation between forecasts and observations for an ensemble forecast, including an adjustment for finite ensemble sizes

`EnsCorr(ens, obs, R.new = NA)`

ens

a N*R matrix representing N time instances of real-valued R-member ensemble forecasts

obs

a numeric vector of length N with real-valued observations

R.new

positive number, can be `Inf`

, ensemble size for which correlation skill should be estimated, default is NA for using the actual size R of the ensemble

A vector with 4 entries:

`cmy`

: Correlation skill of the ensemble mean forecast`cmy_adj`

: Correlation skill of the ensemble mean forecast adjusted to ensemble size R.new`cxx`

: Average correlation between ensemble members`cxy`

: Average correlation between individual ensemble members and observation

Von Storch, Zwiers (2001): Statistical analysis in climate research. Cambridge University Press.

Murphy (1990), Assessment of the practical utility of extended range ensemble forecasts, Q. J. R. Meteorol. Soc., 116, 89-125.

Corr, CorrDiff

# NOT RUN { data(eurotempforecast) EnsCorr(ens, obs, R.new=Inf) # }