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,
scale = c(),
mincol = "yellow",
maxcol = "red",
symbol = "feeding",
symbol.legend = TRUE,
ocean.col = "aliceblue",
land.col = rgb(0, 0.5, 0, 0.6)
)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", "medium" or "large". If the latter, a high-resolution dataset will have to be downloaded using the R package 'rnaturalearthhires'. Since this package is not on CRAN, you will have to download it yourself. Defaults to 'medium' if 'rnaturalearthhires' is not installed, and to 'high' 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).
N_UK <- map.shells(53, -11, 60, 2, scale="medium")
mean(N_UK$dR)
Run the code above in your browser using DataLab