plot.SpdModelTest: Plot result of Monte-Carlo simulation of observed versus...

In rcarbon: Calibration and Analysis of Radiocarbon Dates

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

Description

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.

Usage

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

rcarbon documentation built on Aug. 24, 2023, 5:11 p.m.