EnsQs: Calculate the ensemble-adjusted Quadratic Score (QS) for...

In SpecsVerification: Forecast Verification Routines for Ensemble Forecasts of Weather and Climate

Description

Calculate the ensemble-adjusted Quadratic Score (QS) for categorical forecasts

Usage

EnsQs(ens, obs, R.new = NA)

FairQs(ens, obs)

Arguments

ens	a N*R matrix of integers, representing N time instances of categorical ensemble forecasts; ens[t,r] indicates the category index that the r-th ensemble member forecasts at time t
obs	a vector of length N, obs[t] is the category that occurred at time t
R.new	ensemble size for which the scores should be adjusted

Details

'FairQs(ens, obs)' returns 'EnsQs(ens, obs, R.new=Inf)'

It is assumed that the smallest class index is 1, and the largest class index is calculated by max(c(ens,obs))

Value

numeric vector of length N with the ensemble-adjusted quadratic score values

See Also

EnsBrier, EnsRps, EnsCrps, ScoreDiff, SkillScore

Examples

data(eurotempforecast)
EnsQs(ens.cat, obs.cat, R.new=Inf)

Example output

         1983          1984          1985          1986          1987 
 1.086957e-01  1.326087e+00  2.934783e-01  9.057971e-02  1.307971e+00 
         1988          1989          1990          1991          1992 
 3.260870e-01  7.101449e-01  1.108696e+00  1.086957e-01  8.695652e-01 
         1993          1994          1995          1996          1997 
 1.413043e+00  2.173913e-02  1.159420e-01  7.608696e-01  6.594203e-01 
         1998          1999          2000          2001          2002 
 5.797101e-02  1.086957e-01  1.047101e+00  2.318841e-01  1.557971e+00 
         2003          2004          2005          2006          2007 
 1.608696e+00 -3.035766e-18  3.985507e-01  3.333333e-01  3.985507e-01 
         2008          2009 
 1.086957e-01  5.797101e-02 

SpecsVerification documentation built on March 26, 2020, 7:55 p.m.