rda冗余分析步骤_数量生态学笔记||冗余分析（RDA）

上一节数量生态学笔记||冗余分析(RDA)概述中，我们回顾了RDA的计算过程，不管这个过程我们有没有理解透彻，我希望你能知道的是：RDA是响应变量矩阵与解释变量之间多元多重线性回归的拟合值矩阵的PCA分析。本节我们就是具体来看一个RDA的分析案例，来看看里面的参数以及结果的解读。

# CHAPTER 6 - CANONICAL ORDINATION

# ********************************

# 载入所需程序包

library(ade4)

library(vegan)

library(packfor)# 可在http://r-forge.r-project.org/R/?group_id=195下载，但是好像在R 3.5.1上加载不了，所以这篇我用R3.4来做的。packfor已经不用，函数都搬到adespatial

# 如果是MacOS X系统，packfor程序包内forward.sel函数的运行需要加载

# gfortran程序包。用户必须从"cran.r-project.org"网站内选择"MacOS X"，

# 然后选择"tools"安装gfortran程序包。

rm(list = ls())

setwd("D:\\Users\\Administrator\\Desktop\\RStudio\\数量生态学\\DATA")

library(MASS)