R计算功能多样性— functional diversity (FD)
作者&投稿:荤师 (若有异议请与网页底部的电邮联系)
为研究中包括的每个物种创建一个定性和/或定量性状的矩阵;使用R中的FD包来计算你选择的性状的种间Gower距离。Gower距离矩阵可用于绘制树状图,描述物种间的相似性/不相似性,然后你可以使用R中的FD或BAT或其他软件包,以类似于计算系统发育多样性(PD)的方式,总结连接每个物种库中的物种的分支长度。 参考: https://www.researchgate.net/post/Can_anyone_help_with_measuring_functional_diversity
FD computes different multidimensional FD indices. Tocompute FD indices, a species-by-trait(s) matrix is required (or at least aspecies-by-species distance matrix). gowdis computes the Gower dissimilarityfrom different trait types (continuous, ordinal, nominal, or binary), andtolerates NAs. It can treat ordinal variables as described by Podani (1999),and can handle asymetric binary variables and variable weights. gowdis iscalled by dbFD , the main function of FD.】
dbFD uses principal coordinates analysis (PCoA) toreturn PCoA axes, which are then used as ‘traits’to compute FD. dbFD computesseveral multidimensional FD indices, including the three indices of Villéger etal. (2008): functional richness (FRic), functional evenness (FEve), andfunctional divergence ( FDiv ). It also computesfunctional dispersion (FDis) (Laliberté and Legendre 2010), Rao’s quadraticentropy (Q) (Botta-Dukát 2005), a posteriori functional group richness (FGR),and the community-level weighted means of trait values (CWM), an index offunctional composition. Some of these indices can be weighted by speciesabundances. dbFD includes several options for flexibility.
基本命令: dbFD(x, a, w, w.abun = TRUE, stand.x = TRUE,..)
x: 功能特征的矩阵或数据框架
a: 包括出现在x中的物种丰富度矩阵。The number of species (columns) in a must match the number of species (rows) in x. In addition, thespecies labels in a and x must be identical and in the same order. (x和a两个数据集的物种数目一致和顺序相同)
建立两个数据dataframe, 保持两个数据集中的物种名称顺序保持一致:
情况1:手动录入的数据,或者数据量少的数据可手动调节;
情况2:物种较多,两个数据集包含较多的物种名称不一致情况;
解决办法:使用excel,若X和a中的物种不完全相符,以物种较少的数据物种数目为准,删减x中的数据以对应a。利用数据透视表功能建立数据矩阵; 注意: a中物种名称若出现空格,将不能成功读取,会出现格式变化,建议标准物种名称为“Abies_chensiensis”(加短横线)。
library(“FD”) #加载FD包
traits <- read.csv("D:/ xxx.csv", row.names=1) #x,加载数据
species <- read.csv("D:/xxx.csv", header=T,row.names=1) #a,加载数据
rownames(traits) #查看导入数据的行列名称,是否有不一致情况
colnames(traits) #查看导入数据的行列名称,是否有不一致情况
rownames(species) #查看导入数据的行列名称,是否有不一致情况
colnames(species) #查看导入数据的行列名称,是否有不一致情况
results = dbFD(traits, species, corr="lingoes") #dbFD计算FD值
解决办法:a =dbFD(traits, species, corr="lingoes")
write.csv(results,file = "D:/ xxx-FD-results.csv") #数据导出,结果导出为csv格式数据
## FD一般是FR,FE,FD三个层面,FR指功能多样性丰富度,FE值功能多样性均匀度,FD值功能多样性离散度。
温承克咳: R是一套完整的数据处理、计算和制图软件系统.其功能包括:数据存储和处理系统;数组运算工具(其向量、矩阵运算方面功能尤其强大);完整连贯的统计分析工具;优秀的统计制图功能;简便而强大的编程语言:可操纵数据的输入和输出...
凤县14788846342: matlab与R比较,两者的长处分别在哪? - ?
温承克咳: ①R 功能R是一套完整的数据处理、计算和制图软件系统.其功能包括:数据存储和处理系统;数组运算工具(其向量、矩阵运算方面功能尤e79fa5e98193e78988e69d8331333264663032其强大);完整连贯的统计分析工具;优秀的统计制图...
凤县14788846342: R软件中如何进行群落聚类分析? - ?
温承克咳: 群落按照物种相似形组成进行聚类分析,可以用树状图较好的表现物种的组成关系.受到很多植被学家的重视.这里以R软件实现聚类分析为例. 如果按照物种组成的相似性做聚类分析,那么可以用Jaccard指数(经过转换的).Jaccard指数只...
凤县14788846342: 怎么使用这个.铣床上的电子尺上面有个铣R的功能键?铣床上的电子尺 ?
温承克咳: 1、 进入R功能2、 选择SIMPLE(简易R)3、 选择加工形式(WHICH)① XY (ARC-XY)② XZ (ARC-XZ)③ YZ (ARC-YZ)4、 输入圆弧的半径(RADIUS)5、 输入刀具直径(TL DLA)6、 输入最大切削量(MAX CUT)
凤县14788846342: win10唤醒快捷键 - win10电脑唤醒快捷键 ?
温承克咳: 1. win10电脑唤醒快捷键方法一:适用于休眠后打开屏幕,使用一段时间后.长时间... 可以进行数值计算,又可以进行逻辑计算,还具有存储记忆功能.是能够按照程序运...
凤县14788846342: R型指令需要由op字段和func字段共同决定指令的功能 - 上学吧普法考试?
温承克咳: int fun(int n,int m) {int i,sum;sum=0;for(i=n;i<=m;i++){if(i%2==0)sum=sum+i;}return(sum); } #include<stdio.h> void main() {int n,m,t,h;printf("请输入两个正整数:\n");scanf("%d%d",&n,&m);if(n>m){t=n;n=m;m=t;}h=fun(n,m);printf("%d\n",h); }经过tc和vc++运行后均可出来正确结果
凤县14788846342: R 和 Python 用于统计学分析,哪个更好 - ?
温承克咳: 2012年的时候我们说R是学术界的主流,但是现在Python正在慢慢取代R在学术界的地位.不知道是不是因为大数据时代的到来.Python与R相比速度要快.Python可以直接处理上G的数据;R不行,R分析数据时需要先通过数据库把大数据转化...
凤县14788846342: SAS EVIEWS MATLAB R 这四个软件哪个用处最广??
温承克咳: SAS是统计软件EVIEWS是计量回归软件MATLAB是数值、图形计算,仿真、工程,综合软件R也是统计计量的软件我觉得matlab最广吧
凤县14788846342: 高一生物必修一中的蛋白质怎样好理解? - ?
温承克咳: 你可以这样看,每个氨基酸,都有一个氨基和羧基,都有一个碳和氢.而R基是区别每个蛋白质的方法.R基就是侧链基团…… 就像一群人,都有眼睛,鼻子,嘴巴,手脚,好比蛋白质都有氨基羧基. 但每个人长相不同,好比不同的侧链基团...