Proteomics Data Analysis Functions
Description
Data analysis of proteomics experiments by mass spectrometry is supported by this collection of functions mostly dedicated to the analysis of (bottom-up) quantitative (XIC) data.
Fasta-formatted proteomes (eg from UniProt Consortium ) can be read with automatic parsing and multiple annotation types (like species origin, abbreviated gene names, etc) extracted.
Quantitative proteomics measurements frequently contain multiple NA values, due to physical absence of given peptides in some samples, limitations in sensitivity or other reasons.
The functions provided here help to inspect graphically the data to investigate the nature of NA-values via their respective replicate measurements and to help/confirm the choice of NA-replacement by low random values.
Dedicated filtering and statistical testing using the framework of package 'limma' can be run, enhanced by multiple rounds of NA-replacements to provide robustness towards rare stochastic events.
Multi-species samples, as frequently used in benchmark-tests (eg Navarro et al 2016 , Ramus et al 2016 ), can be run with special options separating the data into sub-groups during normalization and testing.
As example the data-set from Ramus et al 2016 ) is provided quantified by MaxQuant (Tyanova et al 2016 ), ProteomeDiscoverer,
OpenMS () and Proline (Bouyssie et al 2020 ).
Subsequently, ROC curves (Hand and Till 2001 ) can be constructed to compare multiple analysis approaches.