plotmap
Plot Maps
The function takes a list
polygons and draws the corresponding map. Different
colors for each polygon can be used.
- Keywords
- hplot
Usage
plotmap(map, x = NA, id = NULL, select = NULL,
legend = TRUE, names = FALSE, values = FALSE, ...)
Arguments
- map
The map to be plotted, usually an object that inherits from
SpatialPolygons
, but may also be a list of polygons, i.e., each list entry is a matrix with x- and y-coordinates.- x
A vector, data.frame or matrix. In the latter case
x
should have two columns, one column that is the region identifier and one that contains the values to be plotted. In casex
is a data.frame, the function searches for"character"
or"factor"
columns for the region identifier. Ifx
is a matrix, the second column ofx
is supposed to be the region identifier. Ifx = NA
andmap
is an object of class"SpatialPolygonsDataFrame"
only the polygons will be plotted without using the data.- id
If argument
x
is a vector, argumentid
should contain the region identifier vector.- select
Select the column of the data in
x
which should be used for plotting, may be an integer or character with the corresponding column name.- legend
Should a color legend be added to the plot, see also function
colorlegend
.- names
If set to
TRUE
the name for each polygon will also be plotted at the centroids of the corresponding polygons.- values
If set to
TRUE
the corresponding values for each polygon will also be plotted at the centroids of the polygons.- …
Arguments to be passed to
colorlegend
and others, e.g. change theborder
of the polygons anddensity
(mdensity
for missing regions inid
), seepolygon
.
See Also
Examples
# NOT RUN {
## Example from mgcv ?mrf.
## Load Columbus Ohio crime data (see ?columbus for details and credits).
data("columb", package = "mgcv")
data("columb.polys", package = "mgcv")
## Plot the map.
plotmap(columb.polys)
## Plot aggregated data.
a <- with(columb, aggregate(crime,
by = list("district" = district), FUN = mean))
plotmap(columb.polys, x = a$x, id = a$district)
plotmap(columb.polys, x = a$x, id = a$district,
pos = "topleft")
plotmap(columb.polys, x = a$x, id = a$district,
pos = "topleft", side.legend = 2)
plotmap(columb.polys, x = a$x, id = a$district,
pos = "topleft", side.legend = 2, side.ticks = 2)
plotmap(columb.polys, x = a$x, id = a$district,
pos = "topleft", side.legend = 2, side.ticks = 2,
col = heat_hcl, swap = TRUE)
plotmap(columb.polys, x = a$x, id = a$district,
pos = "topleft", side.legend = 2, side.ticks = 2,
col = heat_hcl, swap = TRUE, range = c(10, 50))
plotmap(columb.polys, x = a$x, id = a$district,
pos = "topleft", side.legend = 2, side.ticks = 2,
col = heat_hcl(5), swap = TRUE, range = c(10, 50),
lrange = c(0, 60))
# }