map.shells {rice} | R Documentation |
Plot regional shell-derived dR values
Description
Find the shells that fit within a rectangular region (bounded by N, E, S and W), and plot the dR values and feeding ecology. Uses the marine database downloaded (30 Aug 2024) from calib.org/marine. See Reimer PJ, Reimer RW, 2001. A marine reservoir correction database and on-line interface. Radiocarbon 43:461-3. Expects the coordinates for the map to be provided (starting south, then clockwise as with R axes).
Usage
map.shells(
S = 48,
W = -15,
N = 62,
E = 5,
colour = "dR",
rainbow = FALSE,
size = 2,
mapsize = "large",
mincol = "yellow",
maxcol = "red",
symbol = "feeding",
symbol.legend = TRUE,
ocean.col = "aliceblue",
land.col = rgb(0, 0.5, 0, 0.6)
)
Arguments
S |
The southern limit of the rectangular region. |
W |
The western limit of the rectangular region. |
N |
The northern limit of the rectangular region. |
E |
The eastern limit of the rectangular region. |
colour |
The variable to be plotted as colour. Expects a continuous variable. Defaults to 'dR'. |
rainbow |
Whether or not to use a rainbow scale to plot the variable. |
size |
Size of the symbols. Defaults to 2. |
mapsize |
Resolution of the map. Can be "small", "medium" or "large". If the latter, a high-resolution dataset will have to be downloaded using the R package 'rnaturalearthhires'. Since this package is on github but not on CRAN, you will have to download it yourself (using the command devtools::install_github("ropensci/rnaturalearthhires")). Defaults to 'medium' if 'rnaturalearthhires' is not installed, and to 'high' if it is installed. |
mincol |
Colour for minimum values. |
maxcol |
Colour for maximum values. |
symbol |
The variable to be plotted as symbol. Expects a categoric variable. Defaults to 'feeding'. |
symbol.legend |
Whether or not to plot the legend for the symbols. |
ocean.col |
Colour for the oceans. Defaults to |
land.col |
Colour for the land. Defaults to semi-transparent darkgreen: |
Details
This function uses the 'rnaturalearth' package for country maps. If the high-resolution maps are desired, the 'rnaturalearthhires' package must be installed from GitHub.
Value
A plot and the relevant dR values.
Examples
N_UK <- map.shells(53, -11, 60, 2, mapsize="medium")
mean(N_UK$dR)