NEWS.md
getpaired()/getstd()/globalstd() erroring with “$ operator is invalid for atomic vectors” on datasets where no retention-time cluster contains a single (solo) peak: an empty solo extraction dropped the list element and list$solo$rtg then partial-matched the soloindex logical vector. solo is now always kept as a (possibly empty) data frame.db argument to getchainseq() and data(pmdchain), a built-in database of curated, named multi-step reaction chains (e.g. desaturation-elongation, hydroxylation-glucuronidation, sequential oxidation) for known-pathway screening. getchainseq(list, db = pmdchain) runs every database chain against a feature list and returns chainsearch (which chains are present, each with the best-correlated matched path) plus a keggcount specificity score – the number of compound-connected KEGG reaction paths sharing the chain’s paired-mass-distance signature (lower = rarer = more specific; xenobiotic phase-II conjugations are most specific, common edits like methylation/phosphorylation most degenerate). Each step PMD in data(pmdchain) is the exact monoisotopic mass of its CHNOPS element delta (built by data-raw/pmdchain.R). A match is a relational annotation (co-varying, chromatographically resolved features related by the chain’s mass differences), not compound identification or proof that the pathway occurred.runPMDnet() Shiny app (between-compound / reaction networks) gained a sample-PCA panel for diagnosing whether the dominant variance is biology or a nuisance confound; a targeted reaction-chain search powered by the getchainseq() multi-step PMD grammar (ordered PMDs with wildcards and quantifiers), with downloadable matched chains; and a no-expertise “Screen for known reaction chains” panel that runs getchainseq(db = pmdchain) at one click and returns the named chains present ranked by keggcount specificity (no grammar needed, downloadable).runPMD() Shiny app gained the same no-expertise “Screen for known reaction chains” panel (getchainseq(db = pmdchain), one click, ranked by keggcount specificity, downloadable).pmd.Rmd) so runPMD() launches on case-sensitive filesystems.data.table yielding significant performance gains.getchainseq() to parse structural/reaction-directed MS sequences combining both defined PMDs and wildcards/regex-style quantifiers.gethomolog() specifically designed to extract homologous series in mass spectrometry data based on repeating PMD units.getcorcluster and getcluster to getpseudospectrum and getcorpseudospectrum for pseudo spectrum extractiongetpmddf to show pmd details with given m/z only data and m/z group information(optional, segmentation m/z group for spatial reactomics analysis)getchain
getrda and getpaired to retain reaction related PMDsgetrda for pmd sets and mass defect tablegetchain corcutoff to 0.6