SpatialVx (version 0.1-2)

locmeasures2dPrep: Set up a list object for locmeasures2d

Description

Sets up a list object to run with locmeasures2d.

Usage

locmeasures2dPrep(Fcst.name, Vx.name, thresholds = NULL, k = NULL, alpha = 0.1, bdconst = NULL, p = 2, loc = NULL, qs = NULL, units = NULL)

Arguments

Fcst.name

character giving the name of a forecast field.

Vx.name

character giving the name of the verification field.

thresholds

numeric vector giving the threhsold values over which (inclusive) binary fields are created. If NULL, quantiles for each field are calculated. Specifically the 0, 0.1, 0.25, 0.33, 0.5, 0.66, 0.75, 0.9 and 0.95 quantiles.

k

numeric vector for use with the partial Hausdorff distance. For k that are whole numerics or integers >= 1, then the k-th highest value is returned by locmeasures2d. If 0

alpha

numeric giving the alpha parameter for Pratt's Figure of Merit (FOM). See the help file for locperf for more details.

bdconst

numeric giving the cut-off value for Baddeley's delta metric.

p

numeric vector giving one or more values for the parameter p in Baddeley's delta metric. Usually this is just 2.

loc

(optional) n X 2 matrix giving the lon/lat coordinates of the grid points (currently not used).

qs

character vector giving alternative names for the thresholds. If 'thresholds' is NULL and 'qs' is NULL, then this will give the quantiles used.

units

character describing the units for the thresholds.

Value

• A list object withcomponents:
• Fcst.name, Vx.name
  characters, same as the resp. arguments.
• thresholds
  n X 2 matrix of thresholds.
• qs
  same as argument above.
• loc
  same as argument above.
• units
  same as argument above.
• xdim
  numeric vector of length 2 giving the dimensions of the forecast field (should be the same as the verification field).
• bdconst
  same as argument above.
• p
  same as argument above.
• k
  same as argument above.

See Also

locmeasures2d

Examples

data( geom000)
  data( geom001)
  hold <- locmeasures2dPrep( "geom001", "geom000", thresholds=c(0.1,50.1), k=4, units="in/100")