Re: [R-sig-Geo] nc in Rnews_2005-2.pdf

This works for me. Perhaps you don't have 'sp' loaded? What does sessionInfo() show? Cheers, Mike. library(maptools) fname = system.file("shapes/sids.shp", package="maptools") p4s = CRS("+proj=longlat +datum=NAD27") nc = readShapePoly(fname, proj4string=p4s) plot(nc, axes = TRUE, col

[R-sig-Geo] nc in Rnews_2005-2.pdf

I read Rnews_2005-2.pdf and I can't do the test. For instance, on page 10: maptools and sp are active, > fname = system.file("shapes/sids.shp", package="maptools") > p4s = CRS("+proj=longlat +datum=NAD27") > nc = readShapePoly(fname, proj4string=p4s) > plot(nc, axes = TRUE, col=grey(1-nc$S