mantel_test {metan}

Mantel test

Description

Performs a Mantel test between two correlation/distance matrices. The function calculates the correlation between two matrices, the Z-score that is is the sum of the products of the corresponding elements of the matrices and a two-tailed p-value (null hypothesis: \[r = 0\]).

Usage

mantel_test(mat1, mat2, nboot = 1000, plot = FALSE)

Arguments

mat1, mat2

A correlation matrix or an object of class dist.

nboot

The number of permutations to be used. Defaults to 1000.

plot

if plot = TRUE, plots the density estimate of the permutation distribution along with the observed Z-score as a vertical line.

Value

• mantel_r The correlation between the two matrices.

• z_score The Z-score.

• p-value The quantile of the observed Z-score. in the permutation distribution.

Author(s)

Tiago Olivoto [e-mail]

See Also

pairs_mantel()

Examples

library(metan)
# Test if the correlation of traits (data_ge2 dataset)
# changes between A1 and A2 levels of factor ENV
A1 <- corr_coef(data_ge2 %>% subset(ENV == "A1"))[["cor"]]
A2 <- corr_coef(data_ge2 %>% subset(ENV == "A2"))[["cor"]]
mantel_test(A1, A2, plot = TRUE)