plot.SpdModelTest function

Plot result of Monte-Carlo simulation of observed versus modelled SPDs

The function visualises the observed summed probability distribution of radiocarbon dates along with a simulation envelope for the null model and regions of positive and negative deviation.

## S3 method for class 'SpdModelTest'
plot(
  x,
  calendar = "BP",
  type = "spd",
  ylim = NA,
  xlim = NA,
  col.obs = "black",
  lwd.obs = 0.5,
  xaxs = "i",
  yaxs = "i",
  bbty = "f",
  bbtyRet = TRUE,
  drawaxes = TRUE,
  ...
)

Arguments

• x: A SpdModelTest class object generated using the modelTest function.
• calendar: Either 'BP' or 'BCAD'. Indicate whether the calibrated date should be displayed in BP or BC/AD. Default is 'BP'.
• type: Whether to display SPDs ('spd') or rates of change ('roc'). Default is 'spd'.
• ylim: the y limits of the plot.
• 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).
• col.obs: Line colour for the observed SPD.
• lwd.obs: Line width for the observed SPD.
• xaxs: The style of x-axis interval calculation (see par)
• yaxs: The style of y-axis interval calculation (see par)
• bbty: Display options; one between 'b','n',and 'f'. See details below.
• bbtyRet: Whether details on the intervals of positive and negative deviations are returned. Default is TRUE.
• drawaxes: A logical value determining whether the axes should be displayed or not. Default is TRUE.
• ...: Additional arguments affecting the plot

Details

The argument bbty controls the display options of the Monte-Carlo Test. Default settings (bbty='f') displays the observed SPD (solid black line), the simulation envelope of the fitted model (shaded grey polygon) and regions of significance positive (red semi-transparent rectangle) and negative (blue semi-transparent rectangle) deviation. The option bbty='b' removes the regions of positive/negative deviations, whilst the option bbty='n' displays the simulation envelope on existing plot.

See Also

modelTest