PCA、PCoA、NMDS、CCA、RDA的区别

首先，以上分析本质上都属于排序分析（Ordinationanalysis）。排序(ordination)的过程就是在一个可视化的低维空间(通常是二维)重新排列这些样方,使得样方之间的距离最大程度地反映出平面散点图内样方之间的关系信息。1、只使用物种组成数据的排序称作非限制性排序(unconstrainedordination)(1)主成分分析(principalcomponentsanalysis,PCA)(2)对应分析(correspondenceanalysis,CA)(3)去趋势对应分析(Detrendedcorrespondenceanalysis,DCA)(3)主坐标分析(principalcoordinateanalysis,PCoA)(4)非度量多维尺度分析(non-metricmulti-dimensionalscaling,NMDS)2、同时使用物种和环境因子组成数据的排序叫作限制性排序(constrainedordination)(1)冗余分析(redundancyanalysis,RDA)(2)典范对应分析(canonicalcorrespondenceanalysis,CCA)让我们来仔细看看PCA与PCoA分析：在非限制性排序中，16S和宏基因组数据分析通常用到的是PCA分析和PCoA分析。两者的区别在于：PCA分析是基于原始的物种组成矩阵所做的排序分析，而PCoA分析则是基于由物种组成计算得到的距离矩阵得出的。在PCoA分析中，计算距离矩阵的方法有很多种，例如图1所示的Euclidean,Bray-Curtis,andJaccard，以及图2显示的(un)weightedUnifrac(利用各样品序列间的进化信息来计算样品间距离，其中weighted考虑物种的丰度，unweighted没有对物种丰度进行加权处理)。图1Euclidean,Bray-Curtis,andJaccardPcoA图图2(a)weightedUnifracPCoA分析PCoA1维度将根内样品与根周边和土壤样品很好的分开；而PCoA2能够将不同土壤的样品很好的区分开来；(b)用pairwise-Bray-Curtis相似度对样品进行聚类。再来看看CCA与RDA分析：限制性排序主要有CCA分析和RDA分析。RDA基于线性模型，CCA则是基于单峰模型。一般我们会选择CCA来做直接梯度分析。但是，如果CCA排序的效果不太好，就可以考虑换做用RDA分析。RDA或CCA选择原则：先用species-sample资料做DCA分析，看分析结果中Lengthsofgradient的第一轴的大小，如果大于4.0，就应选CCA；如果在3.0-4.0之间，选RDA和CCA均可；如果小于3.0,RDA的结果要好于CCA。图3是CCA分析图，图中箭头代表不同的环境因子，射线越长表示该环境因子影响越大。环境因子之间的夹角为锐角时表示两个环境因子之间呈正相关关系，钝角时呈负相关关系。图3CCA分析图箭头代表不同的环境因子，射线越长表示该环境因子影响越大。环境因子之间的夹角为锐角时表示两个环境因子之间呈正相关关系，钝角时呈负相关关系。