Learn R Programming

⚠️There's a newer version (0.8.0) of this package.Take me there.

enviGCMS: GC-MS Data Analysis for Environmental Science

enviGCMS provides functions for GC/LC-MS data analysis for environmental sciences.

Installation

You can either use the stable version of enviGCMS from CRAN,

install.packages("enviGCMS")

or the current development snapshot from this GitHub repository:

devtools::install_github("yufree/enviGCMS")

Usage

  • get the mean and RSD of one sample for 5 technique replicate
# enviGCMS use functions in xcms to import the data, just type the path to your single sample
data1 <- enviGCMS:::getmd(‘sample1-1’)
data2 <- enviGCMS:::getmd(‘sample1-2’)
data3 <- enviGCMS:::getmd(‘sample1-3’)
data4 <- enviGCMS:::getmd(‘sample1-4’)
data5 <- enviGCMS:::getmd(‘sample1-5’)
  • get the mean
data <- (data1+data2+data3+data4+data5)/5
  • get the standard deviation
datasd <- sqrt(((data1-data)^2+(data2-data)^2+(data3-data)^2+(data4-data)^2+(data5-data)^2)/4)
  • get the RSD
databrsd <- datasd/data
  • plot the smooth scatter
plotsms(datarsd)
  • plot the heatmap
plotms(data)
  • plot the mz-rt scatter plot
plotmz(data)
  • plot the boundary model
findline(data)

Detailed usage of those functions in Environmental analysis could be found in this paper and the vignettes in this package.

Copy Link

Version

Install

install.packages('enviGCMS')

Monthly Downloads

360

Version

0.5.7

License

GPL-2

Issues

Pull Requests

Stars

Forks

Maintainer

Miao YU

Last Published

August 22nd, 2019

Functions in enviGCMS (0.5.7)

GetIntegration

GetIntegration was mainly used for get the intergration of certain ion's chromatogram data and plot the data
cbmd

Combine two data with similar retention time while different mass range
findline

find line of the regression model for GC-MS
TBBPA

Demo data for TBBPA metabolism in Pumpkin
batch

Get the MIR and related information from the files
Integration

Just intergrate data according to fixed rt and fixed noise area
Getisotopologues

Get the selected isotopologues at certain MS data
enviGCMS-deprecated

Deprecated functions in package enviGCMS.
findmet

Screen metabolites by Mass Defect
Mode

define the Mode function
findohc

Screen organohalogen compounds by retention time, mass defect analysis and isotope relationship modified by literature report. Also support compounds with [M] and [M+2] ratio cutoff.
getareastd

Get the peak information from SCCPs standards
getarea

Get the peak information from sampels for SCCPs detection
getfilter

Filter the data based on row and column index
getdata2

Get XCMSnExp object in one step from structured folder path for xcms 3.
getmzrt2

Get the mzrt profile and group information for batch correction and plot as a list for xcms 3 object
getdata

Get xcmsset object in one step with optimized methods.
getcsv

Convert an list object to csv file.
getmzrt

Get the mzrt profile and group information as a mzrt list and/or save them as csv or rds for further analysis.
getQCraw

get the data of QC compound for a group of data
getpower

Get the index with power restriction for certain study with BH adjusted p-value and certain power.
getformula

Get chemical formula for mass to charge ratio.
getbgremove

Get the peak list with blank samples' peaks removed
getdoe

Filter the data based on DoE, rsd, intensity
getmd

Import data and return the annotated matrix for GC/LC-MS by m/z range and retention time
getmass

Get the exact mass of the isotopologues from a chemical formula or reaction's isotope patterns with the highest abundances
plothist

plot the density of the GC-MS data with EM algorithm to seperate the data into two log normal distribution.
plotgroup

Plot the response group of GC-MS
getbiotechrep

Get the report for biological replicates.
getfeaturesanova

Get the features from anova, with p value, q value, rsd and power restriction
getmassdefect

Get mass defect with certain scaled factor
getupload2

Get the csv files to be submitted to Metaboanalyst
getoverlappeak

Get the overlap peaks by mass and retention time range
getfeaturest

Get the features from t test, with p value, q value, rsd and power restriction
getoverlaprt

Get the overlap peaks by retention time
getupload

Get the csv files from xcmsset/XCMSnExp/list object
getsccp

Quantitative analysis for short-chain chlorinated paraffins(SCCPs)
getrmd

Get the Relative Mass Defect
getmdh

Get the high order unit based Mass Defect
getmdr

Get the raw Mass Defect
list

Demo data
plotmrc

plot the diff scatter plot for one xcmsset objects with threshold between two groups
plotmr

plot the scatter plot for peaks list with threshold
svaplot

Filter the data with p value and q value and show them
getsim

output the similarity of two dataset
svapca

Principal component analysis(PCA) for SVA corrected data and raw data
plotsms

Plot the intensity distribution of GC-MS
getmr

Get the mzrt profile and group information for batch correction and plot as a list directly from path with default setting
plotsub

Plot the backgrond of data
runMDPlot

Shiny application for interactive mass defect plots analysis
qbatch

Get the MIR from the file
plote

plot EIC and boxplot for all peaks and return diffreport
plotrtms

Plot mass spectrum of certain retention time and return mass spectrum vector (MSP file) for NIST search
plotrsd

plot the rsd influnces of data in different groups
ma

filter data by average moving box
plotden

plot the density for multiple samples
plotmz

plot GC/LC-MS data as scatter plot
plotkms

plot the kendrick mass defect diagram
plotrla

Relative Log Abundance (RLA) plots
gifmr

plot scatter plot for rt-mz profile and output gif file for mutiple groups
getupload3

Get the csv files to be submitted to Metaboanalyst
plotintslope

plot the slope information of intergretion
plotridges

Relative Log Abundance Ridge (RLAR) plots
plott

plot GC-MS data as a heatmap for constant speed of temperature rising
getgrouprep

Get the report for samples with biological and technique replicates in different groups
getimputation

Impute the peaks list data
getmzrtcsv

Covert the peaks list csv file into list
plottic

Plot Total Ion Chromatogram (TIC)
svabatch

Plot the influnces of DoE and Batch effects on each peaks
submd

Get the differences of two GC/LC-MS data
plotpca

plot the PCA for multiple samples
plothm

Plot the heatmap of mzrt profiles
gettechrep

Get the report for technique replicates.
runsccp

Shiny application for Short-Chain Chlorinated Paraffins analysis
plotms

plot GC/LC-MS data as a heatmap with TIC
plotint

plot the information of intergretion
plotmsrt

Plot EIC of certain m/z and return dataframe for intergration
gettimegrouprep

Get the time series or two factor DoE report for samples with biological and technique replicates in different groups
getoverlapmass

Get the overlap peaks by mass range
sccp

Short-Chain Chlorinated Paraffins(SCCPs) peaks infomation for quantitative analysis
svacor

Surrogate variable analysis(SVA) to correct the unknown batch effects
svadata

Filter the data with p value and q value
svaupload

Get the corrected data after SVA for metabolanalyst
writeMSP

Write MSP files for NIST search