Paired mass distance (PMD) analysis proposed in Yu, Olkowicz and Pawliszyn (2018) and PMD based reactomics proposed in Yu and Petrick (2020) for gas/liquid chromatography–mass spectrometry (GC/LC-MS) based non-targeted analysis. PMD analysis including GlobalStd algorithm and structure/reaction directed analysis. GlobalStd algorithm could found independent peaks in m/z-retention time profiles based on retention time hierarchical cluster analysis and frequency analysis of paired mass distances within retention time groups. Structure directed analysis could be used to find potential relationship among those independent peaks in different retention time groups based on frequency of paired mass distances. Reactomics analysis could also be performed to build PMD network, assign sources and make biomarker reaction discovery. GUIs for PMD analysis is also included as ‘shiny’ applications.
You can install package from this GitHub repository:
Or find a stable version from CRAN:
PMD analysis paper. This paper proposed structure/reaction directed analysis with PMD.
Reactomics paper. This paper contained the concepts of PMD based reactomics, applications and data mining of reaction database and compounds database.
Slides. This is the slides for ASMS 2020 Reboot and here is the video of presentation. Press “P” and you will see the notes for each slide with details. Another full version of reactomics presentation for one hour presentation could be found here. I will not update the conference presentation while I will add new contents for the full version of reactomics presentation whenever I have new results.
To perform GlobalStd algorithm, use the following code:
library(pmd) data("spmeinvivo") pmd <- getpaired(spmeinvivo) std <- getstd(pmd)
To perform structure/reaction directed analysis, use the following code:
sda <- getsda(std)
To perform GlobalStd algorithem along with structure/reaction directed analysis, use the following code:
result <- globalstd(spmeinvivo)
To use the shiny application within the package, use the following code:
To check the pmd reaction database:
# all reaction data("omics") View(omics) # kegg reaction data("keggrall") View(keggrall) # literature reaction for mass spectrometry data("sda") View(sda)
To check the HMDB pmd database:
To cite related papers: