Hi Rodrigo,
did you consider having a look at package "ggplot2" (thanks to Hadley
Wickham).
I guess you'll find something very close to your need (functions ggplot,
geom_point,...)
Hope this helps. regards. Olivier
Rodrigo Aluizio wrote:>
> 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]]
>
>
> ______________________________________________
> R-help at r-project.org 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.
>
>
--
View this message in context:
http://www.nabble.com/Bubble-Plot-Over-a-Map-tp23426635p23428390.html
Sent from the R help mailing list archive at Nabble.com.
