run_aldex: Perform differential analysis using ALDEx2

In HuaZou/MicrobiomeAnalysis: Data analysis toolkits in metagenomics

Perform differential analysis using ALDEx2

Description

The run_aldex function is used to identify the significant features.

Usage

run_aldex(
  ps,
  group,
  taxa_rank = "all",
  transform = c("identity", "log10", "log10p", "SquareRoot", "CubicRoot", "logit"),
  norm = "none",
  norm_para = list(),
  method = c("t.test", "wilcox.test", "kruskal", "glm_anova"),
  p_adjust = c("none", "fdr", "bonferroni", "holm", "hochberg", "hommel", "BH", "BY"),
  pvalue_cutoff = 0.05,
  mc_samples = 128,
  denom = c("all", "iqlr", "zero", "lvha"),
  paired = FALSE
)

Arguments

ps	a phyloseq::phyloseq object
group	character, the variable to set the group
taxa_rank	character to specify taxonomic rank to perform differential analysis on. Should be one of phyloseq::rank_names(phyloseq), or "all" means to summarize the taxa by the top taxa ranks (summarize_taxa(ps, level = rank_names(ps)[1])), or "none" means perform differential analysis on the original taxa (taxa_names(phyloseq), e.g., OTU or ASV).
transform	character, the methods used to transform the microbial abundance. See transform_abundances() for more details. The options include: "identity", return the original data without any transformation (default). "log10", the transformation is log10(object), and if the data contains zeros the transformation is log10(1 + object). "log10p", the transformation is log10(1 + object). "SquareRoot", the transformation is ⁠Square Root⁠. "CubicRoot", the transformation is ⁠Cubic Root⁠. "logit", the transformation is ⁠Zero-inflated Logit Transformation⁠ (Does not work well for microbiome data).
norm	the methods used to normalize the microbial abundance data. See normalize() for more details. Options include: "none": do not normalize. "rarefy": random subsampling counts to the smallest library size in the data set. "TSS": total sum scaling, also referred to as "relative abundance", the abundances were normalized by dividing the corresponding sample library size. "TMM": trimmed mean of m-values. First, a sample is chosen as reference. The scaling factor is then derived using a weighted trimmed mean over the differences of the log-transformed gene-count fold-change between the sample and the reference. "RLE", relative log expression, RLE uses a pseudo-reference calculated using the geometric mean of the gene-specific abundances over all samples. The scaling factors are then calculated as the median of the gene counts ratios between the samples and the reference. "CSS": cumulative sum scaling, calculates scaling factors as the cumulative sum of gene abundances up to a data-derived threshold. "CLR": centered log-ratio normalization. "CPM": pre-sample normalization of the sum of the values to 1e+06.
norm_para	arguments passed to specific normalization methods
method	test method, options include: "t.test" and "wilcox.test" for two groups comparison, "kruskal" and "glm_anova" for multiple groups comparison.
p_adjust	method for multiple test correction, default none, for more details see stats::p.adjust.
pvalue_cutoff	cutoff of p value, default 0.05.
mc_samples	integer, the number of Monte Carlo samples to use for underlying distributions estimation, 128 is usually sufficient.
denom	character string, specifiy which features used to as the denominator for the geometric mean calculation. Options are: "all", with all features. "iqlr", accounts for data with systematic variation and centers the features on the set features that have variance that is between the lower and upper quartile of variance. "zero", a more extreme case where there are many non-zero features in one condition but many zeros in another. In this case the geometric mean of each group is calculated using the set of per-group non-zero features. "lvha", with house keeping features.
paired	logical, whether to perform paired tests, only worked for method "t.test" and "wilcox.test".

Details

The run_aldex function is used to identify the significant features. It can be applied to both phyloseq::phyloseq and Biobase::ExpressionSet object.

Value

a microbiomeMarker object.

Author(s)

Created by Yang Cao; modified by Hua Zou (5/17/2022 Shenzhen China)

References

Fernandes, A.D., Reid, J.N., Macklaim, J.M. et al. Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome 2, 15 (2014).

See Also

ALDEx2::aldex()

Examples

## Not run: 
data(caporaso)
ps <- phyloseq::subset_samples(
    caporaso,
    SampleType %in% c("gut", "tongue", "right palm")
)
run_aldex(ps, group = "SampleType", method = "kruskal")

## End(Not run)

HuaZou/MicrobiomeAnalysis documentation built on May 13, 2024, 11:10 a.m.