R语言多元数据统计分析在生态环境中的应用

生态环境领域研究中常常面对众多的不同类型的数据或变量，当要同时分析多个因变量（y）时需要用到多元统计分析（multivariate statistical analysis）。多元统计分析内容丰富，应用广泛，是非常重要和实用的多元数据分析方法和统计工具，其中分类(classification)/分组（grouping）和梯度（gradient）/排序（ordination）分析是多元统计分析的核心内容。分类/分组分析主要包括聚类（如层次聚类和k-means聚类等）和辨别分析（如线性辨别分析和mental检验等）；梯度/排序分析分为非约束排序（如PCA和CA等）和约束排序（如RDA和CCA等）两大类（详见课程介绍）。但初学者在应用多元统计分析方法时往往非常困惑，这种困惑主要来自几个方面：首先，多元统计方法众多，分类还是排序？约束排序还是非约束排序？哪种方法或技术更适合我的研究目的或数据？其次，在多元分析中很多术语都有别称，如非约束排序也叫间接梯度分析；再次，多元数据类型包括连续型数据、计数数据、分类数据及混合类型数据，不同数据类型如何选择合适的方法？等等等等。将梳理多元统计中分类/分组和排序/梯度分析各方法不同应用情境，包括生态环境多元数据分析概述、R基础及数据探索、多元数据梯度分析、多元数据分组分析及统计结果作图。课程将通过大量具体案例，利用R语言相关包示范各种方法实现途径，使大家面对多元数据分析时能够驾轻就熟，从容面对。不仅适合R语言和生态环境数据多元统计分析的初学者，同样适合有高阶应用需求的研究生和科研人员。

查看原文

第一篇（基础篇课前预习，Ai尚研修提供预习资料） 生态环境多元数据分析概述 (Working with multivariate data) 

1、多元统计方法概念定义

2、纷繁的多元统计方法应用情景、异同点

3、多元统计方法数据或变量类型和结构

第二篇 R基础及数据探索 (Working with R and Data screening)

专题一：R和Rstudio简介及入门和作图基础

1） R及Rstudio介绍：背景、软件及程序包安装、基本设置等

2） R语言基本操作，包括向量、矩阵、数据框及数据列表等生成和数据提取等

3） R语言数据文件读取、整理（清洗）、结果存储等（含tidverse）

4） R语言基础绘图（含ggplot）：基本绘图、排版、发表质量绘图输出存储

专题二：群落数据准备及探索分析

1) 生物群落数据准备：物种组成、环境变量、物种功能属性、系统发育树等

2) 生物群落数据检查：缺失值和离群值（outliers）等-避免模型错进错出（GIGO）

3) 物种多样性计算：物种多样性（TD）、功能多样性（FD）和系统发育多样性（PD）

4) 物种相似/相异矩阵关联测度介绍

第三篇 多元数据分组分析 (Working with groups)

专题三：聚类(Cluster)分析上：非层次聚类（NHC）

1) 聚类及非层次聚类方法概述

2) 非层次聚类：K均值(K-means)聚类方法（kmeans;pam;clara）

3) 实例鸟类生境数据K均值聚类比较分析：聚类数确定、聚类稳定性、聚类结果评估、复合聚类值构建

专题四：聚类(Cluster)分析下：层次聚类（NHC）

1）层次聚类方法简介：多元聚合层次聚类（PAHC） VS 多元分化层次聚类（PDHC）

2）层次聚类方法（hcluster和agnes）分类结果比较分析：碎石图、轮廓宽度、同表型相关图等

3） 案例1鸟类生境数据的层次聚类分析；案例2鱼类生境数据的层次聚类分析

专题五：辨别(Discrimination)分析上：分组差异性检验（Group Difference Test）

1） 多元群落数据分组差异分析及检验简介

2）（非参数）多元方差分析（NP-MANOVA/ADONIS/PERMANOVA）、多元置换过程（MRPP）、多元相似分析（ANOSIM）、Mantel检验（MANTEL）在多元数据差异性检验的应用

3）多元差异性检验多组数据比较实现方法：MRPP、Mantel

4）Mental方法在群落生态学‘标准’应用：空间取样距离、环境因子及物种组成关系及其偏Mental分析

5） 案例1乌龟适生生境差异性检验；案例2微生物组成数据差异分析；案例3鱼类群落、空间距离及环境因子相互关析分析

专题六：辨别(Discrimination)分析中：线性辨别分析（LDA）

1） 辨别分析（DA）的多面性

2） 线性辨别（LDA）分析基本原理及流程：数据检查、评估假设、样本数量、变量选择、模型确定、结果解读及模型验证

3）其他辨别分析方法介绍（QDA，KNN等）

4） 案例乌龟适生生境辨别及预测

专题七：辨别(Discrimination)分析下：分类回归树（CART）及随机森林模型（RFM）

1） 生物群落数据的分类回归树简介

2） 分类回归树分析（CART）实现：分化准则、节点杂度、gini指数、先验概率效应、误分类代价、分类树裁剪、Monte Carlo检验、变量重要性评估、模型预测等

3） 随机森林模型（RFM）实现：算法流程、模型评估、变量重要性评估及模型分类与回归等

4） 案例1基于分类回归树的乌龟群落生境划分及预测

5）案例2基于随机森林模型的根际微生物群落与植物生长关系及变量重要性评估

第四篇 多元数据梯度分析/(Working with gradients)

专题八：间接梯度分析-非约束排序（Unconstrained ordination）上：PCA

1）生物群落数据的非约束排序简介

2） 主成分分析（PCA）的基本原理：假设条件、数据要求等

3）案例：鱼类生境数据的PCA排序分析实现-数据准备、检查（离群值、多元正态性、线性关系、样本独立性等）、结果验证、排序轴选择（特征根准则、累计解释率、随机断棍准则等）、结果解读、双序图等

专题九：间接梯度分析-非约束排序（Unconstrained ordination）下：PCoA、CA、DCA 及NMDS

1） 其他非约束方法介绍及应用情景： CA、DCA 、PCOA及NMDS

2） 案例1鸟类群落组成数据的对应分析（CA）及去趋势对应分析（DCA）：数据准备、假设条件、总惯量、特征根、排序轴选择、结果解读、偶见种效应/弓形效应等

3） 案例2 基于距离/相似度指数或矩阵的主坐标排序（PCoA）：距离/相似度指数选择、模型假设、负特征根问题、结果解读、排序图等

4） 案例3 NMDS排序方法应用：假设条件、基本分析流程、排序效果评估（应力值）、排序图等

5）案例4 药物对肠道微生物群落影响：PCoA+PERMANOVA

6）案例5 基于随机森林模型的蚂蚁多维属性特征预测：RF+PCA+PCoA+PERMANOVA

专题十：直接梯度分析-约束排序（Constrained ordination）上：RDA

1） 生物群落数据的约束排序简介：非对称约束排序 VS 对称约束排序

2）生物群落数据非对称约束排序的基本流程：响应变量/物种选择（矩阵Y）、响应变量数据预处理（转换或标准化）、分析方法选择（RDA/db-RDA/CCA）、解释变量/约束变量选择（矩阵X）及分析和结果解读、评估及展示

3） 案例景观、斑块及立地条件对森林景观中蛾类群落物种组成的影响

专题十一：直接梯度分析-约束排序（Constrained ordination）中：dbRDA、CCA及对称约束排序方法

1）案例1：基于距离矩阵的冗余分析(dbRDA)：物种组成数据及0，1数据分析

2）案例2：群落物种多度数据的典范对应分析（CCA）：单峰型环境梯度分析方法

3）案例3：对称约束排序方法介绍及物种组成、物种属性及环境变量相关关系第四角分析

专题十二：直接梯度分析-约束排序（Constrained ordination）下：变差分解（Variance Partitioning）

1） 群落数据多元统计分析变差分解简介        

2） 偏回归分析与变差分解

3） 案例景观、斑块及立地条件及空间因素对森林景观蛾类群落物种组成变异的变差分解

第五篇 统计结果作图ggplot (Plotting the results)

1. 群落数据及统计分析结果作图数据准备：结果提取、整理
2. PCA、CA、PCoA及NMDS等非约束排序图：排序图和双序图（biplot）
3. PCoA+PERMANOVA结果图：排序图+分组+PERMANOVA差异显著性+多重比较
4. RDA、db-RDA及CCA等约束排序图：三序图（triplot）和韦恩图（venn）
5. 

推荐：