Hi everybody,
Well here is another doubt! Im elaborating some maps, and I pretend to plot
over it the abundance of species. I could do that using points(), but Id
like to use bubble() from sp package.
When trying to use bubble() I face the fallowing issue: When I call the
bubble function it overwrite everything I previously plotted. It doesnt
happen with points (obviously), only with bubble(). I tried to specify
panel.smooth as panel function to bubble, I tried to use add=T as argument,
but it didnt accept such thing.
The only thing that I thought right now, so I still have to try, is saving
the bubble() as an object and try to manipulate its properties over my map.
I already searched the list posts but wasnt able to find a similar doubt
with its possible solutions.
Bellow is the map code using points, Im trying to substitute the points()
functions with one bubble() function.
for (i in 1:length(colnames(Spp10ml))) {
png(paste(colnames(Spp10ml),'.png',sep='')[i],width=5.2,height=3.6,pointsize
=5,units='in',bg=NA,restoreConsole=T,res=1200)
par(mar=c(5,0,0,0))
map('worldHires','brazil',ylim=c(-5.15,-4.55),xlim=c(-37,-36.1),type='n')
rect(-37.1,-5.25,-36,-4.9,density=NULL,angle=45,col='#dbf4ff',border=F)
plot(Batimetria,ylab='',xlab='',border=F,col=Cor$Cor,add=T)
plot(Municipios,ylab='',xlab='',col=rgb(245,232,140,max=255,alpha=255),add=T
,lwd=0.2)
axis(1,xaxp=c(-37,-36.1,5),cex.axis=1.5)
axis(2,yaxp=c(-5.15,-4.55,5),cex.axis=1.5)
text(Prof$Long,Prof$Lat,rownames(Prof),col='black',cex=1.2,srt=-24,family='A
R',font=2)
text(c(-36.96,-36.96,-36.865,-36.54,-36.37,-36.2),c(-4.96,-5.05,-5.11,-5.14,
-5.14,-5.14),c('Areia
Branca','Serra
do Mel','Porto do
Mangue','Macau','Guamaré','Galinhos'),col='Black',cex=1.3,font=2)
rect(-36.20,-4.56,-36.155,-4.565,density=NULL,angle=45,col='black',border='b
lack',lwd=0.1)
rect(-36.155,-4.56,-36.111,-4.565,density=NULL,angle=45,col='white',border='
black',lwd=0.1)
shadowtext(c(-36.20,-36.155,-36.111,-36.215),c(-4.573,-4.573,-4.573,-4.562),
c('0','5','10','km'),col='white',cex=1)
Norte(-36.155,-4.60,0.02)
box(which='plot',lty='solid')
points(Pontos$long[Spp10ml[,i]>0 & Spp10ml[,i]<=50],Pontos$lat[Spp10ml[,i]>0
& Spp10ml[,i]<=50],pch=21,bg=rgb(70,70,70,max=255,alpha=255),col=
'black',cex=1.5)
points(Pontos$long[Spp10ml[,i]>=51 &
Spp10ml[,i]<=200],Pontos$lat[Spp10ml[,i]>=51 &
Spp10ml[,i]<=200],pch=23,bg='yellow',col='black',cex=1.5)
points(Pontos$long[Spp10ml[,i]>=201 &
Spp10ml[,i]<=400],Pontos$lat[Spp10ml[,i]>=201 &
Spp10ml[,i]<=400],pch=22,bg='orange',col='black',cex=1.5)
points(Pontos$long[Spp10ml[,i]>=401],Pontos$lat[Spp10ml[,i]>=401],pch=24,bg=
'red',col='black',cex=1.5)
dev.off()}
Thanks in advance.
-------------------------------------------------------------
MSc. <mailto:[email protected]> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - Brasil
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
