Excellent suggestions Professors! I really appreciated it. This is what I am using:
library(data.table) > library(rgdal) > library(colourschemes) > library(RColorBrewer) > library(maptools) > library(maps) > library(ggmap) > library(classInt) max_change is my csv file while shapes is my spatial data: max_change <-read.csv ("F:/GCM/max_change.csv") > 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] I am plotting the results now. I am getting different plots regarding county boundaries (most likely because I am doing something wrong). When I use the *plot *command: colours = brewer.pal(6,"PuRd") > 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) I get a "beautiful" plot with the county boundaries clearly separated but for *spplot*; pal = brewer.pal(6,"Reds") > 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")) The county boundaries/FIPS are not defined. What am I doing wrong? Thanks again! Sincerely, Shouro On Wed, May 6, 2015 at 12:23 PM, Roger Bivand <roger.biv...@nhh.no> wrote: > 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. > > > > ______________________________________________ > > R-help <at> r-project.org mailing list -- To UNSUBSCRIBE and more, see > > 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. > > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. > -- Shouro Dasgupta PhD Candidate | Department of Economics Ca' Foscari University of Venezia ------------------------------ Junior Researcher | Fondazione Eni Enrico Mattei (FEEM) Isola di San Giorgio Maggiore, 8 | 30124 Venice, Italy Tel: +39 041 2700 436 [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.