- 博客(18)
- 收藏
- 关注
RSS订阅原创 R语言|β 多样性分解以及三元相图的绘制
物种周转通常是生境过滤、竞争或历史事件影响的结果，而丰富度的差异则可能是由于物种自疏产生的群落嵌套性(Legendre, 2014)。将总beta多样性（β_total）分解为物种周转（β_turnover）和物种嵌套（β_nestedness），并以三元相图的形式展示，是在生态学beta多样性分析中常用的一种数据分析方法。在低版本的R中，ggtern包常常与ggplot包有冲突（我的就是），要注意。beta.sor： total beta（总beta）beta.sim：turnover（物种替换）
2026-03-04 21:07:21
335
原创 生态学|R语言-基于距离的RDA分析（dbRDA）--两种不同方法
普通RDA是获取群落数据与环境因子之间的线性关系区别就在于，RDA是基于群落多度做的，dbRDA是先基于Bray-Curtis距离（或其他距离）进行主坐标分析（Principal Coordinates Analysis，PCoA），随后将特征值大于0的坐标轴与环境因子进行RDA分析，即dbRDA。Eigenvalues, and their contribution to the variance 部分的解释率是指RDA轴对总体（total）的解释率，因此在一些同学的分析中会出现这部分数值非常低的情况。
2026-01-25 23:58:39
600
原创 如何在R语言中实现克里金插值（Kriging）——代码篇
上一篇已经讲述了kriging插值的原理，最重要的就是构建半变异方差函数，并选择拟合效果最佳的函数，进而获得对应的这四个指标也是在R语言中构建半变异方差函数最关键的指标。这里以以500m*500m的样地海拔为例。样地建设之初在许多采样点（+）测定了海拔（小样方的交界处），现在我想要获取每个小样方中心点的海拔数据。当然，可以对四个角的海拔取平均=中心点海拔，这里我们使用kriging插值方法ArcGIS中kriging插值的步骤是：插值→将结果转为raster栅格→按小样方中心点坐标提取数值，
2026-01-20 20:22:48
1314
原创 如何在R语言中实现克里金插值（Kriging）——原理篇
经过几周的摸索，终于把kriging插值原理弄懂了，并且通过对比ArcGIS与R语言的插值结果，我认为在R中可以更方便快捷的进行ordinary kriging ,因此想用两篇的内容详细讲解一下克里金插值。首先尝试的是在ArcGIS中进行普通克里金插值，在网上也可以搜到很多教程，作为一个地理小白，我起初并不太清楚怎么样在ArcGIS中调整插值的各种参数。函数的准则就是：距离近的物体比距离远的更相似，因此半变异方差函数是一个与距离相关的函数（如下图，即下图中的Y轴，X轴为距离。在实际应用中自行选择模型。
2026-01-19 23:46:34
1113
原创 生态学|R语言-NMDS分析检验群落间物种组成差异
PCA基于欧氏距离，适用于连续型数值数据，且要求线性关系。NMDS基于秩次（rank-order）的非参数方法，适用于非连续数据（如分类数据、序数数据）或非线性关系的数据。NMDS（非度量多维尺度分析）和PCA（主成分分析）都是常用的降维技术，但在不同数据特征和研究目标下，选择的方法有所不同。若研究目标是比较样本间的相似性（如群落相似性排序），NMDS的二维/三维图通常更直观。行名为样地-小样方（plot1-subplot1），列名为物种，数值为多度。不同样地的NMDS图。
2026-01-16 22:00:11
401
3
原创 R语言|将多个tif文件批量提取到点并计算生长季长度
R语言|将多个tif文件批量提取到点并计算生长季长度想要基于日均温计算生长季长度时，步骤是以提取2000-2018年的日均温数据为例，得到的数据常常是年份命名的文件夹。每个文件夹内有365个tif数据，在arcgis中提取的话十分繁琐。首先加载所需包读取采样点数据。
2026-01-15 23:55:39
266
原创 R语言|处理nc数据并将气候数据提取到采样点
我们常常需要将nc文件按照所需时间提取为.tif文件以便后续的处理，比如将气候数据提取到多个采样点。nc数据的提取可以在ArcGIS、MATLAB中完成，但对于多个nc文件来说，ArcGIS的提取过于繁琐
2026-01-14 23:59:58
290
原创 生态学|R语言-Beta多样性分解：物种周转和嵌套组分
上一篇讲到如何构建系统发育树，这篇讲群落β多样性以及系统发育多样性的计算与分解。生态学|R语言-β多样性分解：物种周转和嵌套组分。
2026-01-14 23:55:13
402
原创 R语言|单因素方差分析+TukeyHSD多重比较
接上一篇单因素方差分析，这篇介绍TukeyHSD多重比较方法（Tukey Honest Significant Differences），并给出画图代码。控制实验中，对比多个物种在不同处理下的生长量，想要将所有物种的生长量组间差异都绘制在一张图上。R语言|单因素方差分析+TukeyHSD多重比较。
2026-01-13 21:08:12
406
原创 生态学|R语言-单因素方差分析(ANOVA)及LSD多重比较
在R语言中，使用aov函数进行单因素方差分析， 使用agricolae包中的LSD.test函数进行Least significant difference检验。比如进行种源实验时，将样地A、B、C采集的物种带回后在气候箱中培养1周、2周、3周，记录其生长量，对比相同培养时间下，物种的生长率在不同样地间是否存在差异。当收集的数据有多个嵌套分组且有多个变量时，想要在某个组别内对比变量在另一分组间的差异 ，并进行LSD多重比较。下一篇介绍如何进行Tukey多重配对比较并绘制条形图。
2026-01-12 23:58:12
459
原创 生态学|R语言-可视化相关性图
生态学分析中，相关性分析常常是探索数据的第一选择，有时也能够为论文提供支撑性的结果。本文给出了可视化相关性分析的R代码
2026-01-11 20:31:56
373
生态学森林植物群落调查表
2026-01-11
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人