
Dear all,

I am  new to R and even newer to mapping with R. I have installed 
PBSMapping, maptools and maps, and I am struggling to produce a map of 
fish abundance in the Grand Banks. I have tried several approaches to 
plot the bathymetry in my study area:

1) download the data from http://topex.ucsd.edu/cgi-bin/get_data.cgi  
(42N46N313W307W)
     To start with, the map I get is blank between lon -52W -53W. 
Besides, it seems that the shallower isobath is 500 m. Or maybe I am 
wrong here?

2) import a shapefile with the bathymetry, but I get the following error:

     test <- read.shapefile("C:/data")
               Error in seq.default(1, (num.points * 2), by = 2) :   
wrong sign in 'by' argument


3) use a Rdata version of the above mentioned shapefile. This is the map 
I like the most, but when I try to plot the study area (which is >200 nm 
from the nearest shore), I get the error

               map('worldHires',xlim=c(-53,-47),ylim=c(42,46))
               Error in map.poly(database, regions, exact, xlim, ylim, 
boundary, interior,  :   nothing to draw: all regions out of bounds

The next problem will be plotting fish abundance. I was thinking of 
contours or scaled bubbles, but I am not sure I can do any of those with 
PBSMapping. Am I right? Any advice?

Thanks for your help,

Elena



	[[alternative HTML version deleted]]

