Calculate the distribution of age estimates of single or multiple depths.
Bacon.hist(d, set = get("info"), BCAD = set$BCAD, yr.lab = c(),
yr.lim = c(), hist.lab = "Frequency", hist.lim = c(), draw = TRUE,
prob = set$prob, hist.col = grey(0.5), hist.border = grey(0.2),
range.col = "blue", med.col = "green", mean.col = "red")
The depth or depths for which a histogram and age ranges should be provided. If multiple depths are given, then just the age ranges, median and means (no graphs) are provided for each depth.
Detailed information of the current run, stored within this session's memory as variable info.
The calendar scale of graphs and age output-files is in cal BP
by default, but can be changed to BC/AD using BCAD=TRUE
.
The labels for the calendar axis (default yr.lab="cal BP"
or "BC/AD"
if BCAD=TRUE
).
Minimum and maximum calendar age ranges, calculated automatically by default (yr.lim=c()
).
The y-axis is labelled ylab="Frequency"
by default.
Limits of the y-axis.
Draw a plot or not. Defaults to draw=TRUE
, however no plots are made if more than one depth d
is provided.
If draw=FALSE
, then the age ranges, median and mean are given for each depth (as four columns).
Age ranges are given as quantiles, e.g., 2.5% and 97.5% for the default of 95% confidence limits (prob=0.95
)).
Colour of the histogram. Default grey, hist.col=grey(0.5)
.
Colour of the histogram's outline. Default dark grey, hist.border=grey(0.2)
.
Colour of confidence ranges. Defaults to range.col="blue"
.
Colour of the median. Defaults to med.col="green"
.
Colour of the mean. Defaults to mn.col="red"
.
A plot with the histogram and the age ranges, median and mean, or just the age ranges, medians and means if more than one depth d
is given.
Age estimates of specific depths can also be plotted.
Blaauw, M. and Christen, J.A., Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6 (2011), no. 3, 457--474. https://projecteuclid.org/download/pdf_1/euclid.ba/1339616472
# NOT RUN {
Bacon(run=FALSE, coredir=tempfile())
agedepth(yr.res=50)
Bacon.hist(20)
Bacon.hist(20:30)
# }
Run the code above in your browser using DataLab