plot.SpdModelTest {rcarbon}

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