Bchronproxyplot {Bchron}R Documentation

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 = 10, col=c("red","red","red"),smooth=FALSE,...)

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 Bchronmenu for more details.
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.
col Colours for the mean and 95pc prediction intervals.
smooth Whether kernel smoothing should be run on the proxy time series.
... 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), together with an estimate of the 95 uncertainty in the proxy series. As the uncertainty is in the horizontal axis, uncertainty bounds can sometimes cross each other if there is an inflexion point.

Author(s)

Andrew Parnell <Andrew.Parnell@ucd.ie>

References

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

See Also

Bchronplotevent

Examples

# 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)

[Package Bchron version 3.1 Index]