paleoMAS (version 2.0-1)

chron: Age-depth computation

Description

This function computes and returns interpolated ages for specified sample depths based on a matrix with depths and calibrated ages.

Usage

chron(age, depths, max.depth, type = "l", linear = TRUE, method = "fmm", ci = TRUE, dates = "points" , length = 0.05, digits = 2)

Arguments

age
A numeric matrix with three columns: depths of dated intervals, calibrated ages and standard deviation. The date or suspected date of the uppermost sample MUST be provided (usually 0 cm are present).
depths
A vector with a list of depths for which ages will be interpolated. The depths MUST be in ascendant order. Repeated depths are not allowed (see details).
max.depth
Maximum depth of the stratigraphic sequence (including analyzed and dated depths). It is used to extrapolate ages when the bottom of the sequence goes beyond the deepest dated interval.
type
Graphic parameter for the output plot. See plot for details.
linear
Interpolation method desired. If TRUE, ages are calculated using linear interpolation. If FALSE, a spline is fitted to the dataset.
method
Argument of the function spline. This must be one of "fmm", "natural", "periodic" or "monoH.FC". Only relevant if linear=FALSE. See spline for details.
ci
A logical value indicating whether confidence intervals are desired in the output plot.
dates
Graphical argument. "points" to display dates as red points, or "bars" to display dates as error bars spanning the provided standard deviation of calibrated ages.
length
Graphical parameter. Length in inches for the horizontal bar heads.
digits
Number of digits desired for interpolated ages.

Value

A list containing:
base
A matrix with the initial chronology.
chronology
A matrix with interpolated ages up to max depth and with deviation if ci=TRUE. The deviation is calculated according to the error provided by the user, therefore it does not necessarily correspond to 0.95 intervals.

Details

chron calculates interpolated ages using linear interpolation or spline as method. It is based on a three-column matrix with depth of dated interval, calibrated age, and calibrated standard deviation. The interpolated ages are calculated for a set of depths (defined by the user in depths). Whenever the stratigraphic sequence goes beyond the last dated point, the age for the bottom samples is extrapolated. max.depth is the maximum depth of the stratigraphic sequence, including dates and sampled points. Confidence intervals correspond to the provided chronology deviation. Therefore, it does not necessarily correspond to 0.95 interval.

If there is need to repeat a depth (e.g the floor and roof of a hiatus), the two numbers must be differentiated by at least decimal positions.

See Also

See plot for details on graphic parameters, approx for details on linear interpolation, and spline for details on non-linear interpolation.

Examples

Run this code
data(quexilchron,quexildepths)
#linear interpolation
chron(quexilchron,quexildepths,max.depth=1957,dates="points")
# Spline estimation
chron(quexilchron,quexildepths,max.depth=1957,dates="points")

Run the code above in your browser using DataCamp Workspace