ggplot Fill SpatialPolygonsDataFrame












0















I am trying to plot the U.S. counties SPDF and shade the counties that appear in another vector (matching by GEOID). This is the code I currently have, but none of the counties are shaded.



I have downloaded the ShapeFile of U.S. county borders from this link.



library(ggplot2)
library(raster)
US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")
counties_in_radius <- c("12011", "12015", "12021", "12027", "12043", "12049", "12051", "12055", "12061")
base_plot <- ggplot(data = US_counties) +
coord_fixed(ratio = 1.3) +
geom_polygon(data=US_counties, aes(x = long, y = lat, group = group), size = 0.1, colour = "gray",
fill= (if(US_counties$GEOID %in% counties_in_radius) "blue" else NA))
base_plot


If I just do GEOID, I get the error object 'GEOID' not found, whereas when I do US_counties$GEOID, I get "the condition has length > 1 and only the first element will be used".



How can I properly fill the counties that are in counties_in_radius?










share|improve this question


















  • 1





    Fill needs to be assigned inside your aes, and you should do it with bare column names, no $ notation

    – camille
    Nov 28 '18 at 20:07
















0















I am trying to plot the U.S. counties SPDF and shade the counties that appear in another vector (matching by GEOID). This is the code I currently have, but none of the counties are shaded.



I have downloaded the ShapeFile of U.S. county borders from this link.



library(ggplot2)
library(raster)
US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")
counties_in_radius <- c("12011", "12015", "12021", "12027", "12043", "12049", "12051", "12055", "12061")
base_plot <- ggplot(data = US_counties) +
coord_fixed(ratio = 1.3) +
geom_polygon(data=US_counties, aes(x = long, y = lat, group = group), size = 0.1, colour = "gray",
fill= (if(US_counties$GEOID %in% counties_in_radius) "blue" else NA))
base_plot


If I just do GEOID, I get the error object 'GEOID' not found, whereas when I do US_counties$GEOID, I get "the condition has length > 1 and only the first element will be used".



How can I properly fill the counties that are in counties_in_radius?










share|improve this question


















  • 1





    Fill needs to be assigned inside your aes, and you should do it with bare column names, no $ notation

    – camille
    Nov 28 '18 at 20:07














0












0








0








I am trying to plot the U.S. counties SPDF and shade the counties that appear in another vector (matching by GEOID). This is the code I currently have, but none of the counties are shaded.



I have downloaded the ShapeFile of U.S. county borders from this link.



library(ggplot2)
library(raster)
US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")
counties_in_radius <- c("12011", "12015", "12021", "12027", "12043", "12049", "12051", "12055", "12061")
base_plot <- ggplot(data = US_counties) +
coord_fixed(ratio = 1.3) +
geom_polygon(data=US_counties, aes(x = long, y = lat, group = group), size = 0.1, colour = "gray",
fill= (if(US_counties$GEOID %in% counties_in_radius) "blue" else NA))
base_plot


If I just do GEOID, I get the error object 'GEOID' not found, whereas when I do US_counties$GEOID, I get "the condition has length > 1 and only the first element will be used".



How can I properly fill the counties that are in counties_in_radius?










share|improve this question














I am trying to plot the U.S. counties SPDF and shade the counties that appear in another vector (matching by GEOID). This is the code I currently have, but none of the counties are shaded.



I have downloaded the ShapeFile of U.S. county borders from this link.



library(ggplot2)
library(raster)
US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")
counties_in_radius <- c("12011", "12015", "12021", "12027", "12043", "12049", "12051", "12055", "12061")
base_plot <- ggplot(data = US_counties) +
coord_fixed(ratio = 1.3) +
geom_polygon(data=US_counties, aes(x = long, y = lat, group = group), size = 0.1, colour = "gray",
fill= (if(US_counties$GEOID %in% counties_in_radius) "blue" else NA))
base_plot


If I just do GEOID, I get the error object 'GEOID' not found, whereas when I do US_counties$GEOID, I get "the condition has length > 1 and only the first element will be used".



How can I properly fill the counties that are in counties_in_radius?







r ggplot2






share|improve this question













share|improve this question











share|improve this question




share|improve this question










asked Nov 28 '18 at 19:53









AlanAlan

8910




8910








  • 1





    Fill needs to be assigned inside your aes, and you should do it with bare column names, no $ notation

    – camille
    Nov 28 '18 at 20:07














  • 1





    Fill needs to be assigned inside your aes, and you should do it with bare column names, no $ notation

    – camille
    Nov 28 '18 at 20:07








1




1





Fill needs to be assigned inside your aes, and you should do it with bare column names, no $ notation

– camille
Nov 28 '18 at 20:07





Fill needs to be assigned inside your aes, and you should do it with bare column names, no $ notation

– camille
Nov 28 '18 at 20:07












0






active

oldest

votes












Your Answer






StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53527107%2fggplot-fill-spatialpolygonsdataframe%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Stack Overflow!


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53527107%2fggplot-fill-spatialpolygonsdataframe%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Contact image not getting when fetch all contact list from iPhone by CNContact

count number of partitions of a set with n elements into k subsets

A CLEAN and SIMPLE way to add appendices to Table of Contents and bookmarks