| Title: |
Omnivariant generalized least squares regression: Theory, geochronological applications, and making the case for reconciled Δ$_{47}$ calibrations |
| Authors: |
Daëron, Mathieu; Vermeesch, Pieter |
| Contributors: |
Paléocéanographie (PALEOCEAN); Laboratoire des Sciences du Climat et de l'Environnement Gif-sur-Yvette (LSCE); Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Direction de Recherche Fondamentale (CEA) (DRF (CEA)); Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Direction de Recherche Fondamentale (CEA) (DRF (CEA)); Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA); University College London UCL (UCL) |
| Source: |
ISSN: 0009-2541 ; Chemical Geology ; https://hal.science/hal-04372546 ; Chemical Geology, 2024, 647, pp.121881. ⟨10.1016/j.chemgeo.2023.121881⟩. |
| Publisher Information: |
CCSD; Elsevier |
| Publication Year: |
2024 |
| Collection: |
Université de Versailles Saint-Quentin-en-Yvelines: HAL-UVSQ |
| Subject Terms: |
[SDU.STU.GC]Sciences of the Universe [physics]/Earth Sciences/Geochemistry; [STAT.AP]Statistics [stat]/Applications [stat.AP]; [STAT.TH]Statistics [stat]/Statistics Theory [stat.TH] |
| Description: |
International audience ; Least-squares regression methods are mathematically powerful, conceptually and computationally simple, and widely used in many fields. However, none of the commonly-used flavors of least-squares regression, such as York regression or Generalized Least Squares (GLS), take into account the full set of covariances between all observed (x, y) values. Here we describe the Omnivariant Generalized Least Squares (OGLS) method to fit a model of the form y = f(x), accounting for the full error correlation structure of the (x, y) data, based on a first-order linear propagation of the uncertainties in all variables into errors in y residuals, followed by minimizing the vector of y residuals with respect to the Mahalanobis norm defined by its covariance matrix. This approach may be described as a generalization of both York regression and GLS. It is mathematically exact for straight-line fits, and is also suitable for many non-linear models. Here we describe the principles of OGLS regression and discuss its properties, caveats, and practical use, and provide two consistent open-source implementations in Python and R. To illustrate how various fields of geochronology and stable-isotope geochemistry may benefit from this new method, we discuss how OGLS may specifically apply to $^{40}$Ar/$^{39}$Ar dating and how it provides robust mathematical evidence that Δ$_{47}$ carbonate calibrations in the recently defined I-CDES metrological scale are statistically indistinguishable, effectively solving long-standing methodological discrepancies. |
| Document Type: |
article in journal/newspaper |
| Language: |
English |
| DOI: |
10.1016/j.chemgeo.2023.121881 |
| Availability: |
https://hal.science/hal-04372546; https://hal.science/hal-04372546v1/document; https://hal.science/hal-04372546v1/file/666936.pdf; https://doi.org/10.1016/j.chemgeo.2023.121881 |
| Rights: |
https://about.hal.science/hal-authorisation-v1/ ; info:eu-repo/semantics/OpenAccess |
| Accession Number: |
edsbas.76379C5B |
| Database: |
BASE |