createDatabase: Create database

In msPurity: Automated Evaluation of Precursor Ion Purity for Mass Spectrometry Based Fragmentation in Metabolomics

Description

** General **

Create and SQLite database of an LC-MS(/MS) experiment (replaces the create_database function).

Schema details can be found here.

Example LC-MS/MS processing workflow

• Purity assessments

  • (mzML files) -> purityA -> (pa)

• XCMS processing

  • (mzML files) -> xcms.xcmsSet -> xcms.merge -> xcms.group -> xcms.retcor -> xcms.group -> (xset)

• Fragmentation processing

  • (xset, pa) -> frag4feature -> filterFragSpectra -> averageAllFragSpectra -> createDatabase -> spectralMatching -> (sqlite spectral database)

Usage

createDatabase(
  pa,
  xset,
  xsa = NULL,
  outDir = ".",
  grpPeaklist = NA,
  dbName = NA,
  metadata = NA
)

Arguments

pa	purityA object; Needs to be the same used for frag4feature function
xset	xcms object; Needs to be the same used for frag4feature function (this will be ignored when using xsa parameter)
xsa	CAMERA object (optional); if CAMERA object is used, we ignore the xset parameter input and obtain all information from the xset object nested with the CAMERA xsa object. Adduct and isotope information will be included into the database when using this parameter. The underlying xset object must be the one used for the frag4feature function
outDir	character; Out directory for the SQLite result database
grpPeaklist	dataframe (optional); Can use any peak dataframe. Still needs to be derived from the xset object though
dbName	character (optional); Name of the result database
metadata	list; A list of metadata to add to the s_peak_meta table

Value

path to SQLite database and database name

Examples

#msmsPths <- list.files(system.file("extdata", "lcms", "mzML",
# package="msPurityData"), full.names = TRUE, pattern = "MSMS")
#xset <- xcms::xcmsSet(msmsPths)
#xset <- xcms::group(xset)
#xset <- xcms::retcor(xset)
#xset <- xcms::group(xset)

#pa  <- purityA(msmsPths)
#pa <- frag4feature(pa, xset)
#pa <- filterFragSpectra(pa, allfrag=TRUE)
#pa <- averageAllFragSpectra(pa)
#dbPth <- createDatabase(pa, xset, metadata=list('polarity'='positive',
#'instrument'='Q-Exactive'))

# Run from previously generated data:
pa <- readRDS(system.file("extdata", "tests", "purityA",
          "9_averageAllFragSpectra_with_filter_pa.rds", package="msPurity"))
xset <- readRDS(system.file("extdata","tests", "xcms",
                "msms_only_xset.rds", package="msPurity"))
msmsPths <- list.files(system.file("extdata", "lcms", "mzML",
               package="msPurityData"), full.names = TRUE, pattern = "MSMS")
pa@fileList[1] <- msmsPths[basename(msmsPths)=="LCMSMS_1.mzML"]
pa@fileList[2] <- msmsPths[basename(msmsPths)=="LCMSMS_2.mzML"]
xset@filepaths[1] <- msmsPths[basename(msmsPths)=="LCMSMS_1.mzML"]
xset@filepaths[2] <- msmsPths[basename(msmsPths)=="LCMSMS_2.mzML"]
td <- tempdir()
db_pth = createDatabase(pa, xset, outDir = td)

msPurity documentation built on Jan. 14, 2021, 2:44 a.m.