Quilt Plot and a Matrix Plot.

Share:

Description

Function to make a quilt plot and a matrix plot for a matrix whose rows represent neighborhood lengths, and whose columns represent different threshold choices.

Usage

1

Arguments

x

l by q numeric matrix.

args

list object with components: threshold (numeric vector giving the threshold values), qs (optional numeric vector giving the quantiles used if the thresholds represent quantiles rather than hard values), levels (numeric giving the neighborhood lengths (in grid squares) used, units (optional character giving the units for the thresholds)

...

optional arguments to image and image.plot functions. May not include xaxt, yaxt, lab, lab, col, or legend.only

Details

Used by plot.hoods2d, but can be useful for other functions. Generally, however, this is an internal function that should not be called by the user. However, it might be called instead of plot.hoods2d in order to make a subset of the available plots.

Value

No value is returned. A plot is created.

Author(s)

Eric Gilleland

References

Becker, R. A., Chambers, J. M. and Wilks, A. R. (1988) The New S Language. Wadsworth & Brooks/Cole.

See Also

matplot, image, image.plot, plot.hoods2d, hoods2d

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
x <- y <- matrix( 0, 50, 50)
x[ sample(1:50,10), sample(1:50,10)] <- rexp( 100, 0.25)
y[ sample(1:50,20), sample(1:50,20)] <- rexp( 400)
hold <- make.SpatialVx(x, y, thresholds=c(0.1, 0.5), field.type="random")
look <- hoods2d(hold, which.methods=c("multi.event", "fss"),
    levels=c(1, 3, 20))
hoods2dPlot( look$multi.event$hk, args=hold,
    main="Hanssen Kuipers Score (Multi-Event Cont. Table)")
## Not run: 
data(geom001)
data(geom000)
data(ICPg240Locs)
hold <- make.SpatialVx(geom000, geom001, thresholds=c(0.01,50.01),
    loc=ICPg240Locs, projection=TRUE, map=TRUE, loc.byrow = TRUE,
    field.type="Precipitation", units="mm/h",
    data.name=c("Geometric ICP Test Cases", "geom000", "geom001"))
look <- hoods2d(hold, levels=c(1, 3, 5, 17, 33, 65), verbose=TRUE)
par(mfrow=c(1,2))
hoods2dPlot(look$pragmatic$bss, args=attributes(hold))

## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.