我想使用dismo
(R包)下载物种分布分析的gbif物种出现数据。
我已阅读指南,并在大多数情况下成功进行了分析。
我在这里使用了代码。
gbif("genus","species", geo=FALSE)
但是,我想将某些亚种分开进行分析。 我是否可以通过任何方式从某些亚种中提取信息?
答案 0 :(得分:0)
根据gbif
的文档和dismo
的插图,我们可以将*
附加到物种名称,以获取该物种名称下的所有变体。以下示例来自小插曲,下载 Solanum acaule 下的所有记录。
acaule <- gbif("solanum", "acaule*", geo = FALSE)
结果数据框acaule
包含115个变量,大约有7000条记录。在这些变量中,scientificName
列包含该记录的全名。我们可以使用该列来识别与某个亚种相关的记录。
例如,我们可以先使用table
函数查看每个物种名称的数量。
table(acaule$scientificName)
Solanum acaule Bitter
5608
Solanum acaule subsp. acaule
1444
Solanum acaule subsp. punae (Juz.) Hawkes & Hjert.
14
Solanum acaule var. punae (Juz.) Hawkes
9
Solanum acaule var. subexinterruptum Bitter
2
Solanum schreiteri Bukasov
2
Solanum uyunense Cárdenas
2
现在我们可以根据scientificName
中的关键词对数据框进行子集化,以找到我们想要的记录。假设我们只想在subsp.
中使用scientificName
的记录,我们可以执行以下操作:
acaule_sub <- subset(acaule, grepl("subsp.", scientificName, fixed = TRUE))
现在acaule_sub
是一个包含1458条记录的数据框,这些记录都是科学名称与“subsp。”相关联的记录。
不确定这是否是您所需要的,但我希望它有所帮助!