modelPlot: Plot Growth Models

In nimbleCarbon: Bayesian Analyses of Radiocarbon Dates with NIMBLE

Plot Growth Models

Description

Plots growth models based on user provided parameters for prior and posterior predictive checks.

Usage

modelPlot(
  model,
  a,
  b,
  params,
  type = c("spaghetti"),
  nsample = NULL,
  interval = 0.9,
  calendar = "BP",
  col = "lightgrey",
  alpha = 0.1,
  ylim = NULL,
  xlim = NULL,
  xlab = NULL,
  ylab = NULL,
  add = FALSE,
  lwd = 1,
  ...
)

Arguments

model	growth model.
a	lower (earliest) limit of the distribution (in BP).
b	upper (latest) limit of the distribution (in BP).
params	a list of vectors containing model parameters. The names attribute of each vector should match growth model parameters.
type	either a 'spaghetti' plot or a quantile based 'envelope' plot. Default is 'spaghetti'.
nsample	number of samples to be used. Default is the length of the parameter vectors supplied in the argument params.
interval	quantile interval used for the envelope plot. Ignored when type is set to 'spaghetti'.
calendar	either 'BP' or 'BCAD'. Indicate whether the calibrated date should be displayed in BP or BC/AD. Default is 'BP'.
col	fill color for the quantile envelope (when type=='envelope') or line colour (when type=='spaghetti').
alpha	transparency value for each line in the spaghetti plot or the fill color in the 'envelope' plot. Default is 1.
ylim	the y limits of the plot.
xlim	the x limits of the plot (in Cal BP).
xlab	a label for the x axis. Default is 'Years cal BP','Years BC/AD','Years BC', or 'Years AD' depending on data range and settings of calendar.
ylab	a label for the y axis. Default is 'Probability'.
add	whether or not the new graphic should be added to an existing plot.
lwd	line width. Default is 1.
...	additional arguments affecting the plot

Value

None.

Examples

params = list(k=runif(100,0.01,0.02),r=runif(100,0.003,0.004))
modelPlot(model=dLogisticGrowth,a=5000,b=2000,params=params,type=c('spaghetti'),alpha=0.5)

nimbleCarbon documentation built on Aug. 14, 2023, 5:08 p.m.