draw.ccurve: Draw a calibration curve.

In rice: Radiocarbon Equations

Draw a calibration curve.

Description

Draw one or two of the calibration curves, or add a calibration curve to an existing plot.

Usage

draw.ccurve(
  cal1 = c(),
  cal2 = c(),
  cc1 = "IntCal20",
  cc2 = NA,
  cc1.postbomb = FALSE,
  cc2.postbomb = FALSE,
  BCAD = FALSE,
  timescale = "C14",
  as.F = FALSE,
  as.pMC = FALSE,
  as.Delta = FALSE,
  timescale2 = c(),
  cal.lab = c(),
  cal.rev = FALSE,
  c14.lab = c(),
  cc2.lab = c(),
  c14.lim = c(),
  c14.rev = FALSE,
  ka = FALSE,
  add.yaxis = FALSE,
  cc1.col = rgb(0, 0, 1, 0.5),
  cc1.fill = rgb(0, 0, 1, 0.2),
  cc2.col = rgb(0, 0.5, 0, 0.5),
  cc2.fill = rgb(0, 0.5, 0, 0.2),
  add = FALSE,
  bty = "l",
  mar = c(),
  mgp = c(),
  cc.dir = NULL,
  legend = "topleft",
  ...
)

Arguments

cal1	First calendar year for the plot. Defaults to 0 cal BP.
cal2	Last calendar year for the plot. Defaults to 55,000 cal BP.
cc1	Name of the calibration curve. Can be "IntCal20", "Marine20", "SHCal20", or for the previous curves "IntCal13", "Marine13" or "SHCal13". Can also be "nh1", "nh2", "nh3", "sh1-2", "sh3", "nh1_monthly", "nh1_monthly", "nh2_monthly", "nh3_monthly", "sh1-2_monthly", "sh3_monthly", "Kure", "LevinKromer" or "Santos" for postbomb curves.
cc2	Optional second calibration curve to plot. Can be "IntCal20", "Marine20", "SHCal20", or for the previous curves "IntCal13", "Marine13" or "SHCal13". Defaults to nothing, NA.
cc1.postbomb	Use postbomb=TRUE to get a postbomb calibration curve for cc1 (default cc1.postbomb=FALSE).
cc2.postbomb	Use postbomb=TRUE to get a postbomb calibration curve for cc2 (default cc2.postbomb=FALSE).
BCAD	The calendar scale of graphs and age output-files is in cal BP (calendar or calibrated years before the present, where the present is AD 1950) by default, but can be changed to cal BC/AD using BCAD=TRUE.
timescale	Which timescale of radiocarbon to use. Defaults to timescale="C14" but can also be set to timescale="F14C", timescale="pMC" or timescale="Delta14C". Can be shorted to, respectively, "C", "F", "P" or "D" (or their lower-case equivalents). Alternatively, the timescale can be defined using 'as.F=TRUE', 'as.pMC=TRUE' or 'as.Delta=TRUE'.
as.F	Plot as F14C values. Defaults to as.F=FALSE.
as.pMC	Plot as pMC values. Defaults to as.pMC=FALSE.
as.Delta	Plot as Delta14C values. Defaults to as.Delta=FALSE.
timescale2	Which timescale to use for the second calibration curve (if used). Defaults to timescale2="C14" but can also be set to timescale2="F14C", timescale2="pMC" or timescale2="Delta14C". Can be shorted to, respectively, "C", "F", "P" or "D" (or their lower-case equivalents).
cal.lab	The labels for the calendar axis (default age.lab="cal BP" or "BC/AD" if BCAD=TRUE), or to age.lab="kcal BP" etc. if ka=TRUE.
cal.rev	Reverse the calendar axis.
c14.lab	Label for the C-14 axis. Defaults to 14C BP (or 14C kBP if ka=TRUE).
cc2.lab	Label for the righthand axis (if present). Defaults to the chosen timescale.
c14.lim	Axis limits for the C-14 axis. Calculated automatically by default.
c14.rev	Reverse the C-14 axis.
ka	Use kcal BP (and C14 kBP).
add.yaxis	Whether or not to plot the second calibration. Defaults to add.yaxis=FALSE.
cc1.col	Colour of the calibration curve (outline).
cc1.fill	Colour of the calibration curve (fill).
cc2.col	Colour of the calibration curve (outline), if activated (default cc2=NA).
cc2.fill	Colour of the calibration curve (fill), if activated (default cc2=NA).
add	Whether or not to add the curve(s) to an existing plot. Defaults to FALSE, which draws a new plot.
bty	Draw a box around a box of a certain shape. Defaults to bty="l".
mar	Plot margins (amount of white space along edges of axes 1-4). Defaults to give more white space if a second y-axis is to be plotted.
mgp	Axis text margins (where should titles, labels and tick marks be plotted).
cc.dir	Directory of the calibration curves. Defaults to where the package's files are stored (system.file), but can be set to, e.g., cc.dir="curves".
legend	Location of the legend (only activated if more than one curve is plotted). Plotted in the topleft corner by default. Use legend=c() to leave empty
...	Any additional optional plotting parameters.

Value

A plot of the calibration curve

Examples

draw.ccurve()
draw.ccurve(1000, 3000, cc2="Marine20")
draw.ccurve(1800, 2020, BCAD=TRUE, cc2="nh1", cc2.postbomb=TRUE)
draw.ccurve(1800, 2010, BCAD=TRUE, cc2="nh1", add.yaxis=TRUE)

rice documentation built on Sept. 9, 2025, 5:42 p.m.