Bchronproxyplot: Plot proxy data with associated chronological unceratinties

Description

Function produces plots of a proxy data time series, with associated error bands to indicate chronological uncertainty.

Usage

Bchronproxyplot(Bchrondata, proxy, title = NULL, xlabel = "Age (k cal yrs BP)", ylabel = "Proxy", num = 1000, col = "blue", meancol="red",transp=0.1,...)

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

proxy

A vector of proxy observations of the same length as the number of output depths in the core. See Bchronload for more information on the output depths.

title

A text string giving the title for the plot.

xlabel

A text string giving the x-axis label for the plot

ylabel

A text string giving the y-axis label for the plot

num

The number of sample chronologies to plot, default set at 1000

col

The colour of the lines drawn

meancol

The colour of the line drawn to represent the mean chronology

transp

The degree of transparancy of individual chronology lines drawn

...

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

Details

The plots produces show a number of sample proxy series based on different chronological samples (the number of these is controlled by num).

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")
# x <- sin(seq(0,pi,length=150))
# Bchronproxyplot(myrun,x)