plot.CalGrid {rcarbon}

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