DA_corncob {benchdamic} | R Documentation |
Fast run for corncob differential abundance detection method.
DA_corncob( object, pseudo_count = FALSE, formula, phi.formula, formula_null, phi.formula_null, test, boot = FALSE, coefficient = NULL, norm = c("TMM", "TMMwsp", "RLE", "upperquartile", "posupperquartile", "none", "ratio", "poscounts", "iterate", "TSS", "CSSmedian", "CSSdefault"), verbose = TRUE )
object |
phyloseq object. |
pseudo_count |
add 1 to all counts if TRUE (default
|
formula |
an object of class |
phi.formula |
an object of class |
formula_null |
Formula for mean under null, without response |
phi.formula_null |
Formula for overdispersion under null, without response |
test |
Character. Hypothesis testing procedure to use. One of
|
boot |
Boolean. Defaults to |
coefficient |
The coefficient of interest as a single word formed by the variable name and the non reference level. (e.g.: 'ConditionDisease' if the reference level for the variable 'Condition' is 'control'). |
norm |
name of the normalization method used to compute the
normalization factors to use in the differential abundance analysis. If
|
verbose |
an optional logical value. If |
A list object containing the matrix of p-values 'pValMat', the matrix of summary statistics for each tag 'statInfo', and a suggested 'name' of the final object considering the parameters passed to the function.
bbdml
and
differentialTest
for differential abundance and
differential variance evaluation.
set.seed(1) # Create a very simple phyloseq object counts <- matrix(rnbinom(n = 60, size = 3, prob = 0.5), nrow = 10, ncol = 6) metadata <- data.frame("Sample" = c("S1", "S2", "S3", "S4", "S5", "S6"), "group" = as.factor(c("A", "A", "A", "B", "B", "B"))) ps <- phyloseq::phyloseq(phyloseq::otu_table(counts, taxa_are_rows = TRUE), phyloseq::sample_data(metadata)) # No use of scaling factors ps_NF <- norm_edgeR(object = ps, method = "none") # The phyloseq object now contains the scaling factors: scaleFacts <- phyloseq::sample_data(ps_NF)[, "NF.none"] head(scaleFacts) # Differential abundance DA_corncob(object = ps_NF, formula = ~ group, phi.formula = ~ group, formula_null = ~ 1, phi.formula_null = ~ group, coefficient = "groupB", norm = "none", test = "Wald")