read.biom: Extracts counts, metadata, taxonomy, and phylogeny from a...

Description Usage Arguments Value Examples

View source: R/read.biom.r

Description

Extracts counts, metadata, taxonomy, and phylogeny from a biom file.

Usage

1
read.biom(src, tree = "auto", progressbar = NULL)

Arguments

src

Input data as either a file path, URL, or JSON string. read.biom can read BIOM files formatted according to both the version 1.0 (JSON) and 2.1 (HDF5) specifications as well as classical tabular format. URLs must begin with http://, https://, ftp://, or ftps://. JSON files must have { as their first non-whitespace character. Compressed (gzip or bzip2) biom files are also supported.

tree

The default value of auto will read the tree from the biom file specified in src, if present. The value TRUE will do the same, but will generate an error message if a tree is not present. Setting tree=FALSE will return a BIOM object without any tree data. You may also provide a file path, URL, or Newick string to load that tree data into the final BIOM object.

progressbar

An object of class Progress.

Value

A BIOM class object containing the parsed data. This object can be treated as a list with the following named elements:

counts

A numeric slam sparse matrix of observation counts. Taxa (OTUs) as rows and samples as columns.

metadata

A data frame containing any embedded metadata. Row names are sample IDs.

taxonomy

Character matrix of taxonomic names, if given. Row names are taxa (OTU) IDs. Column rows are named Kingdom, Phylum, Class, Order, Family, Genus, Species, and Strain, or TaxLvl.1, TaxLvl.2, ... , TaxLvl.N when more than 8 levels of taxonomy are encoded in the biom file.

phylogeny

An object of class phylo defining the phylogenetic relationships between the taxa. Although the official specification for BIOM only includes phylogenetic trees in BIOM version 2.1, if a BIOM version 1.0 file includes a phylogeny entry with newick data, then it will be loaded here as well. The ape package has additional functions for working with phylo objects.

info

A list of other attributes defined in the BIOM file, such as id, type, format, format_url, generated_by, date, matrix_type, matrix_element_type, and shape

metadata, taxonomy, and phylogeny are optional components of the BIOM file specification and therefore will be empty in the returned object when they are not provided by the BIOM file.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
    library(rbiom)

    infile <- system.file("extdata", "hmp50.biom", package = "rbiom")
    biom <- read.biom(infile)

    summary(biom)

    # Taxa Abundances
    as.matrix(biom$counts[1:4,1:4])

    top5 <- names(head(rev(sort(slam::row_sums(biom$counts))), 5))
    biom$taxonomy[top5,c('Family', 'Genus')]
    as.matrix(biom$counts[top5, 1:6])

    # Metadata
    table(biom$metadata$Sex, biom$metadata$Body.Site)
    sprintf("Mean age: %.1f", mean(biom$metadata$Age))

    # Phylogenetic tree
    tree <- biom$phylogeny
    top5.tree <- rbiom::subtree(tree, top5)
    ape::plot.phylo(top5.tree)

cmmr/rbiom documentation built on Nov. 21, 2018, 10:04 a.m.