Bchronplotevent: Histogram of the probability distribution for an event in a core.

Description

Provides a histogram and summary of the probability distribution for a given depth, slice, or event in a core.

Usage

Bchronplotevent(Bchrondata, depth = NULL, slice = NULL, event = NULL, eventname = NULL, nbreaks = 30, histcolour = "light blue", ...)

Arguments

Bchrondata

A Bchron data object with fields for the the input and output files, the name of the core being used, the details of the calibration curve, and details of the model run required. See Bchronload and

depth

One of the output depths (in cm) in the core (see Bchronload for more information on output depths).

slice

A slice number for one of the depths in the core.

event

An event name. Must have an associated corenameeventdepthseventname.txt file associated with it.

eventname

The full name of the event (eg "Alnus rise").

nbreaks

The number of breaks for the histogram.

histcolour

The colour of the histogram.

...

Other arguments, such as the width and height of the plot.

Details

This function can be used to identify the age of a certain part of the core according to the chronology model. It does not necessarily have to be an age that has a radiocarbon (or non-radiocarbon) date. Where an event name is given, a file should be created giving the depths at which the event occurred. See Bchronpredictevent for more information.

References

Haslett and Parnell (2008, JRSSC, 57:399-418) and Parnell et al (2008, QSR, 27(19-20):1872-1885)

Examples

Run this code

# For the example data, the following should run (when the files are placed in the default directory in windows; c:/Bchron)
# myrun <- Bchronload("Glendalough")
# Bchroncalibrate(myrun)
# BchronMCMC(myrun)
# Bchronpredictevent(myrun,event="Alnus")
# Bchronplotevent(myrun,event="Alnus",eventname="Alnus rise")