getPERMANOVA: Calculate PERMANOVA (Permutational Multivariate Analysis of...
In FelixErnst/mia: Microbiome analysis

getPERMANOVA

R Documentation

Calculate PERMANOVA (Permutational Multivariate Analysis of Variance)

Description

These functions perform PERMANOVA to assess the significance of group differences based on a specified dissimilarity matrix. The results can be returned directly or added to metadata in an object of class TreeSummarizedExperiment.

Usage

getPERMANOVA(x, ...)

addPERMANOVA(x, ...)

## S4 method for signature 'SingleCellExperiment'
getPERMANOVA(x, ...)

## S4 method for signature 'SummarizedExperiment'
getPERMANOVA(x, assay.type = "counts", formula = NULL, col.var = NULL, ...)

## S4 method for signature 'ANY'
getPERMANOVA(x, formula, data, method = "bray", test.homogeneity = TRUE, ...)

## S4 method for signature 'SummarizedExperiment'
addPERMANOVA(x, name = "permanova", ...)

Arguments

`x`	`TreeSummarizedExperiment`.
`...`	additional arguments passed to `vegan::adonis2`. `by`: `Character scalar`. Specifies how significance is calculated. (Default: `"margin"`) `na.action`: `function`. Action to take when missing values for any of the variables in `formula` are encountered. (Default: `na.fail`) `full` `Logical scalar`. should all the results from the homogeneity calculations be returned. When `FALSE`, only summary tables are returned. (Default: `FALSE`) `homogeneity.test`: `Character scalar`. Specifies the significance test used to analyse `vegan::betadisper` results. Options include 'permanova' (`vegan::permutest`), 'anova' (`stats::anova`) and 'tukeyhsd' (`stats::TukeyHSD`). (Default: `"permanova"`)
`assay.type`	`Character scalar`. Specifies which assay to use for calculation. (Default: `"counts"`)
`formula`	`formula`. If `x` is a `SummarizedExperiment` a formula can be supplied. Based on the right-hand side of the given formula `colData` is subset to `col.var`. `col.var` and `formula` can be missing, which turns the CCA analysis into a CA analysis and dbRDA into PCoA/MDS.
`col.var`	`Character scalar`. When `x` is a `SummarizedExperiment`,`col.var` can be used to specify variables from `colData`.
`data`	`data.frame` or coarcible to one. The covariance table including covariates defined by `formula`.
`method`	`Character scalar`. A dissimilarity metric used in PERMANOVA and group dispersion calculation. (Default: `"bray"`)
`test.homogeneity`	`Logical scalar`. Should the homogeneity of group dispersions be evaluated? (Default: `TRUE`)
`name`	`Character scalar`. A name for the results that will be stored to metadata. (Default: `"permanova"`)

Details

PERMANOVA is a non-parametric method used to test whether the centroids of different groups (as defined by the formula or covariates) are significantly different in terms of multivariate space.

PERMANOVA relies on the assumption of group homogeneity, meaning the groups should be distinct and have similar variances within each group. This assumption is essential as PERMANOVA is sensitive to differences in within-group dispersion, which can otherwise confound results. This is why the functions return homogeneity test results by default.

The functions utilize vegan::adonis2 to compute PERMANOVA. For homogeneity testing, vegan::betadisper along with vegan::permutest are utilized by default, which allow testing for equal dispersion across groups and validate the homogeneity assumption.

PERMANOVA and distance-based redundancy analysis (dbRDA) are closely related methods for analyzing multivariate data. PERMANOVA is non-parametric, making fewer assumptions about the data. In contrast, dbRDA assumes a linear relationship when constructing the ordination space, although it also employs a PERMANOVA-like approach to test the significance of predictors within this space. dbRDA offers a broader scope overall, as it provides visualization of the constrained ordination, which can reveal patterns and relationships. However, when the underlying data structure is non-linear, the results from the two methods can differ significantly due to dbRDA's reliance on linear assumptions.

Value

getPERMANOVA returns the PERMANOVA results or a list containing the PERMANOVA results and homogeneity test results if test.homogeneity = TRUE. addPERMANOVA adds these results to metadata of x.

Examples


data(GlobalPatterns)
tse <- GlobalPatterns

# Apply relative transformation
tse <- transformAssay(tse, method = "relabundance")
# Perform PERMANOVA
tse <- addPERMANOVA(
    tse,
    assay.type = "relabundance",
    method = "bray",
    formula = x ~ SampleType,
    permutations = 99
    )
# The results are stored to metadata
metadata(tse)[["permanova"]]

# Calculate dbRDA
rda_res <- getRDA(
    tse, assay.type = "relabundance", method = "bray",
    formula = x ~ SampleType, permutations = 99)
# Significance results are similar to PERMANOVA
attr(rda_res, "significance")

FelixErnst/mia documentation built on Nov. 18, 2024, 5:02 a.m.