介绍
Mantel test,由Nathan Mantel在1976年提出,是一种用于检验两个矩阵相关关系的统计方法。它特别适用于分析不同样本的变量在不同类别中的相关性,比如生态学中样本间的两两距离矩阵的相关性分析。与传统的相关系数不同,相关系数主要处理两列数据之间的相关性,而Mantel test能够处理两个矩阵之间的相关性问题。
Mantel test的原假设是两个矩阵之间没有相关关系。检验过程包括将两个矩阵展开为变量两列,计算相关系数(常用Pearson相关系数),然后通过置换测试来评估相关性。具体来说,会进行成千上万次的置换,每次置换后重新计算相关系数,并形成相关系数的分布,最后比较实际观测到的相关系数在该分布中的位置。
在实际应用中,Mantel test可以分析例如微生物群落与植被群落之间的关系,通过比较微生物的Unifrac矩阵和植物的Bray-Curtis距离矩阵来进行相关性分析。此外,Mantel test的结果会受到相关方法的影响,例如使用不同的相关系数计算方法(如Pearson、Spearman等)可能会导致不同的显著性结果。
加载R包
knitr::opts_chunk$set(message = FALSE, warning = FALSE)
library(tidyverse)
library(linkET)
library(vegan)
# rm(list = ls())
options(stringsAsFactors = F)
options(future.globals.maxSize = 10000 * 1024^2)