plot.CalGrid function

Plot a summed probability distribution (from a CalGrid object)

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.

## 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

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

See Also

spd; plot.CalSPD