微生物分析-排序方法

约束性/非约束性排序

主成分分析(PCA)、对应分析(CA)(加权平均法迭代运算中)、去趋势对应分析(DCA)和非度量多维尺度分析(NMDS),冗余分析(RDA)约束的PCA 和典范对应分析(CCA)

1. PCA和PCoA

主成分分析PCA(Principal component analysis)是一种研究数据相似性或差异性的可视化方法,采取降维的思想,PCA 可以找到距离矩阵中最主要的坐标,把复杂的数据用一系列的特征值和特征向量进行排序后,选择主要的前几位特征值 ,来表示样品之间的关系。通过PCA 可以观察个体或群体间的差异。PCA后面的百分数表示对应特征向量对数据的解释量,此值越大越好;

PCoA的计算原理与PCA相同,只是不象 PCA 只用欧氏距离方程计算点间距离,它可以用各种距离系数,从上面的分析过程看,PCoA与PCA有很大的相似性。可以证明,如果使用欧氏距离矩阵 PcoA与使用种类化中心的PCA正分析相同

对应分析(Correspondence analysis, CA)(Hill 1974)也叫做相互平均法,但它有一个重大缺点,就是CA/RA的第二排序轴在许多情况下是第一轴的二次变形,即所谓的“弓形效应”(Arch effect)或者“马蹄形效应”
The CA distortion is called the arch effect, which is not as serious as the horseshoe effect of PCA because the ends of the gradients are not incurved. Nevertheless, the distortion is prominent enough to seriously impair ecological interpretation.(注:PCA存在“马蹄形效应”,搞出了CA,CA存在“弓形效应”,搞出了DCA)

马蹄形效应

弓形分析

2. PCoA和NMDS

一般的排序方法都要求P(种)×N(样方)的原始数据矩阵,再从原始数据出发去进行排序。但是还有一类排序方法不是从原始数据出发,而是以样方间相异距离矩阵为起点,这一距离矩阵可以是基于某个相异系数而计算出来的,也可以是通过观察直接得到的样方间的某种差异。这样排序的目的是将N个样方排列在一定的空间,使得样方间的空间差异与原始距离矩阵保持一致,这类排序方法称作多维标定排序(Multi-dimensional scaling)。如果排序依赖于相异系数的数量值,就叫有度量多维标定法(Metric multi-dimensional scaling);如果排序仅仅决定于相异系数的大小顺序,则称为无度量多维标定法(Non-Metric Multi-Dimensional Scaling;NM-MDS)。前面讲到的主坐标分析就是一个有度量多维标定法(张金屯 1995,余世孝 1995)。

在排序空间给样方布局,让样方间距离的远近度量代表它们之间物种组成的相似度。这类排序方法称作多维标定排序(Multidimensional scaling) 如果排序依赖于相异系数的数量值,就叫有度量多维标定法(metric mult-dimensional scaling)PCOA;如果排序仅仅决定于相异系数的大小顺序,则称为无度量多维标定法(Non-Metric Multi-Dimensional Scaling;NM-MDS)

NMDS:试图在一个低维的空间中建立一个样品分布地图,是样品间的欧式距离的等级顺序与其相似性或非相似性等级顺序保持一致。等级相似性rank order similarity:将相似性矩阵相似性按大小排序,最高的相似性给与最低的等级1,最低的相似性给与最低的等级你(n-1)/2。

A measure of ‘stress’ (mismatch between the rank order of distances in the data, and the rank order of distances in the ordination) is calculated,在无度量多维标定法,我们不用影射规则,而是相异系数的大小顺序,通过数量算法球得样方间的距离。
排序后样方的空间差异性与原距离矩阵差异最小或者
要使得N个样方在m维排序空间中的距离与预先已知的距离保持一致

非约束排序也叫间接梯度分析,而约束排序也叫直接梯度排序分析。
非约束排序可以下面两种方式构建:

间接梯度分析的目标就是发现这样的坐标轴,让群落中的样方或是物种的最大变化量能够在坐标轴上体现出来。换句话说,让尽可能多的变化量能够在尽可能少的轴上展示出来,并且让样方或物种在排序图能够可视化展示出来。当然,我们会经常期望这些轴能够代表一些潜在的环境变量。而约束排序(constrained ordination)的目的就是发现物种在环境梯度上的变化情况。

所有排序方法都是基于一定的模型之上,这种模型反映植物种和环境之间的关系以及在某一环境梯度上的种间关系。最常用的关系模型有两种:一种是线形模型(linear model),另一种是单峰模型(unimodal model)。线形模型的含义表示某个植物种随着某一环境因子的变化而呈线性变化或叫线性响应(linear response)。单峰模型的含义是某个植物种的个体数随某个环境因子值的增加而增加。当环境因子增加到某一值时,植物种的个体数达到最大值,此时的环境因子值称为该种的最适值(optimum);随后当环境因子值继续增加时,种的个体数逐渐下降。为了简化单峰模型,我们经常假设单峰曲线以峰值的为中心,两边是对称的(见Hutchinson 1957,资源梯度与物种最适值)。在排序之前,我们首先决定用那种模型。大部分情况下,两种模型都仅仅是实际数据的近似拟合,所以我们要决定的那种模型对我们所要分析的拟合更好。其实,单峰模型也是理想化的模型,在现实中,物种对环境才响应很少以最适值为中心两边对称的,更复杂的响应曲线也经常出现(比如双峰模型)。

进行排序分析之前,首先要判断是选择线性模型(PCA 和RDA)还是单峰模型(CA和CCA)的排序方法。一般来说,如果物种分布变化大,选择单峰模型效果比较好,反之,线性模型也是不错。可以通过DCA分析来判断,如果DCA排序前4个轴中最大值超过4,选择单峰模型排序更合适。如果是小于3,则选择线性模型更好(Lepx & Smilauer 2003)。如果介于3-4之间,单峰模型和线性模型都可行。


环境梯度与模型选择
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,588评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,456评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,146评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,387评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,481评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,510评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,522评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,296评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,745评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,039评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,202评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,901评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,538评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,165评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,415评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,081评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,085评论 2 352