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]] ______________________________________________ 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.
