Bubble Plot Over a Map

Hi everybody,

Well here is another doubt! I?m elaborating some maps, and I pretend to
plot
over it the abundance of species. I could do that using points(), but I?d
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 doesn?t
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 didn?t 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 wasn?t able to find a similar doubt
with its possible solutions.

 

Bellow is the map code using points, I?m 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:r.aluizio at gmail.com> 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]]