IsoplotR (version 0.6)

isochron: Calculate and plot isochrons

Description

Plots cogenetic $^{40}$Ar/$^{39}$Ar data as X-Y scatterplots, fits an isochron curve through them using the yorkfit function, and computes the corresponding isochron age, including decay constant uncertainties.

Usage

isochron(x, ...)
"isochron"(x, xlim = NA, ylim = NA, alpha = 0.05, sigdig = 2, show.numbers = FALSE, ellipse.col = rgb(0, 1, 0, 0.5), line.col = "red", lwd = 2, title = TRUE, ...)
"isochron"(x, xlim = NA, ylim = NA, alpha = 0.05, sigdig = sigdig, show.numbers = FALSE, ellipse.col = rgb(0, 1, 0, 0.5), inverse = TRUE, line.col = "red", lwd = 2, plot = TRUE, ...)

Arguments

x

EITHER a list or a matrix with the following vectors:

X: the x-variable

Y: the y-variable

sX: the standard error of X

sY: the standard error of Y

rXY: the correlation coefficient of X and Y

OR an object of class ArAr

...

optional arguments

xlim

2-element vector with the plot limits of the x-axis

ylim

2-element vector with the plot limits of the y-axis

alpha

confidence cutoff for the error ellipses

sigdig

the number of significant digits of the numerical values reported in the title of the graphical output

show.numbers

boolean flag (TRUE to show grain numbers)

ellipse.col

background colour of the error ellipses

line.col

colour of the isochron line

lwd

line width

title

add a title to the plot?

inverse

if TRUE, plots $^{36}$Ar/$^{40}$Ar vs. $^{39}$Ar/$^{40}$Ar. If FALSE, plots $^{40}$Ar/$^{36}$Ar vs. $^{39}$Ar/$^{36}$Ar.

plot

if FALSE, suppresses the graphical output

Value

if plot=FALSE, returns a list with the following items:a: the intercept of the straight line fit and its standard errorb: the slope of the fit and its standard errory0: the atmospheric $^{40}$Ar/$^{36}$Ar ratio and its standard errorage: the $^{40}$Ar/$^{39}$Ar age and its standard error

Examples

data(examples)
isochron(examples$ArAr)