mantel_test function

Mantel test \loadmathjax

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: \mjsdeqn r = 0).

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.

Returns

• 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.

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)

See Also

pairs_mantel()

Author(s)

Tiago Olivoto tiagoolivoto@gmail.com