ora_est {mseapca} | R Documentation |
ORA adjusting for undetected metabolites
Description
This function performs metabolite set enrichment analysis using over-representation analysis (ORA), incorporating point estimates to adjust for potentially undetected metabolites. It supports three estimation methods: naive, weighted, and shrinkage-based adjustment.
Usage
ora_est(SIG, DET, M, method = "naive", lambda = 5)
Arguments
SIG |
Character vector of significant metabolites |
DET |
Character vector of detected metabolites (background set) |
M |
A named list, where each element is a metabolite set (e.g., pathway) containing character vectors of metabolites |
method |
A character string specifying the estimation method to use. One of |
lambda |
A numeric value used in the |
Details
This function estimates the impact of undetected metabolites on enrichment results. It builds upon the ORA results from detected metabolites, then adjusts the contingency tables by estimating how many undetected metabolites might be significant, based on a specified method.
Value
A list containing:
-
Result of MSEA (ORA with adjustment)
: A matrix with raw p-values and adjusted q-values (BH correction) for each metabolite set. -
significant metabolites
: A list of significant metabolites overlapping with each metabolite set. -
Contingency tables
: A list of 2×2 contingency tables used in Fisher's exact tests.
Author(s)
Hiroyuki Yamamoto
References
Yamamoto H, Fujimori T, Sato H, Ishikawa G, Kami K, Ohashi Y, Statistical hypothesis testing of factor loading in principal component analysis and its application to metabolite set enrichment analysis. BMC Bioinformatics, (2014) 15(1):51.
Yamamoto H. Probabilistic Over-Representation Analysis for Metabolite Set Enrichment Analysis Considering Undetected Metabolites", Jxiv, (2024).
Examples
# Example: Metabolome data
data(fasting_mseapca)
SIG <- fasting_mseapca$SIG
DET <- fasting_mseapca$DET
M <- fasting_mseapca$pathway
# Perform ORA using detected metabolites only
B <- ora_est(SIG, DET, M)
B$`Result of MSEA (ORA with adjustment)`