Statistical Toolbox for Radiometric Geochronology

Plots U-Pb data on Wetherill and Tera-Wasserburg concordia diagrams. Calculates concordia and discordia ages. Performs linear regression of measurements with correlated errors using 'York', 'Titterington', 'Ludwig' and Omnivariant Generalised Least-Squares ('OGLS') approaches. Generates Kernel Density Estimates (KDEs) and Cumulative Age Distributions (CADs). Produces Multidimensional Scaling (MDS) configurations and Shepard plots of multi-sample detrital datasets using the Kolmogorov-Smirnov distance as a dissimilarity measure. Calculates 40Ar/39Ar ages, isochrons, and age spectra. Computes weighted means accounting for overdispersion. Calculates U-Th-He (single grain and central) ages, logratio plots and ternary diagrams. Processes fission track data using the external detector method and LA-ICP-MS, calculates central ages and plots fission track and other data on radial (a.k.a. 'Galbraith') plots. Constructs total Pb-U, Pb-Pb, Th-Pb, K-Ca, Re-Os, Sm-Nd, Lu-Hf, Rb-Sr and 230Th-U isochrons as well as 230Th-U evolution plots.

IsoplotR

IsoplotR is a free and open-source substitute for Kenneth Ludwig's popular Isoplot add-in to Microsoft Excel. This Github repository contains the code for IsoplotR's core commands. A graphical user interface is provided in a separate repository at https://github.com/pvermees/IsoplotRgui.

Installation

You must have R installed on your system (see https://www.r-project.org/). The most recent stable version of IsoplotR is available from CRAN at https://cran.r-project.org/package=IsoplotR and can be installed on your system as follows:

install.packages('IsoplotR')

Alternatively, the current development version of IsoplotR can be installed from Github with the remotes package:

install.packages('remotes')
remotes::install_github('pvermees/IsoplotR')

Example

Once installed, IsoplotR can be loaded into memory by entering the following code at the R command prompt:

library(IsoplotR)

Now we can issue IsoplotR commands to R. For example:

setwd(system.file(package='IsoplotR')) # navigate to the built-in data files
RbSr <- read.data('RbSr1.csv',method='Rb-Sr',format=1)  
isochron(RbSr)

Further information

See https://isoplotr.es.ucl.ac.uk/home/index.html

Vermeesch, P., 2018, IsoplotR: a free and open toolbox for geochronology. Geoscience Frontiers, v.9, p.1479-1493, doi: 10.1016/j.gsf.2018.04.001.

Author

Pieter Vermeesch

License

This project is licensed under the GPL-3 License