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).
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),
legend.loc = c(0.95, 0.02),
legend.size = c(0.05, 0.2),
padding = 0.1,
warn = TRUE
)A plot and the relevant dR values.
The southern limit of the rectangular region.
The western limit of the rectangular region.
The northern limit of the rectangular region.
The eastern limit of the rectangular region.
The variable to be plotted as colour. Expects a continuous variable. Defaults to 'dR'.
Whether or not to use a rainbow scale to plot the variable.
Size of the symbols. Defaults to 2.
Resolution of the map. Can be "small" 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 "small" if 'rnaturalearthhires' is not installed, and to "large" if it is installed.
Colour for minimum values.
Colour for maximum values.
The variable to be plotted as symbol. Expects a categoric variable. Defaults to 'feeding'.
Whether or not to plot the legend for the symbols.
Colour for the oceans. Defaults to ocean.col="aliceblue".
Colour for the land. Defaults to semi-transparent darkgreen: land.col=rgb(0, 0.5, 0, 0.6).
Location of the legend, if using a basic plot. Defaults to the bottom right corner based on par("usr"), legend.loc=c(0.95, 0.02)
Size of the legend, if using a basic plot. Defaults to legend.size=c(0.05, 0.2)
Area around the map if using a basic plot. Avoids strange line features. Defaults to padding=0.1.
Whether or not to warn if some recommended are not available.
This function uses the `rnaturalearth` package for country maps. If the high-resolution maps are desired, the `rnaturalearthhires` package must be installed from GitHub.
N_UK <- map.shells(53, -11, 60, 2, mapsize="small")
mean(N_UK$dR)
Run the code above in your browser using DataLab