Description Usage Arguments Details Value Author(s) See Also Examples
A list object containing the verification sets of spatial verification and forecast fields with pertinent information.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20  make.SpatialVx(X, Xhat, thresholds = NULL, loc = NULL, projection =
FALSE, subset = NULL, time.vals = NULL, reg.grid =
TRUE, map = FALSE, loc.byrow = FALSE, field.type = "",
units = "", data.name = "", obs.name = "X", model.name
= "Xhat", q = c(0, 0.1, 0.25, 0.33, 0.5, 0.66, 0.75,
0.9, 0.95), qs = NULL)
## S3 method for class 'SpatialVx'
hist(x, ..., time.point = 1, obs = 1, model = 1,
threshold.num = NULL)
## S3 method for class 'SpatialVx'
plot( x, ..., time.point = 1, obs = 1, model = 1,
col = c( "gray", tim.colors( 64 ) ), zlim, mfrow = c(1, 2) )
## S3 method for class 'SpatialVx'
print(x, ...)
## S3 method for class 'SpatialVx'
summary(object, ...)

X 
An n X m matrix or n X m X T array giving the verification field of interest. If an array, T is the number of time points. 
Xhat 
An n X m matrix or n X m X T array giving the forecast field of interest, or a list of such matrices/arrays with each component of the list an n X m matrix or n X m X T array defining a separate forecast model. 
thresholds 
single numeric, numeric vector, or Nu X Nf matrix, where Nu are the number of thresholds and Nf the number of forecast models plus one (for the verification) giving the threshold values of interest for the verification set or components of the set. If NULL (default), then thresholds will be calculated as the quantiles (defined through argument 
loc 
If lon/lat coordinates are available, then this is an n * m X 2 matrix giving the lon/lat coordinates of each grid point or location. Should follow the convention used by the maps package. 
projection 
logical, are the grids projections onto the globe? If so, when plotting, it will be attempted to account for this by using the 
subset 
vector identifying which specific grid points should be included (if not all of them). This argument may be ignored by most functions and is included for possible future functionality. 
time.vals 
If more than one time point is available in the set (i.e., the set is of n X m X T arrays, with T > 1), then this argument can be used to define the time points. If missing, the default will yield the vector 
reg.grid 
logical, is the verification set on a regular grid? This is another feature intended for possible future functionality. Most functions in this package assume the set is on a regular grid. 
map 
logical, should the plot function attempt to place a map onto the plot? Only possible if the 
field.type, units 
character used for plot labelling and printing information to the screen. Describes what variable and in what units the field represents. 
data.name, obs.name, model.name 
character vector describing the verification set overall, the observation(s) and the model(s), resp. 
q 
numeric vector giving the values of quantiles to be used for thresholds. Only used if 
qs 
character vector describing the quantiles used. Again, only used if 
time.point 
numeric or character indicating which time point from the “SpatialVx” verification set to select for analysis. May also be a function name, in which case the function is applied at each grid point individually across time. 
obs, model 
numeric indicating which observation/forecast model to select for the analysis. 
col, zlim 
optional arguments to 
mfrow 
optional argument to change the mfrow argument for the graphic device. Default is one row with two plots (obs and model). If null, then the mfrow argument will not be changed. 
x, object 
list object of class “SpatialVx”. 
loc.byrow 
logical determining whether to set up the location matrices using 
threshold.num 
If not null, then the threshold index to apply a threshold to the fields before creating the histogram. 
... 

This function merely describes a spatial verification set that includes the actual data as well as numerous attributes that are used by several of the subsequent functions that might be employed. In many cases, the attribute information may be passed on to output from other functions for plot labelling and printing purposes (e.g., in order to identify the verification set, time point(s), etc.).
All (or perhaps most) subsequent functions in this package utilize objects of this class and the information contained in the attributes. This function simply gathers information and data sets into a particular form.
The plot method function attempts to create an image plot of each field in the set (at each time point). If projection is TRUE, then it will attempt to preserve the projection (via poly.image
of package fields). It will also add white contour lines showing the thresholds. If map is TRUE and loc
was supplied, then a map will also be added, if possible.
A list object with two (unnamed) components:
1 
matrix or array (same as input argument) giving the observation 
2 
Either a matrix or array (same as input argument) or a list of such objects if more than one forecast model. 
Several attributes are also included among the following:
xdim 
numeric of length 2 or 3 giving the dimensions of the verification set (i.e., m, n and T, if relevant). 
time 
vector giving the time values 
thresholds 
matrix giving the thresholds for each field. If there is more than one forecast, and they use the same threshold, this matrix may have only two columns. 
udim 
the dimensions of the thresholds matrix. 
loc 
nm X 2 matrix giving the locations. If loc was not given, this will be c(rep(1:n, m), rep(1:m, each=n)). 
subset 
If given, this is a numeric vector describing a subset of loc to be used. 
data.name, obs.name, model.name 
character vector giving the names of the data sets (same as input arguments). 
nforecast 
single numeric giving the number of different forecast models contained in the object. 
field.type, units 
character strings, same as input arguments. 
projection 
logical, is the grid a projection? 
reg.grid 
logical, is the grid a regular grid? 
map 
logical, should a map be added to image plots of the data? 
qs 
character vector giving the names of the threshold quantiles. 
msg 
A message involving the data name, field type and units for adding info to plots, etc. 
Eric Gilleland
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124  data( "UKobs6" )
data( "UKfcst6" )
data( "UKloc" )
hold < make.SpatialVx( UKobs6, UKfcst6, thresholds = c(0.01, 20.01),
loc = UKloc, field.type = "Precipitation", units = "mm/h",
data.name = "Nimrod", obs.name = "Observations 6", model.name = "Forecast 6",
map = TRUE)
hold
plot( hold )
hist( hold )
hist( hold, threshold.num = 2 )
## Not run:
# Stage 2 Observation files from
# NSSL/NCEP Spring 2005 Forecast Experiment.
data( "obs0426" )
data( "obs0513" )
data( "obs0514" )
data( "obs0518" )
data( "obs0519" )
data( "obs0525" )
data( "obs0601" )
data( "obs0603" )
data( "obs0604" )
st2 < array(c(c(obs0426), c(obs0513), c(obs0514), c(obs0518),
c(obs0519), c(obs0525), c(obs0601), c(obs0603), c(obs0604)),
dim=c(601, 501, 9))
rm(obs0426, obs0513, obs0514, obs0518, obs0519, obs0525, obs0601,
obs0603, obs0604)
# wrf2caps
data( "wrf2caps0425" )
data( "wrf2caps0512" )
data( "wrf2caps0513" )
data( "wrf2caps0517" )
data( "wrf2caps0518" )
data( "wrf2caps0524" )
data( "wrf2caps0531" )
data( "wrf2caps0602" )
data( "wrf2caps0603" )
# wrf4ncar
data( "wrf4ncar0425" )
data( "wrf4ncar0512" )
data( "wrf4ncar0513" )
data( "wrf4ncar0517" )
data( "wrf4ncar0518" )
data( "wrf4ncar0524" )
data( "wrf4ncar0531" )
data( "wrf4ncar0602" )
data( "wrf4ncar0603" )
# wrf4ncep
data( "wrf4ncep0425" )
data( "wrf4ncep0512" )
data( "wrf4ncep0513" )
data( "wrf4ncep0517" )
data( "wrf4ncep0518" )
data( "wrf4ncep0524" )
data( "wrf4ncep0531" )
data( "wrf4ncep0602" )
data( "wrf4ncep0603" )
wrf2caps < array(c(c(wrf2caps0425), c(wrf2caps0512),
c(wrf2caps0513), c(wrf2caps0517), c(wrf2caps0518),
c(wrf2caps0524), c(wrf2caps0531), c(wrf2caps0602),
c(wrf2caps0603)),
dim=c(601, 501, 9))
wrf4ncar < array(c(c(wrf4ncar0425), c(wrf4ncar0512),
c(wrf4ncar0513), c(wrf4ncar0517), c(wrf4ncar0518),
c(wrf4ncar0524), c(wrf4ncar0531), c(wrf4ncar0602),
c(wrf4ncar0603)),
dim=c(601, 501, 9))
wrf4ncep < array(c(c(wrf4ncep0425), c(wrf4ncep0512),
c(wrf4ncep0513), c(wrf4ncep0517), c(wrf4ncep0518),
c(wrf4ncep0524), c(wrf4ncep0531), c(wrf4ncep0602),
c(wrf4ncep0603)),
dim=c(601, 501, 9))
rm(wrf2caps0425, wrf2caps0512, wrf2caps0513, wrf2caps0517,
wrf2caps0518, wrf2caps0524, wrf2caps0531, wrf2caps0602,
wrf2caps0603, wrf4ncar0425, wrf4ncar0512, wrf4ncar0513,
wrf4ncar0517, wrf4ncar0518, wrf4ncar0524, wrf4ncar0531,
wrf4ncar0602, wrf4ncar0603, wrf4ncep0425, wrf4ncep0512,
wrf4ncep0513, wrf4ncep0517, wrf4ncep0518, wrf4ncep0524,
wrf4ncep0531, wrf4ncep0602, wrf4ncep0603)
fcst < list(wrf2caps, wrf4ncar, wrf4ncep)
rm(wrf2caps, wrf4ncar, wrf4ncep)
# Now, create the object.
data( "ICPg240Locs" )
ICPreal < make.SpatialVx( st2, fcst, thresholds = c(0.1, 20.1),
loc = ICPg240Locs, projection = TRUE, loc.byrow = TRUE,
time.vals = c(2005042600, 2005051300, 2005051400, 2005051800,
2005051900, 2005052500, 2005060100, 2005060300, 2005060400),
map = TRUE, field.type = "Precipitation", units = "mm/h",
data.name = "ICP Real Test Cases",
obs.name = "Stage II Analysis",
model.name = c( "WRF 2 CAPS", "WRF 4 NCAR", "WRF 4 NCEP"))
plot( ICPreal )
plot( ICPreal, time.point = mean )
plot( ICPreal, model = 2 )
plot( ICPreal, model = 3 )
plot( ICPreal, time.point = 2, model = 2 )
## End(Not run)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.