Turns a vector of numbers into an object of class <code>KDE</code> using
a combination of the Botev (2010) bandwidth selector and the
Abramson (1982) adaptive kernel bandwidth modifier.

Bundles a number of established statistical methods to facilitate the visual interpretation of large datasets in sedimentary geology. Includes functionality for adaptive kernel density estimation, principal component analysis, correspondence analysis, multidimensional scaling, generalised procrustes analysis and individual differences scaling using a variety of dissimilarity measures. Univariate provenance proxies, such as single-grain ages or (isotopic) compositions are compared with the Kolmogorov-Smirnov, Kuiper, Wasserstein-2 or Sircombe-Hazelton L2 distances. Categorical provenance proxies such as chemical compositions are compared with the Aitchison and Bray-Curtis distances,and count data with the chi-square distance. Varietal data can either be converted to one or more distributional datasets, or directly compared using the multivariate Wasserstein distance. Also included are tools to plot compositional and count data on ternary diagrams and point-counting data on radial plots, to calculate the sample size required for specified levels of statistical precision, and to assess the effects of hydraulic sorting on detrital compositions. Includes an intuitive query-based user interface for users who are not proficient in R.

Pieter Vermeesch

provenance

Statistical Toolbox for Sedimentary Provenance Analysis

KDE function

<dl><dt>x</dt>
<dd>a vector of numbers</dd>
<dt>from</dt>
<dd>minimum age of the time axis. If <code>NULL</code>, this is
set automatically</dd>
<dt>to</dt>
<dd>maximum age of the time axis. If <code>NULL</code>, this is set
automatically</dd>
<dt>bw</dt>
<dd>the bandwidth of the KDE. If NULL, <code>bw</code> will be
calculated automatically using <code>botev()</code></dd>
<dt>adaptive</dt>
<dd>boolean flag controlling if the adaptive KDE
modifier of Abramson (1982) is used</dd>
<dt>log</dt>
<dd>transform the ages to a log scale if <code>TRUE</code></dd>
<dt>n</dt>
<dd>horizontal resolution of the density estimate</dd>
<dt>...</dt>
<dd>optional arguments to be passed on to <code>density</code></dd></dl>

Arguments

Create a kernel density estimate — KDE

<dl>

<dt>x</dt>
<dd>a vector of numbers</dd>


<dt>from</dt>
<dd>minimum age of the time axis. If <code>NULL</code>, this is
set automatically</dd>


<dt>to</dt>
<dd>maximum age of the time axis. If <code>NULL</code>, this is set
automatically</dd>


<dt>bw</dt>
<dd>the bandwidth of the KDE. If NULL, <code>bw</code> will be
calculated automatically using <code>botev()</code></dd>


<dt>adaptive</dt>
<dd>boolean flag controlling if the adaptive KDE
modifier of Abramson (1982) is used</dd>


<dt>log</dt>
<dd>transform the ages to a log scale if <code>TRUE</code></dd>


<dt>n</dt>
<dd>horizontal resolution of the density estimate</dd>


<dt>...</dt>
<dd>optional arguments to be passed on to <code>density</code></dd>

</dl>

KDE: Create a kernel density estimate

Description

Usage

Value

Arguments

See Also

Examples