利用R

利用R,r,mapping,r-factor,R,Mapping,R Factor,我有一个小数据框(见下文),包含两个物种“Aganope”和“Brienope”(因子格式),以及它的地理坐标(Aganope有19行记录,Brienope有12行记录)。我想绘制这个物种的地图,对于第一个物种,我的代码是: library(dismo) str(data) data data<-data[1:19,] str(data) data acgeo <- subset(data, !is.na(EST) & !is.na(SOUTH)) dim(acgeo) ac

我有一个小数据框(见下文),包含两个物种“Aganope”和“Brienope”(因子格式),以及它的地理坐标(Aganope有19行记录,Brienope有12行记录)。我想绘制这个物种的地图,对于第一个物种,我的代码是:

library(dismo)
str(data)
data
data<-data[1:19,]
str(data)
data
acgeo <- subset(data, !is.na(EST) & !is.na(SOUTH))
dim(acgeo)
acgeo[1:19, c(1:3,1:3)]
# Placer tous les lieux de récolte de Aganope sur le fond de carte bioclimate 5
points(acgeo$EST, acgeo$SOUTH, col="red", pch=21, bg="black", cex=1)
对!

请注意,在未调用绘图之前,无法调用点。代码中使用的符号和表中的符号还有一些其他问题

data <- read.csv() #Your raw data
subset<-data[data$Species_name == "Aganope",]
plot(subset$East, subset$South, col="red", pch=21, bg="black", cex=1)

数据此处未包含
数据的样本值,因此您的问题不存在。请阅读该链接,了解有关共享示例数据的提示,以便我们可以运行与您相同的代码。为什么在这里包含
library(dismo)
?您似乎没有使用该软件包中的任何函数。我猜你可能只是想要
acgeo[acgeo[,1]==“Aganope”,]
假设物种在第一列。你是否尝试过
newdata作为预付款谢谢你的回复。我在上一条消息中发送了我的数据帧。
data <- read.csv() #Your raw data
subset<-data[data$Species_name == "Aganope",]
plot(subset$East, subset$South, col="red", pch=21, bg="black", cex=1)