|
if(!require(devtools)) install.packages("devtools")
devtools::install_github("helixcn/plantlist")
library(plantlist)
# 读入大样地数据
xueming <- read.csv("E:/guiides/species.csv",header = T)
# 把数据中是拉丁学名的一行或者一列,提取出来
species_names<-xueming$scientific.name
# 统计物种数量
species_number <- length(unique(species_names))
# 依据上面的拉丁学名,函数TPL()查询目、科、属、以及科的编号
taxa_genus_family <- TPL(species_names)
# 统计所有物种所在的属数
genus_number <- length(unique(taxa_genus_family$POSSIBLE_GENUS))
# 统计所有物种所在的科数
family_number <- length(unique(taxa_genus_family$FAMILY))
# 统计所有物种所在的目数
order_number <- length(unique(taxa_genus_family$ORDER))
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-4-24 23:39
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社