使用CANOCO进行CCA或RDA

1、如何用如何用CANOCO 4.5 进行进行CCA/RDA分析分析 郑有用 由于制作者水平有限，本篇不涉及蒙特卡洛检验。进一步学习应参 考Canoco for Windows 4.5 中文简明教程 与Multivariate Analysis of Ecological Data Using Canoco。 本篇内容参考上述资料。 数据文件预处理 WCanoImp是将excel数据转换为CANOCO可以处理的数据格式的程序。 使用方法： 1、打开excel文件，将数据组织为“以行代表样方，以列代表变量”的格 式。 2、选中这一部分，并复制（此时这些数据在粘贴板上，WCanoImp将 直接从粘贴板

2、上读取这些数据）。 3、打开WCanoImp，点击保存（如果数据格式以如前所述在excel文 件中组织好，则不需要点击任何选框，直接保存即可）。保存文件后 缀为“dat”。 4、将环境变量和物种变量皆以上述步骤处理，得到两个dat文件，这 是我们进行排序的原材料。 RDA排序分析 CANOCO for Windows 4.5是我们进行排序操作的程序。 使用方法： 1、打开该程序，点击New project开始一次新的排序。 2、程序将出现一个名为”Available Data的对话框，选框意义如下图。选择 完毕后点选下一步（下面以进行RDA分析为例）。 同时有物种数据与环境数据时 选择（进行C

3、CA、RDA必须有 这两种数据）。 直接梯度分析，包括RDA、 CCA。进行这些分析时选择。 这里藏着一行，是间接梯度 分析，包括PCA、CA等。 3、进入“Data Files”对话框，用于指明将要进行处理的数据文件。指 定后点击下一步。 分别选择前面处理后的环境 与物种数据文件（dat文件）。 选择结果输出文件。 程序没有自动建立结果输出 文件的功能，因此需要我们 手动建立。建立方法是：建 立一个空的文本文件，然后 把后缀改为“sol”即可。建 立后再在本对话框中指定它。 4、进入“Type of Analysis”对话框，用于指明需要的分析类型。此处 以RDA为例。指定后点击下一步。 前

4、面若在”Available Data对话框 中选择了间接梯度排序，在此处 就可以选择这几种方法 因为我们前面在”Available Data对话框 中选择了直接梯度排序，在此处只能选 择几种直接梯度排序的方法，包括CCA、 RDA。 5、进入“Scaling： Linear Methods”对话框，此处参数的选择见 Multivariate Analysis of Ecological Data Using Canoco的P55，选择后 对图解释方式的影响见P152-153。指定后点击下一步。 关注样方 情况时选 择这个。 关注物种 时选择。 物种箭头长度代表物种 在排序空间内的变异性。 选择

5、这个时，物种 箭头长度代表了该 物种在多大程度上 被排序图解释。 样方情况 与物种情 况都想关 注时选择。 由于对原理的不了解，Species scores的选择对箭头长度的影响 我还不能很好地理解。 此外，这两个选项的选择还将对解释有下面影响： 前者：物种变量箭头在另一个物种变量箭头线上的垂直投影代表 了两者的相关性。 后者：两个变量箭头之间cosine值代表了两者之间的相关性。 但是，两种Species scores的在物种相关性方面解释的结果是差 不多的。 CCA与RDA类似，但由于模型的不同，有个别参数设置不同。在此 只介绍不同的部分 前面选择了CCA分析后，进入“Scaling： U

6、nimodal Methods”对 话框。 与RDA 部分相同 biplot scaling 适用于物 种比较集中，梯度较短 类型。 Hills scaling选项适 用于物种梯度较长尺度 类型。 一般选择前者就行，它 能为排序图提供一个更 量化的解释。 CCA中该对话框的选择 6、进入“Transformation of Species Data”对话框，此处参数的选择 见Multivariate Analysis of Ecological Data Using Canoco的P13。指定 后点击下一步。 一般进行Log转换，默认 值A、B为1的情况适应大 多数情况。但是当物种数 据的值太

7、小（比如0.1以下） 时，B=1会导致结果受参 数B影响过大。此时应该 改变A的值，比如说改为 10。 一般程序默认值在物种数 据大约在0100这个规模 下结果良好。 7、进入“Centering and Standardization”对话框，只有线性模型的分 析（PCA或RDA）才涉及此参数设置。此处参数的选择见Multivariate Analysis of Ecological Data Using Canoco的 P51。指定后点击下一步。 物种数据经过中心化和标准化后，各列变量的平均值均为0 ，方差均为1 对于样方，中心化和标准化并不总是必要的；对于物种数据，中心化是必要 的，标准

8、化不一定是必要的。 但是当物种数据量纲不同时，该数据一定要进行标准化! 而当物种数据量纲相同时，不标准化物种数据可能比较好，因为当不同物种 的平均丰度差别较大时，标准化会使得罕见物种的权重变大。 物种数据量纲相同时选 择这个。 物种数据量纲不相同时 选择这个。 8、进入“Data Editing Choices”对话框。这个对话框是用于对数据进 行选择操作的。指定后点击下一步。 用于在分析中 排除不想要分 析的环境变量 在分析中删除 不想要分析的 物种数据 如果数据都是想要分析的，则直接选择下一步即可。 9、进入“Forward selection of ”对话框。这个对话框是用于对环 境数据

9、进行预先选择操作的。指定后点击下一步。 并不是所有的解释变量对于响应变量的解释都有显著贡献。环境变量预 选的目的就是筛选出有显著解释性的环境变量。 手动预选，可以看 到各个环境变量的 解释量，从中选择 自动预选的参数 自动预选 与蒙特卡洛检 验相关，这部 分不理解。 10、此对话框为最终确认对话框，点选完成后参数将保存在一个con文 件中，选择好保存的文件即可。 11、出现Project对话框，点选Analyze即开始分析。 12、前面如果选择了手动预选，则此时进入Forward Selection Step对 话框。 总解释量 变量Salinity的解释量 当前已选变量所能解 释的解释量 用

10、Include variable按钮将选中的变量加入分析中，用Stop selection按钮结束 选择。 完成后得到分析结果。 13、分析结果显示在log对话框中，该log文件可以保存，在需要时还可 以再次打开。 这一部分数据有些文献会列出来，需要注意。 至此，我们得到了一个log文件和一个con文件，后者用于作图。 作图 CanoDraw for Windows程序用于作图。 使用方法： 1、打开分析后得到的con文件，并按提示保存作图结果在cdw文件中。 2、点选CreatBiplot and Joint PlotsSpecies and env. variables （根据希望做的图点

11、选即可，此处用物种与环境变量的叠图为例） 3、美化排序图。右键点击希望改变的部分，选择Select Suchlike可以 同时选中类似的部分，然后按下F5对属性进行改变。 4、在FileExport中选择希望导出的格式，即可得到成品图。 RDA/CCA方法的选择 两种方法的原理相同，但是使用的模型不同，需要先确定待分析数据适合 哪一种分析方法。方法是先进行DCA分析。 选择这个 在这里 使用方法： 1、选择DCA分析。指定数据文件等步骤与前文相同。 一般选择这个 即可 2、“Detrending Method”对话框。 3、转换。此部分参数设置与前文排序介绍中相同。 4、一路直接点击“下一步”

12、。 5、在“Project”对话框中直接点击分析 6、在“log”对话框中找到下图所示的一行 如果四个轴中梯度最大值超过4，选择单峰模型（CCA）比 较合适；如果小于3，线性模型（RDA）比较合适；介于 34 之间，两种模型均适合。 物种数据的量纲不同时不适合做单峰模型排序，即便梯度长度大 于 4。 结果解读 1、log文件的解读此处直接粘贴了中文简明教程，log文件的其他部分请参考该 资料。 2、如何解释排序图。 RDA排序图的解释 此处只给出一般较关注的物种与环境变量的关系图的解读，其他 部分参考中文教程。 物种箭头之间的关系 一般而言，箭头夹角的余弦值 代表两个物种的相关性，如近 直角则代表几乎不相关。 具体参考PPT的P8， “Scaling： Linear Methods” 对话框的选择。 物种和环境变量之间的关系 空心箭头代表环境变量，实心 箭头为物种。相应的余弦值代 表了之间的相关性。 CCA排序图的解释 物种和环境变量之间的关系 空心箭头代表环境变量，三角 形为物种。物种垂直投影于环 境变量延长线上，物种2距离 变量