数量生态学笔记||典范排序概述

起了个大早做这张导图,不知道会不会吓到大家。有一周没有更新我们的《数量生态学笔记》了,因为这一章还是蛮有难度的,书中就是说这一章我们来讲典范排序,然后就直接讲RDA了。看RDA一会排序一会回归,书中的每一个汉字都认识,放在一起的一个句子就不懂了。这基本上就是我的一个状态,所以就没更新。

本来以为只要多查些资料就好了,多看嘛,总会理解的。看张金屯老师的《数量生态学》有知道典范排序(约束排序)是在非约束排序分析的每一步与环境因子进行回归,再将回归系数结合到下一步的排序值得计算之中。其中与环境因子的结合方式如下:

\begin{equation} y_j=b_0+b_1z_{1j}+b_2z_{2j}+b_3z_{3j}+...+b_qz_{qj} \text {,(1)} \label{eq:Sample} \end{equation}

式中,y_j为第j个样方的排序值;b_0为截距;b_1为第i个环境因子的回归系数(i=1,2,3..q 为环境因子的数目)可用多元线性回归求得;z_{ij}为第i个环境因子的观测值。

可以看出,典范排序不再是被动加入解释变量,而是从排序开始直接加入解释变量进行计算。与非约束排序不同的是,典范排序只是提出和展示与解释变量有关的数据结构,并可以通过统计检验方法检验解释变量与响应变量之间关系的显著性

不管怎么说我们讲典范排序就离不开多元回归模型,因为

回归 + 非约束排序 = 典范排序

线性多元回归模型

所以我找来同样是赖老师翻译的《基于 CANOCO 的生态学数据的多元统计分析 》第五章5约束排序与置换检验(Constrained Ordination and
Permutation tests),在这里和大家分享一下其中的部分内容。

首先,我们必须回顾一下传统的线性回归模型,因为这对于我们理解“直接梯度分析”(约束排序)相当重要。图5-1展示的是最简单的线性回归模型,线性模型可以模拟响应变量Y依赖自变量X的程度。图5-1中不仅有拟合回归线,也展示了模拟值和实测值之间的差别。模拟值\hat{Y_i} (回归线上的值)与实测值Yi之间的差值叫做回归残差(regression residual),用e表示。所有的统计模型(statistical models,包括回归模型)有个重要的特征是它们都有两个主要的部分构成:系统组成部分(systematic component)表示响应变量中能被一个或更多的解释变量(模型)解释的部分,这部分用带参数的函数表示。另外一部分就是随机部分(stochastic
component),表示不能被目前解释变量(模型)所能解释的部分。随机部分通常用概率和分布特性来定义。

图 5-1

我们通常通过响应变量有多少能够被系统组成部分解释来评估拟合模型的好坏。也经常将能被解释和未被解释的部分进行比较。目的在于,尽力去构建一个最简约的回归模型来解释最多的变化量,让所有的自变量对于响应变量的解释都有显著贡献。我们可以通过逐步迭代(回归)(stepwise selection)的方式来选择解释变量(环境变量)的子集合,在排序术语中往往叫预选(forward selection)。预选变量的过程是无响应变量的零和假设开始,零和假设是响应变量中没有可以被解释变量预测,而仅仅由随机变量解释。当我们选择一个解释变量(环境变量)进入分析,可以导致回归模型能解释一部分响应变量。可以根据所加入的变量所能解释部分的大小来确定是否选择的该环境变量。另外,需要用随意置换(randomly swap)环境
因子的值来检验这种解释量是偶然的,还是真的为环境变量所解释? 如果被检验的变量所能解释的部分被证明是非随机的(统计显著), 我们就可以接受这个变量。可以重复这个过程,再进一步从剩下的变量中选择另外更好的变量,直到选择具有足够的变量为止。

约束排序模型 (Constrained ordination model)

非约束的排序(PCA 和CA)被定义为寻找潜在的梯度代表最优的解释变量(预测器)来拟合物种的回归模型。

约束排序跟非约束的排序有一个很大的区别,非约束排序是虚拟的(潜在的)梯度、而约束排序的梯度是明确给出的。这些梯度(排序轴)是参与排序的环境变量的线性组合。因此我们通过合成变量(排序轴)尽力解释物种的多度变化,这些合成变量是实测环境变量的线性组合。

因此,约束排序方法(RDA和CCA)类似于多元多重回归。但在多元多重回归中,如果有m个响应变量,p个环境因子,我们必须估计出m×p个的参数(回归系数)(每个方程需要p个参数,m个方程自然是需要m×p个的参数)。然而,在约束排序里面,不必这么麻烦,环境因子对于响应变量的影响被集中在几个合成的梯度(排序轴),也叫典范轴(canonical axes)。典范轴的数量是跟独立解释变量的数量一样多,但是我们经常使用前面几轴。如果环境因子和物种的数量很少的时候,是不必做约束排序的,做普通回归即可,mp比较多的时候,约束排序的优势才能展示出来。

如果有协变量(covariables),我们经常使用偏分析(partial analyses)。有协变量情况,表示我们要将这些协变量的所能解释的部分先剔除出去。协变量在方差分析中也有相同的用法,通常是把量化的协变量作为一种因子处理。而在传统的回归中,协变量的概念是不常用的,协变量与真变量没有什么不同,叫法不同而已。

别怕,这些概念都会在我们的典范排序的笔记中给到大家的。不管怎么排序,就记住一条就好:在排序空间内,离得越近的关系越近。

参考:
Ordination Methods - an overview
Canonical ordination
基于 CANOCO 的生态学数据的多元统计分析|第五章
Multiple Linear Regression - MLR
Multiple Linear Regression

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

推荐阅读更多精彩内容