Wasserstein.diss

Wasserstein.diss.default

Wasserstein.diss.distributional

Wasserstein.diss.varietal

Returns the Wasserstein distance between two samples

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

Wasserstein.diss function

<dl><dt>x</dt>
<dd>the first sample as a vector</dd>
<dt>...</dt>
<dd>optional arguments to the
<code>transport::wasserstein()</code> or
<code>T4transport::wasserstein()</code> functions. Warning: the
latter function is very slow.</dd>
<dt>y</dt>
<dd>the second sample as a vector</dd>
<dt>log</dt>
<dd>logical. Take the lograthm of the data before
calculating the distances?</dd>
<dt>package</dt>
<dd>the name of the package that provides the 2D
Wasserstein distance. Currently, this can be either
<code>'transport'</code> or <code>T4transport</code>.</dd>
<dt>verbose</dt>
<dd>logical. If <code>TRUE</code>, gives progress updates
during the construction of the dissimilarity matrix.</dd></dl>

Arguments

The default S3 method was written by Pieter Vermeesch,
 using modified code from Dominic Schuhmacher's <code>transport</code>
 package (<code>transport1d</code> function), as implemented in
 <code>IsoplotR</code>.

Author

Wasserstein distance — Wasserstein.diss

<dl>

<dt>x</dt>
<dd>the first sample as a vector</dd>


<dt>...</dt>
<dd>optional arguments to the
<code>transport::wasserstein()</code> or
<code>T4transport::wasserstein()</code> functions. Warning: the
latter function is very slow.</dd>


<dt>y</dt>
<dd>the second sample as a vector</dd>


<dt>log</dt>
<dd>logical. Take the lograthm of the data before
calculating the distances?</dd>


<dt>package</dt>
<dd>the name of the package that provides the 2D
Wasserstein distance. Currently, this can be either
<code>'transport'</code> or <code>T4transport</code>.</dd>


<dt>verbose</dt>
<dd>logical. If <code>TRUE</code>, gives progress updates
during the construction of the dissimilarity matrix.</dd>

</dl>

Wasserstein.diss: Wasserstein distance

Description

Usage

Value

Arguments

Author

Examples