I'll retract my request for help. I managed to solve the problem by
calling map followed by a call to points after rescaling the cex. Its
the kind of solution that deserves a head-slap.
library(Hmisc)
library(maps)
sample.data <- data.frame(lat = c(12.1667, 14.6333, -6.874, 2.6167, 14.6833,
11.2, 3.2667, 11.4167, -13.8623, 13.1667), lon = c(-126.25, -103.4667,
-88.4572, -93.65, -97.7, -88.65, -111.0167, -119.7333, -100.0748,
-108.0333), exp.index = c(23.6266244576543, 2.06824648038330,
0, 1.46378849121688, 24.1824663424805, 0, 0.364600074527829,
4.468039274638, 0.543833744085446, 18.2301918845932))
xlim <- c(-150, -80)
ylim <- c(-5, 30)
cex.range <- c(1, 5)
exp.range <- range(sample.data$exp.index)
exp.cex <- diff(cex.range) * sample.data$exp.index / exp.range[2] +
cex.range[1]
map("world", fill = TRUE, col = "gray", xlim = xlim, ylim = ylim)
points(sample.data$lon, sample.data$lat, pch = 21, cex = exp.cex)
eric.archer wrote:
I'm using xYplot to create a bubble plot of values that I'd like to
visualize on top of a filled-in map of the coast, but I'm too much of
a lattice (which I understand xYplot is built on) and mapping newbie
to figure out how to begin to make this happen. Below is some sample
code that doesn't work but illustrates my goal. Any pointers anyone
has would be much appreciated. I'm using R v2.9.0 on Windows with
Hmisc v3.6-0 and maps v2.1-0.
Cheers,
eric
library(Hmisc)
library(maps)
sample.data <- data.frame(lat = c(12.1667, 14.6333, -6.874, 2.6167,
14.6833,
11.2, 3.2667, 11.4167, -13.8623, 13.1667), lon = c(-126.25, -103.4667,
-88.4572, -93.65, -97.7, -88.65, -111.0167, -119.7333, -100.0748,
-108.0333), exp.index = c(23.6266244576543, 2.06824648038330,
0, 1.46378849121688, 24.1824663424805, 0, 0.364600074527829,
4.468039274638, 0.543833744085446, 18.2301918845932))
xlim <- c(-150, -80)
ylim <- c(-5, 30)
print(xYplot(lat ~ lon, data = sample.data, size = sample.data$exp.index,
xlim = xlim, ylim = ylim
))
map("world", fill = TRUE, col = "gray", xlim = xlim, ylim = ylim, add
= TRUE)
--
Eric Archer, Ph.D.
Southwest Fisheries Science Center
National Marine Fisheries Service
3333 North Torrey Pines Court
La Jolla, CA 92037-1022
858-546-7121 (work)
858-546-7003 (FAX)
ETP Cetacean Assessment Program: http://swfsc.noaa.gov/prd-etp.aspx
Population ID Program: http://swfsc.noaa.gov/prd-popid.aspx
"Innocence about Science is the worst crime today."
- Sir Charles Percy Snow
"Lighthouses are more helpful than churches."
- Benjamin Franklin
"...but I'll take a GPS over either one."
- John C. "Craig" George
______________________________________________
R-help@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.
