Masserstein: Linear regression of mass spectra by optimal transport
Autor: | Grzegorz Skoraczynski, Dirk Valkenborg, Michał Aleksander Ciach, Michał Startek, Szymon Majewski, Anna Gambin, Błażej Miasojedow |
---|---|
Rok vydání: | 2021 |
Předmět: |
Chemistry
010401 analytical chemistry Organic Chemistry Python (programming language) 01 natural sciences Regression 0104 chemical sciences Analytical Chemistry Data set ComputingMethodologies_PATTERNRECOGNITION Linear regression Mass spectrum Preprocessor Computational problem Linear combination computer Algorithm Spectroscopy computer.programming_language |
Zdroj: | Rapid Communications in Mass Spectrometry. |
ISSN: | 1097-0231 0951-4198 |
DOI: | 10.1002/rcm.8956 |
Popis: | RATIONALE The linear regression of mass spectra is a computational problem defined as fitting a linear combination of reference spectra to an experimental one. It is typically used to estimate the relative quantities of selected ions. In this work, we study this problem in an abstract setting to develop new approaches applicable to a diverse range of experiments. METHODS To overcome the sensitivity of the ordinary least-squares regression to measurement inaccuracies, we base our methods on a non-conventional spectral dissimilarity measure, known as the Wasserstein or the Earth Mover's distance. This distance is based on the notion of the cost of transporting signal between mass spectra, which renders it naturally robust to measurement inaccuracies in the mass domain. RESULTS Using a data set of 200 mass spectra, we show that our approach is capable of estimating ion proportions accurately without extensive preprocessing of spectra required by other methods. The conclusions are further substantiated using data sets simulated in a way that mimics most of the measurement inaccuracies occurring in real experiments. CONCLUSIONS We have developed a linear regression algorithm based on the notion of the cost of transporting signal between spectra. Our implementation is available in a Python 3 package called masserstein, which is freely available at https://github.com/mciach/masserstein. |
Databáze: | OpenAIRE |
Externí odkaz: |