1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
| library(ggplot2)
library(ggmap)
library(maptools)
library(maps)
library(scatterpie)
geodata<-read.table('R_Data_Map_Chrm3.tsv',sep="\t")
geodata2.0<-geodata[,-9]
names(geodata2.0)[1]<-"long"
names(geodata2.0)[2]<-"lat"
names(geodata2.0)[3]<-"antibiotic target replacement"
names(geodata2.0)[4]<-"antibiotic efflux"
names(geodata2.0)[5]<-"antibiotic inactivation"
names(geodata2.0)[6]<-"antibiotic efflux reduced permeability to antibiotic"
names(geodata2.0)[7]<-"antibiotic target protection"
names(geodata2.0)[8]<-"antibiotic target alteration"
ggplot(geo, aes(x = long, y = lat, group = group)) + geom_polygon(fill="lightgray", colour = "white") + geom_scatterpie(data = geodata2.0, aes(x=long, y=lat), cols = c("antibiotic target replacement","antibiotic efflux","antibiotic inactivation","antibiotic efflux reduced permeability to antibiotic","antibiotic target protection","antibiotic target alteration"),color=NA, alpha=1) |
Partager