plot.CalGrid: Plot a summed probability distribution (from a CalGrid...

In ahb108/rcarbon: Calibration and Analysis of Radiocarbon Dates

Plot a summed probability distribution (from a CalGrid object)

Description

Plot a summed radiocarbon probability distribution. This is a basic function for plotting SPDs that have been constructed manually or by calibrating a summed or otherwise irregular CRA grid. In most instances, it is sensible to use plot.CalSPD instead.

Usage

## S3 method for class 'CalGrid'
plot(
  x,
  runm = NA,
  calendar = "BP",
  fill.p = "grey50",
  border.p = NA,
  xlim = NA,
  ylim = NA,
  cex.lab = 0.75,
  cex.axis = cex.lab,
  mar = c(4, 4, 1, 3),
  add = FALSE,
  ...
)

Arguments

x	A "CalGrid" class object of summed probabilities per calendar year BP.
runm	A number indicating the window size of the moving average to smooth the SPD. If set to NA no moving average is applied. Default is NA
calendar	Either 'BP' or 'BCAD'. Indicate whether the calibrated date should be displayed in BP or BC/AD. Default is 'BP'.
fill.p	Fill colour of the polygon depicting the summed probability distribution.
border.p	Border colour of the polygon depicting the summed probability distribution.
xlim	the x limits of the plot. In BP or in BC/AD depending on the choice of the parameter calender. Notice that if BC/AD is selected BC ages should have a minus sign (e.g. c(-5000,200) for 5000 BC to 200 AD).
ylim	Adjust y axis limits (otherwise sensible default).
cex.lab	Size of label text.
cex.axis	Size of axis text.
mar	Adjust margins around plot.
add	Whether or not the new graphic should be added to an existing plot.
...	Additional arguments affecting the plot

See Also

spd; plot.CalSPD

Examples

data(euroevol)
mycaldates <- calibrate(euroevol[1:10,"C14Age"], euroevol[1:10,"C14SD"], normalised=FALSE)
myspd <- spd(mycaldates, timeRange=c(8000,2000))
plot(myspd) #ordinary plot using \code{plot.CalSPD}
plot(myspd$grid) #working plot using the internal CalGrid object

ahb108/rcarbon documentation built on Aug. 29, 2023, 9:21 a.m.