run_metabolomeDA: Perform differential analysis on metabolomic data

In HuaZou/MicrobiomeAnalysis: Data analysis toolkits in metagenomics

Perform differential analysis on metabolomic data

Description

Differential expression analysis contains foldchange analysis, VIP and t-test on metabolomic data.

Usage

run_metabolomeDA(
    object_raw = NULL,
    object_norm = NULL,
    variable,
    variable_name,
    DA_method = c("all", "fc", "vip", "t"),
    method_VIP = c("PLS", "OPLS"),
    cutoff_prev = 0.1,
    p_adjust = c("none", "fdr", "bonferroni", "holm",
                 "hochberg", "hommel", "BH", "BY"))

Arguments

object_raw	(Optional). a SummarizedExperiment::SummarizedExperiment object with raw intensity values (default: NULL).
object_norm	(Optional). a SummarizedExperiment::SummarizedExperiment object with normalized intensity values (default: NULL).
variable	(Required). character. grouping variable for test.
variable_name	(Required). two characters. variable' names.
DA_method	(Optional). character. method for differential analysis (default: "all"). options include: "all": three methods. "fc": foldchange analysis. "vip": VIP analysis. "t": t-test analysis.
method_VIP	(Optional). character. method for VIP (default: "PLS"), options include: "PLS": Partial Least Squares Discriminant Analysis. "OPLS": Orthogonal Partial Least Square Discriminant Analysis.
cutoff_prev	(Optional). Numeric. the Prevalence threshold (default: 0.1).
p_adjust	(Optional). character. method to adjust p-values by. Options include "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". See stats::p.adjust() for more details (default: "BH").

Details

To identify the potential significant metabolites is important in metabolomics. Here, we use three methods to obtain the results:

1. FoldChange: the raw intensity values;

2. VIP (Variable important in projection) by PLS-DA: the normalized intensity values;

3. T-test: the normalized intensity values.

Combining all the results and providing more choice for users to determine the differential metabolites (Recommend: Log2FoldChange and AdjustedPvalue).

Value

A data frame of the differential results.

Author(s)

Created by Hua Zou (12/2/2023 Shenzhen China)

Examples

data("Zeybel_2022_protein")
Zeybel_2022_protein_imput <- impute_abundance(
      Zeybel_2022_protein,
      group = "LiverFatClass",
      method = "knn")
Zeybel_2022_protein_norm <- scale_variables(
      Zeybel_2022_protein_imput,
      method == "zscore")
DA_results <- run_metabolomeDA(
  object_raw = Zeybel_2022_protein,
  object_norm = Zeybel_2022_protein_norm,
  variable = "LiverFatClass",
  variable_name = c("None", "Severe"))

## Not run: 


## End(Not run)

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