Abstract:Chemical analyses of powdered rocks by different laboratories often yield varying results, requiring estimation of the rock's true composition and associated uncertainty. Challenges arise from the peculiar nature of geochemical data. Traditionally, major and trace elements have been measured using different methods, resulting in chemical analyses where the sum of the parts fluctuates around 1 rather than precisely totaling 1. Additionally, all chemical analyses contain an undisclosed mass fraction representing undetected chemical elements. Because of this undisclosed and unknown mass fraction, geochemical data represent a particular kind of compositional data in which closure to unity is not guaranteed. We argue that chemical analyses exist in the hypercube while being sampled from a true composition residing in the simplex. Therefore, we propose an algorithm that generates random chemical analyses by simulating the data acquisition protocol in geochemistry. Using the algorithm's output, we measure the bias and mean squared error (MSE) of various estimators of the true mean composition. Additionally, we explore the impact of missing values on estimator performance. Our findings reveal that the optimized binary log-ratio mean, a new estimator, exhibits the lowest MSE and bias. It performs well even with up to 70% missing values, in contrast to other classical estimators such as the arithmetic mean or the geometric mean. Applying our approach to the GeoPT database, which contains replicate analyses of igneous rocks from numerous geochemical laboratories, we introduce an outlier detection technique based on the Mahalanobis distance between a laboratory's logit coordinates and the optimized mean estimate. This enables a probabilistic ranking of laboratories based on the atypicality of their performance. Finally, we offer an accessible R implementation of our findings through the GitHub repository linked to this paper [subject classification numbers: 10 (compositions) 85 (statistics)].

Error Propagation in Isometric Log-ratio Coordinates for Compositional Data: Theoretical and Practical Considerations

Spatial Hyperspheric Models for Compositional Data

Units Recovery Methods in Compositional Data Analysis

Classical and Robust Regression Analysis with Compositional Data

Three approaches to supervised learning for compositional data with pairwise logratios

High-dimensional log contrast models with measurement errors

Estimating Rock Composition from Replicate Geochemical Analyses: Theory and Application to Magmatic Rocks of the GeoPT Database

How to choose a proper representation of compositional data for mineral exploration?

Principal balances of compositional data for regression and classification using partial least squares

GeoCoDA: Recognizing and Validating Structural Processes in Geochemical Data. A Workflow on Compositional Data Analysis in Lithogeochemistry

On the distribution of isometric log-ratio transformations under extra-multinomial count data

Modeling the Probability Distribution of Positional Errors Incurred by Residential Address Geocoding.

Regression Models for Compositional Data: General Log-Contrast Formulations, Proximal Optimization, and Microbiome Data Applications

A comprehensive workflow for compositional data analysis in archaeometry, with code in R

Amalgamations are valid in compositional data analysis, can be used in agglomerative clustering, and their logratios have an inverse transformation

A novel, divergence based, regression for compositional data

Spatial Autoregressive Model on a Dirichlet Distribution

Likelihood analysis for a class of spatial geostatistical compositional models

Identifying Important Pairwise Logratios in Compositional Data with Sparse Principal Component Analysis

Exploring geochemical data using compositional techniques: A practical guide

Perturbation-based Effect Measures for Compositional Data