set_Maaslin2: set_Maaslin2

In mcalgaro93/benchdamic: Benchmark of differential abundance methods on microbiome data

set_Maaslin2

Description

Set the parameters for Maaslin2 differential abundance detection method.

Usage

set_Maaslin2(
  assay_name = "counts",
  normalization = c("TSS", "CLR", "CSS", "NONE", "TMM"),
  transform = c("LOG", "LOGIT", "AST", "NONE"),
  analysis_method = c("LM", "CPLM", "ZICP", "NEGBIN", "ZINB"),
  correction = "BH",
  random_effects = NULL,
  fixed_effects = NULL,
  contrast = NULL,
  reference = NULL,
  expand = TRUE
)

Arguments

assay_name	the name of the assay to extract from the TreeSummarizedExperiment object (default assayName = "counts"). Not used if the input object is a phyloseq.
normalization	The normalization method to apply.
transform	The transform to apply.
analysis_method	The analysis method to apply.
correction	The correction method for computing the q-value.
random_effects	The random effects for the model, comma-delimited for multiple effects.
fixed_effects	The fixed effects for the model, comma-delimited for multiple effects.
contrast	character vector with exactly, three elements: a string indicating the name of factor whose levels are the conditions to be compared, the name of the level of interest, and the name of the other level.
reference	The factor to use as a reference for a variable with more than two levels provided as a string of 'variable,reference' semi-colon delimited for multiple variables.
expand	logical, if TRUE create all combinations of input parameters (default expand = TRUE).

Value

A named list containing the set of parameters for DA_Maaslin2 method.

See Also

DA_Maaslin2

Examples

# Set some basic combinations of parameters for Maaslin2
base_Maaslin2 <- set_Maaslin2(normalization = "TSS", transform = "LOG",
    analysis_method = "LM", fixed_effects = "group",
    contrast = c("group", "B", "A"))
many_Maaslin2 <- set_Maaslin2(normalization = c("TSS", "CLR", "CSS", "TMM", 
    "NONE"), transform = c("LOG", "NONE"),
    analysis_method = c("LM", "NEGBIN"), fixed_effects = "group",
    contrast = c("group", "B", "A"))

mcalgaro93/benchdamic documentation built on Nov. 23, 2024, 5:06 p.m.