Skip to content

An R interface to the omixer-rpm for turning metagenomic functional profiles into pathway/module profiles

License

Notifications You must be signed in to change notification settings

omixer/omixer-rpmR

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

27 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

omixer-rpmR

An R interface to omixer-rpm, the tool for metabolic module profiling of microbiome samples

Dependencies

R and Java8 (Docker users, please make sure Java8 is part of the R image)

Installation

Binaries for Linux only (tested on 16.04.1-Ubuntu)

Download the latest binary (omixeRpm_x.y.z.tar.gz) from the release page, then install it as follows after replacing x.y.z by the correct version

R CMD INSTALL omixeRpm_x.y.z.tar.gz

From source

Download the latest source from the release page, then install it as follows after replacing x.y.z by the correct version

R CMD INSTALL omixer-rpmR-x.y.z.tar.gz

Docker Shiny [experimental]

Start the docker container

docker run --rm -ti -p 3838:3838 -v $PWD:/workspace omixer/shinyrpm:0.1

then open the browser at http://localhost:3838/sample-apps/rpm/

Usage

Mapping example

Download the example matrix.tsv (raw link) form the test directory.

library(omixerRpm)
# read a functional profile matrix into R or create it inside R. Please note that row.names should not be used while reading the matrix. 
dat <- read.table("matrix.tsv", header=T, sep="\t")
# Run the module mapping on the loaded table.
mods <- rpm(dat, minimum.coverage=0.3, annotation = 1)

# alternatively run the mapping without loading the table into R.
mods <- rpm("matrix.tsv", minimum.coverage=0.3, annotation = 1)

# Load the default mapping database
db <- loadDefaultDB()
# get the name of the first predicted module
getNames(db, mods@annotation[1,])

# get the abundance|coverage as a data.frame with module id and description
coverage <- asDataFrame(mods, "coverage")
Using an alternative database, several options are available
  1. load one of the bundled databases. Type listDB() to check the list of available databases
db <- loadDB("GBMs.v1.0")
  1. load an external database. Please refer to this module.list and module.names for examples
db <- ModuleDB(directory = "/path/to/moduledb/", modules = "module.list", module.names.file="module.names")

Bundled databases

  1. Gut Brain Modules, Valles-Colomer et al. 2019, The neuroactive potential of the human gut microbiota in quality of life and depression, Nature Microbiology 2019.
  2. Gut Metabolic Modules, Vieira-Silva et al. 2016, Species-function relationships shape ecological properties of the human gut microbiome, Nature Microbiology 2016.

Citing omixer-rpmR

omixer-rpmR was developed as part of GOmixer. If you use omixer-rpmR in your work please cite:

Youssef Darzi, Gwen Falony, Sara Silva, Jeroen Raes. Towards biome-specific analysis of meta-omics data, The ISME journal, 2015.

License

GNU General Public License v3.0. The bundled omixer-rpm.jar is licensed under an Academic Non-commercial Software License Agreement

About

An R interface to the omixer-rpm for turning metagenomic functional profiles into pathway/module profiles

Topics

Resources

License

Stars

Watchers

Forks

Packages

No packages published