fss2dPlot: Create Several Graphics for List Objects Returned from...

In SpatialVx: Spatial Forecast Verification

Create Several Graphics for List Objects Returned from hoods2d

Description

Creates several graphics for list objects returned from hoods2d. Mostly quilt and matrix plots for displaying results of smoothing fields over different neighborhood lengths and thresholds.

Usage

fss2dPlot(x, ..., matplotcol = 1:6, mfrow = c(1, 2), add.text = FALSE)

upscale2dPlot(object, args, ..., type = c("all",
                 "gss", "ts", "bias", "rmse"))

Arguments

x	list object with components fss, fss.random and fss.uniform. Effectively, it does the same thing as hoods2dPlot, but adds the fss.random and fss.uniform horizontal lines to the matrix plot.
object	list object with named components: rmse (numeric vector), ets, ts and bias all matrices whose rows represent neighborhood lengths, and whose columns represent thresholds.
args	list object passed to hoods2dPlot, see its help file for more details.
mfrow	mfrow parameter (see help file for par). If NULL, then the parameter is not re-set.
add.text	logical, if TRUE, FSS values will be added to the quilt plot as text (in addition to the color).
type	character string stating which plots to make (default is “all”).
...	Optional arguments to image and image.plot for fss2dPlot, and optional arguments to hoods2dPlot for upscale2dPlot
matplotcol	col argument to function matplot.

Details

makes quilt and matrix plots for output from hoods2d.

Value

No value is returned. A series of plots are created. It may be useful to use this function in conjunction with pdf in order to view all of the plots. See the help file for hoods2dPlot to plot individual results.

Author(s)

Eric Gilleland

See Also

hoods2dPlot, matplot, image, image.plot, hoods2d, pdf

Examples

##
## This is effectively an internal function, so the example is commented out
## in order for R's check to run faster.
##
## Not run: 
data( "geom001" )
data( "geom000" )
data( "ICPg240Locs" )

hold <- make.SpatialVx( geom000, geom001, thresholds = c(0.01,50.01),
    loc = ICPg240Locs, map = TRUE, projection = TRUE, loc.byrow = TRUE,
    units = "mm/h", data.name = "Geometric", obs.name = "observation",
    model.name = "case 1" )

look <- hoods2d(hold, levels=c(1, 3, 5, 33, 65),
    verbose=TRUE)
plot( look)

## End(Not run)

SpatialVx documentation built on May 29, 2024, 9:31 a.m.