agespectrum: Plot a (40Ar/39Ar) release spectrum

Description Usage Arguments Details Value See Also Examples

View source: R/agespectrum.R

Description

Produces a plot of boxes whose widths correspond to the cumulative amount of ^{39}Ar (or any other variable), and whose heights express the analytical uncertainties. Only propagates the analytical uncertainty associated with decay constants and J-factors after computing the plateau composition.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
agespectrum(x, ...)

## Default S3 method:
agespectrum(x, alpha = 0.05, plateau = TRUE,
  random.effects = TRUE, levels = NA, clabel = "",
  plateau.col = c("#00FF0080", "#FF000080"),
  non.plateau.col = "#00FFFF80", sigdig = 2, line.col = "red",
  lwd = 2, title = TRUE, show.ci = TRUE,
  xlab = "cumulative fraction", ylab = "age [Ma]", hide = NULL, ...)

## S3 method for class 'ArAr'
agespectrum(x, alpha = 0.05, plateau = TRUE,
  random.effects = TRUE, levels = NA, clabel = "",
  plateau.col = c("#00FF0080", "#FF000080"),
  non.plateau.col = "#00FFFF80", sigdig = 2, exterr = TRUE,
  line.col = "red", lwd = 2, i2i = FALSE, hide = NULL, ...)

Arguments

x

a three-column matrix whose first column gives the amount of ^{39}Ar in each aliquot, and whose second and third columns give the age and its uncertainty.

OR

an object of class ArAr

...

optional parameters to the generic plot function

alpha

the confidence level of the error bars/boxes and confidence intervals.

plateau

logical flag indicating whether a plateau age should be calculated. If plateau=TRUE, the function will compute the weighted mean of the largest succession of steps that pass the Chi-square test for age homogeneity. If TRUE, returns a list with plateau parameters.

random.effects

if TRUE, computes the weighted mean using a random effects model with two parameters: the mean and the dispersion. This is akin to a ‘model-3’ isochron regression.

if FALSE, attributes any excess dispersion to an underestimation of the analytical uncertainties. This akin to a ‘model-1’ isochron regression.

levels

a vector with additional values to be displayed as different background colours of the plot symbols.

clabel

label of the colour legend

plateau.col

a vector of two fill colours of the rectangles used to mark the steps belonging to the age plateau. If levels=NA, then only the first colour is used. If levels is a vector of numbers, then bg is used to construct a colour ramp.

non.plateau.col

if plateau=TRUE, the steps that do NOT belong to the plateau are given a different colour.

sigdig

the number of significant digits of the numerical values reported in the title of the graphical output (only used if plateau=TRUE).

line.col

colour of the average age line

lwd

width of the average age line

title

add a title to the plot?

show.ci

show a 100(1-α)% confidence interval for the plateau age as a grey band

xlab

x-axis label

ylab

y-axis label

hide

vector with indices of aliquots that should be removed from the plot.

exterr

propagate the external (decay constant and calibration factor) uncertainties?

i2i

‘isochron to intercept’: calculates the initial (aka ‘inherited’, ‘excess’, or ‘common’) ^{40}Ar/^{36}Ar ratio from an isochron fit. Setting i2i to FALSE uses the default values stored in settings('iratio',...)

Details

IsoplotR defines the ‘plateau age’ as the weighted mean age of the longest sequence (in terms of cumulative ^{39}Ar content) of consecutive heating steps that pass the modified Chauvenet criterion (see weightedmean). Note that this definition is different (and simpler) than the one used by Isoplot (Ludwig, 2003). However, it is important to mention that all definitions of an age plateau are heuristic by nature and should not be used for quantitative inference.

Value

If plateau=TRUE, returns a list with the following items:

mean

a 3-element vector with:

x: the plateau mean

s[x]: the estimated standard deviation of x

ci[x]: the width of a 100(1-α)% confidence interval of t

disp

a 3-element vector with:

w: the overdispersion, i.e. the standard deviation of the Normal distribution that is assumed to describe the true ages.

ll: the width of the lower half of a 100(1-α)% confidence interval for the overdispersion

ul: the width of the upper half of a 100(1-α)% confidence interval for the overdispersion

df

the degrees of freedom for the weighted mean plateau fit

mswd

the mean square of the weighted deviates of the plateau

p.value

the p-value of a Chi-square test with df=n-2 degrees of freedom, where n is the number of steps in the plateau and 2 degrees of freedom have been removed to estimate the mean and the dispersion.

fract

the fraction of ^{39}Ar contained in the plateau

plotpar

plot parameters for the weighted mean (see weightedmean), which are not used in the age spectrum

i

indices of the steps that are retained for the plateau age calculation

See Also

weightedmean

Examples

1
2

IsoplotR documentation built on Dec. 9, 2018, 1:04 a.m.