DA.bay: baySeq

In Russel88/DAtest: Comparing Differential Abundance/Expression Methods

baySeq

Description

Implementation of baySeq for DAtest

Usage

DA.bay(data, predictor, allResults = FALSE, ...)

Arguments

data	Either a matrix with counts/abundances, OR a phyloseq object. If a matrix/data.frame is provided rows should be taxa/genes/proteins and columns samples
predictor	The predictor of interest. Factor, OR if data is a phyloseq object the name of the variable in sample_data(data) in quotation
allResults	If TRUE will return raw results from the getLikelihoods function
...	Additional arguments to the getPriors.NB and getLikelihoods functions

Value

A data.frame with with results.

Examples

# Creating random count_table and predictor
set.seed(4)
mat <- matrix(rnbinom(100, size = 0.1, mu = 500), nrow = 50, ncol = 10)
rownames(mat) <- 1:50
pred <- c(rep("Control", 5), rep("Treatment", 5))

# Running baySeq
res <- DA.bay(data = mat, predictor = pred)

Russel88/DAtest documentation built on March 24, 2022, 3:50 p.m.