plot.mcsimres: Plot Monte-Carlo simulation results on aoristic data

Description

Plot Monte-Carlo simulation based percentile intervals on frequency or rate of change of events.

Usage

# S3 method for mcsimres
plot(
  x,
  interval = 0.9,
  changexpr = expression((t1 - t0)/r),
  minortick = NULL,
  ylim = NULL,
  xlab = NULL,
  ylab = NULL,
  calendar = "BP",
  col = "black",
  lwd = 1,
  lty = 1,
  col.fill = "lightblue",
  pch = 20,
  type = "sum",
  plot.legend = TRUE,
  legend.arg = NULL,
  ...
)

Value

No return value (plot function)

Arguments

x: A mcsimres class object generated using the mcsim() function.
interval: A value between 0 and 1 defining the percentile interval. Default is 0.9.
changexpr: An expression for calculating the rate of change between abutting time-blocks. Available input options are t1 (the focal time-block), t0 (the previous time-block), r (the distance between t0 and t1, i.e. the time-block resolution), and any other standard constants and mathematical operators. Default is expression((t1-t0)/r). A possible alternative could be expression(log(t1/t0)/r).
minortick: Interval for minor ticks in the x-axis label. Default is estimated based on timescale.
ylim: Limits of the y-axis. Default estimated from posterior ranges.
xlab: Label for the x-axis. Default based on calendar.
ylab: Label for the y-axis. Default is "Probability Mass".
calendar: Either 'BP' or 'BCAD'. Indicate whether the x-axis should be displayed in BP or BC/AD. Default is 'BP'.
col: Color of Monte-Carlo simulation mean. Default is black.
lwd: Line width of Monte-Carlo mean. Default is 1.
lty: Line type Monte-Carlo mean. Default is 1.
col.fill: Fill color for the first (inner) percentile interval. Default is 'lightblue'.
pch: Point symbol used to display mean posteriors. Default is 20.
type: Determine whether to display total number of events (if set to 'sum') or the rate of change ('roc'), computed as (t0/t1)^(1/r)-1, where t0 is the number of events in given time-block t, t1 is the number of events of the next time-block t+1, and r is the size (in years) of the time-blocks. Defaults is 'sum'.
plot.legend: Logical indicating whether to display a legend or not (default is TRUE).
legend.arg: List containing arguments to be directed to the legend() function.
...: Additional arguments affecting the plot.