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