Bubbleplot due to Finnish method
This function plots multivariate data with respect to the value. The size of the bubble represents the value of the datapoint.
bubbleFIN(x, y, z, radi = 10000, S = 9, s = 0.9, wa = 0, wb = 0.95, wc = 0.05, plottitle = "BubblePlot", legendtitle = "Legend", text.cex = 1, legtitle.cex = 1, backgr = "kola.background", leg = TRUE, ndigits = 1)
| x | x coordinates | 
| y | y coordinates | 
| z | measured value at point (x,y) | 
| radi | scaling for the map | 
| S, s | control the size of the largest and smallest bubbles | 
| wa, wb, wc | factors which defines the shape of the exponential function | 
| plottitle | the titel of the plot | 
| legendtitle | the titel of the legend | 
| text.cex | multiplier for the size of the labels | 
| legtitle.cex | multiplier for the size of the legendtitle | 
| backgr | which background should be used | 
| leg | if TRUE the bubbles are plotted to the legend | 
| ndigits | how much digits should be plotted at the legend | 
The smallest bubbles represent the 10% quantile and the biggest bubbles represent the 99
Plots bubbles in the existing plot.
Peter Filzmoser <P.Filzmoser@tuwien.ac.at> http://cstat.tuwien.ac.at/filz/
C. Reimann, P. Filzmoser, R.G. Garrett, and R. Dutter: Statistical Data Analysis Explained. Applied Environmental Statistics with R. John Wiley and Sons, Chichester, 2008.
data(kola.background)
data(ohorizon)
el=ohorizon[,"Mg"]
X=ohorizon[,"XCOO"]
Y=ohorizon[,"YCOO"]
plot(X,Y,frame.plot=FALSE,xaxt="n",yaxt="n",xlab="",ylab="",type="n") #plot bubbles with background
plotbg(map.col=c("gray","gray","gray","gray"),add.plot=TRUE)
bubbleFIN(X,Y,el,S=9,s=2,plottitle="",legendtitle="Mg [mg/kg]", text.cex=0.63,legtitle.cex=0.80)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.