simulate_many_poly_spectra: simulate_many_poly_spectra

In MGMS2: 'MGMS2' for Polymicrobial Samples

Description

The function creates simulated mass spectra in pdf file and returns simulated mass spectra (m/z and intensity values of peaks).

Usage

simulate_many_poly_spectra(
  mono.info,
  nsim = 10000,
  file = NULL,
  mixture.ratio,
  mixture.missing.prob.peak = 0.05,
  noise.peak.ratio = 0.05,
  snr.basepeak = 500,
  noise.cv = 0.25,
  mz.range = c(1000, 2200),
  mz.tol = 0.5
)

Arguments

mono.info	A list output of gather_summary or gather_summary_file.
nsim	The number of simulated spectra. (Default: 10000)
file	An output file name. (By default, file=NULL. No pdf file will be generated.)
mixture.ratio	A list of bacterial mixture ratios for given bacterial species in sim.template.
mixture.missing.prob.peak	A real value. The missing probability caused by mixing multiple bacteria species. (Default: 0.05)
noise.peak.ratio	A ratio between the numbers of noise and signal peaks. (Default: 0.05)
snr.basepeak	A (base peak) signal to noise ratio. (Default: 5000)
noise.cv	A coefficient of variation of noise peaks. (Default: 0.25)
mz.range	A range of m/z values. (Default: c(1000,2200))
mz.tol	m/z tolerance. (Default: 0.5)

Value

A list of data frames. A list of simulated mass spectra (data frames) that contains m/z values of peaks, normalized intensities of peaks, species names, and strain names. This function also creates pdf files which contain simulated spectra.

Examples

spectra.processed.A <- process_monospectra(
   file=system.file("extdata", "listA.txt", package="MGMS2"),
   mass.range=c(1000,2200))
spectra.processed.B <- process_monospectra(
   file=system.file("extdata", "listB.txt", package="MGMS2"),
   mass.range=c(1000,2200))
spectra.processed.C <- process_monospectra(
   file=system.file("extdata", "listC.txt", package="MGMS2"),
   mass.range=c(1000,2200))
spectra.mono.summary.A <- summarize_monospectra(
   processed.obj=spectra.processed.A,
   species='A', directory=tempdir())
spectra.mono.summary.B <- summarize_monospectra(
   processed.obj=spectra.processed.B,
   species='B', directory=tempdir())
spectra.mono.summary.C <- summarize_monospectra(
   processed.obj=spectra.processed.C,
   species='C', directory=tempdir())
mono.info=gather_summary(c(spectra.mono.summary.A, spectra.mono.summary.B, spectra.mono.summary.C))
mixture.ratio <- list()
mixture.ratio['A']=1
mixture.ratio['B']=0.5
mixture.ratio['C']=0
insilico.spectra <- simulate_many_poly_spectra(mono.info, mixture.ratio=mixture.ratio, nsim=10)

MGMS2 documentation built on April 24, 2020, 9:06 a.m.