Ecology letter coordinate value check

read data

data_smith<-fread("DataforPayne_Smith_EL.csv")
data_sunday<-fread("sd01.csv")

篩掉“_sp“,”_sp1“,”_sp2“等以上未明確定義之物種

data_sunday_filter<-data_sunday %>% 
                    filter(!grepl("_sp",data_sunday$species))

data_smith與data_sunday_filter join

data.merge<-merge(data_smith[,-(14:16)],data_sunday_filter[,c("species","lat","longitude")],by.x= "Species",by.y="species")

rename columns

data.merge<-as.data.frame(data.merge)
str(data.merge)
colnames(data.merge)[which(names(data.merge) == "lat")]<-"latitude"

colnames(data.merge)[which(names(data.merge) == "Absolute latitude")]<-"Absolute_latitude"

add new column {“absolute latitude” - “lat”} for check the Coordinate accuracy

data.merge.accuracy<-data.merge %>% 
                   mutate(accuracy = Absolute_latitude-abs(latitude)) %>%
                   select(Species,Absolute_latitude,Reference,latitude,longitude,accuracy)%>%
                   filter(accuracy >10 ) 

data.merge.accuracy[,c(1,2,4)]

##                   Species Absolute_latitude  latitude
## 1 Craugastor_fleischmanni              31.4   9.93000
## 2     Dendrobates_auratus              26.0   9.94000
## 3    Liolaemus_multicolor              40.4 -24.30527
## 4 Sceloporus_occidentalis              45.8  33.95000