Plot by FIPS Code using Shapefiles

library(rgdal)
library(colourschemes)
library(RColorBrewer)
library(maptools)
library(maps)
library(ggmap)
library(classInt)

max_change <- as.data.table(max_change)
max_change$FIPS <- sprintf("%05d",as.numeric(max_change$FIPS))

max_change <- max_change[with(max_change, order(FIPS)), ]

shapes=readOGR(dsn="F:/GCM//tl_2014_us_county", "tl_2014_us_county")
shapes$FIPS = paste(shapes$STATEFP,shapes$COUNTYFP,sep="")

m = match(shapes$GEOID,max_change$FIPS)

shapes$change = max_change$change[m]

sd = data.frame(col=colours,values=c(-2.5,-2,-1.5,-1,-.5,0))
sc = nearestScheme(sd)

plot(c(-129,-61),c(21,53),type="n",axes=FALSE,xlab="",ylab="")
title(main = "Mid-Century Projections (GISS-ER-2) using Max GAM")
plot(shapes,col=sc(shapes$change),add=TRUE,border="white",lwd=0.2,
colorkey=T)

brks.eq = classIntervals(shapes$change, style = "jenks")
spplot(shapes, "change",xlim = c(-129,-61), ylim = c(21,53),
at=brks.eq$brks,col.regions=pal, col="transparent",
       main = list(label="Mid-Century Projections (GISS-ER-2) using Max
GAM"))

Corey Sparks <corey.sparks <at> utsa.edu> writes:

Joining data the way you're doing it is dangerous, Roger Bivand and

others

describes a standard way to do this process here:

http://r-sig-geo.2731867.n2.nabble.com/Merging-shapefiles-and-csv-td7586839.html


Quite right - the chunks Corey is referring to are:

Please do refer to the vignette in the maptools package, and to previous
threads which have advised that merge() should not be used, and that the
row.names of the data frames be used as ID keys. Typically using match() on
the row.names of the two objects will show which are not correctly aligned.

and

Beware that the data from the objects may be jumbled - never use merge,
always use match() on the row.names vectors of the objects to ensure that
the key-IDs agree. Jumbled data happens, it is important not to think
"shapefile" but to think DBMS with the ID key your way of staying sane.

The maptools vignette is at:


http://cran.r-project.org/web/packages/maptools/vignettes/combine_maptools.pdf

or:

library(maptools)
vignette("combine_maptools")

Here I also suspect that you'll find that there are non-unique FIPS in the
county polygons file, so may need to go through
maptools::unionSpatialPolygons() first.

Roger

And I do an example using US Census data here, using merge():
http://spatialdemography.org/wp-content/uploads/2013/04/9.-Sparks.pdf
<http://spatialdemography.org/wp-content/uploads/2013/04/9.-Sparks.pdf>

look at page 134 of that pdf.

Hope this helps

-----
Corey Sparks, PhD
Assistant Professor
Department of Demography
University of Texas at San Antonio
501 West C?sar E. Ch?vez  Blvd
Monterey Building 2.270C
San Antonio, TX 78207
210-458-3166
corey.sparks 'at' utsa.edu
coreysparks.weebly.com
--
View this message in context:

http://r.789695.n4.nabble.com/Plot-by-FIPS-Code-using-Shapefiles-tp4706830p4706840.html

Sent from the R help mailing list archive at Nabble.com.