Plot grey-scale representation of modelled accumulation rates over a core's depth. Each section of the core (see Bacon's option "thick"
) will have modelled accumulation rates.
accrate.depth.ghost(set = get("info"), d = set$elbows, d.lim = c(),
acc.lim = c(), d.lab = c(), cmyr = FALSE, acc.lab = c(),
dark = 1, grey.res = 100, prob = 0.95, plot.range = TRUE,
range.col = grey(0.5), range.lty = 2, plot.mean = TRUE,
mean.col = "red", mean.lty = 2, rotate.axes = FALSE,
rev.d = FALSE, rev.acc = FALSE)
Detailed information of the current run, stored within this session's memory as variable info.
The depths for which the accumulation rates are to be calculated. Default to the entire core.
Axis limits for the depths.
Axis limits for the accumulation rates.
Label for the depth axis.
Accumulation rates can be calculated in cm/year or year/cm. By default cmyr=FALSE
and accumulation rates are calculated in year per cm. Axis limits are difficult to calculate when cmyr=TRUE
, so a manual adaptation of acc.lim
might be a good idea.
Axis label for the accumulation rate.
The darkest grey value is dark=1 by default; lower values will result in lighter grey but values >1 are not advised.
Grey-scale resolution. Default grey.res=100
.
Probability ranges. Defaults to prob=0.95
.
If plot.range=TRUE
, the confidence ranges (two-tailed; half of the probability at each side) are plotted.
Colour of the confidence ranges.
Line type of the confidence ranges.
If plot.mean=TRUE
, the means are plotted.
Colour of the mean accumulation rates.
Type of the mean lines.
The default is to plot the accumulation rates horizontally and the depth vertically (rotate.axes=FALSE
). Change rotate.axes value to rotate axes.
The direction of the depth axis can be reversed from the default (rev.d=TRUE
.
The direction of the accumulation rate axis can be reversed from the default (rev.acc=TRUE
).
A grey-scale plot of accumulation rate against core depth.
This plot shows the modelled accumulation rates in grey-scales, where darker grey indicates more likely accumulation rates.
Axis limits for accumulation rates are estimated automatically, however upper limits can be very variable (and thus hard to predict)
if calculated in cm/yr; therefore you might want to manually adapt the axis limits after plotting with default settings (e.g., acc.lim=c(0,1)
). See also accrate.age.ghost
, accrate.depth
and accrate.age
.
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/euclid.ba/1339616472
# NOT RUN {
Bacon(run=FALSE, coredir=tempfile())
agedepth(yr.res=50, d.res=50, d.by=10)
layout(1)
accrate.depth.ghost()
# }
Run the code above in your browser using DataLab