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

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

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

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

`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))

EnsBrier, EnsRps, EnsCrps, ScoreDiff, SkillScore

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