不使用渐变在R中制作合唱曲

不使用渐变在R中制作合唱曲,r,map,R,Map,我对R很陌生,我正试图用adm1(或州)分区绘制阿根廷蚂蚁物种的分布图 我已经从GADM网站下载了数据,我创建了一个csv文件,其中包含说明每个adm1中是否存在该物种的信息。 即使我没有梯度,我还能做一个合唱团吗?如果没有,我可以使用哪些其他类型的地图 我看过几个网站,包括,和Choropleth地图挑战赛,它们确实很有帮助,但它们都有数字数据,我使用的是当前(1)或不存在(0)列。我尝试过的不同软件包有sp(使用RColorBrewer)、ggplot2、rgeos、和maptools 以下

我对R很陌生,我正试图用adm1(或州)分区绘制阿根廷蚂蚁物种的分布图

我已经从GADM网站下载了数据,我创建了一个csv文件,其中包含说明每个adm1中是否存在该物种的信息。 即使我没有梯度,我还能做一个合唱团吗?如果没有,我可以使用哪些其他类型的地图

我看过几个网站,包括,和Choropleth地图挑战赛,它们确实很有帮助,但它们都有数字数据,我使用的是当前(1)或不存在(0)列。我尝试过的不同软件包有
sp
(使用
RColorBrewer
)、
ggplot2
rgeos
、和
maptools

以下是我目前掌握的代码:

library(sp)
library(RColorBrewer)
write.csv(atr, "atr_data.csv")
atr_data<-read.csv("atr_data.csv", header=TRUE)
  spcode country_code adm1_code newcol
1    atr          VEN     VE.AR      0
2    atr          PRY     PY.CE      0
3    atr          PAN     PA.CL      0
4    atr          PAN     PA.CL      0
5    atr          PAN     PA.PN      0
6    atr          PAN     PA.PN      0
库(sp)
图书馆(RColorBrewer)
write.csv(atr,“atr_data.csv”)

atr_数据在您提供的示例中出现了几个小问题

首先,可以使用
gadm@data$COLU no
而不是gadm$COLU no
。填写完状态/缺勤表后,通过直接访问插槽或使用
maptools
包中的
spCbind
可以将状态/缺勤数据添加到空间多边形数据框中

第二,如果您的
col_no
系数中只有2个级别,则必须将
mypalete
子集为2种颜色,因为Brewer调色板只能使用至少3个级别

library(sp)
library(RColorBrewer)
con <- url("http://gadm.org/data/rda/ARG_adm1.RData")
print(load(con))
close(con)

# Randomly assigning presence/absence data for display purposes only
gadm@data$col_no <- as.factor(rbinom(n = 24, size = 1, prob = 0.5))

myPalette <- brewer.pal(3, "Purples")

# col.regions is limited to 2 colors below with the middle color dropped.
spplot(gadm, zcol = "col_no", colorkey = TRUE, col.regions = myPalette[-2],
    main="Distribution of Atratus in Argentina")
库(sp)
图书馆(RColorBrewer)

con在您提供的示例中,出现了几个小问题

首先,可以使用
gadm@data$COLU no
而不是gadm$COLU no
。填写完状态/缺勤表后,通过直接访问插槽或使用
maptools
包中的
spCbind
可以将状态/缺勤数据添加到空间多边形数据框中

第二,如果您的
col_no
系数中只有2个级别,则必须将
mypalete
子集为2种颜色,因为Brewer调色板只能使用至少3个级别

library(sp)
library(RColorBrewer)
con <- url("http://gadm.org/data/rda/ARG_adm1.RData")
print(load(con))
close(con)

# Randomly assigning presence/absence data for display purposes only
gadm@data$col_no <- as.factor(rbinom(n = 24, size = 1, prob = 0.5))

myPalette <- brewer.pal(3, "Purples")

# col.regions is limited to 2 colors below with the middle color dropped.
spplot(gadm, zcol = "col_no", colorkey = TRUE, col.regions = myPalette[-2],
    main="Distribution of Atratus in Argentina")
库(sp)
图书馆(RColorBrewer)
骗局