afc.mm {afc} | R Documentation |
Routine to calculate the Generalized Discrimination Score (aka Two-Alternatives Forced Choice Score 2AFC) for the situation of polychotomous observations (ordinal) and polychotomous forecasts (ordinal)
afc.mm(obsv, fcst, mv = 3, mf = 3)
obsv |
vector with polychotomous observations (values in 1,..,mv) |
fcst |
vector of same length as obsv with polychotomous forecasts (values in 1,..,mf) |
mv |
number of observation categories (default = 3) |
mf |
number of forecast categories (default = 3) |
This routine applies Eq.14 of Mason and Weigel (2009) to calculate the 2AFC.
p.afc |
Value of Generalized Discrimination (2AFC) Score |
Andreas Weigel, Federal Office of Meteorology and Climatology, MeteoSwiss, Zurich, Switzerland
S.J. Mason and A.P. Weigel, 2009. A generic verification framework for administrative purposes. Mon. Wea. Rev., 137, 331-349
#Forecasts and observations of Nino-3.4 index
#Load set of polychotomous observations (4 categories) and polychotomous forecasts (4 categories)
data(cnrm.nino34.mm)
obsv = cnrm.nino34.mm$obsv
fcst = cnrm.nino34.mm$fcst
#Calculate skill score
afc.mm(obsv,fcst,4,4)