[R] map at fips level using multiple variables
jones at reed.edu
Thu Dec 8 19:11:35 CET 2011
since each county gets a single color, you would need to combine two
color schemes. I am skeptical that this will work well, but you could
try using rgb(). For example - code one variable using the red
component, the other using the blue component.
On Wed, Dec 07, 2011 at 11:14:27PM -0500, bby2103 at columbia.edu wrote:
> Hi David,
> Sorry it sounds vague.
> Here is my current code, which gives the distribution of family size
> at US county level. You will see on a US map family size
> distribution represented by different colors.
> Now if i have another variable income, which has 3 categories(<50k,
> 50k-80k,>80k). How do I show 5x3 categories on the map? I guess I
> could always create a third variable to do this. Just wondering
> maybe there is a function to do this readily?
> Thank you!
> Bonnie Yuan
> hist(y,nclass=15) # histogram of y
> # specify the cut-off point of y
> y.colorBuckets=as.numeric(cut(y, c(1,2, 3,6)))
> # legend showing the cut-off points.
> map("county", col = colors[colorsmatched], fill = TRUE, resolution = 0)
> map("state", col = "white", fill = FALSE, add = TRUE, lty = 1, lwd = 0.2)
> title("Family Size")
> legend("bottom", legend.txt, horiz = TRUE, fill = colors, cex=0.7)
> Quoting David Winsemius <dwinsemius at comcast.net>:
> >On Dec 7, 2011, at 6:12 PM, bby2103 at columbia.edu wrote:
> >>Hi, I just started playing with county FIPS feature in maps
> >>package which allows geospatial visualization of variables on
> >>US county level. Pretty cool.
> >Got code?
> >>I did some search but couldn't find answer to this question--how
> >>can I map more than 2 variables on US map?
> >"2 variables" is a bit on the vague side for programming purposes.
> >>For example, you can map by the breakdown of income or family
> >>size. How do you further breakdown based on the values of both
> >>variables and show them on the county FIPS level?
> >"breakdown" suggests a factor construct. If so, then :
> >But the "show" part of the question remains very vague.
> > Can't you be a bit more specific? What DO you want?
> >David Winsemius, MD
> >West Hartford, CT
> R-help at r-project.org mailing list
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
jones at reed.edu
More information about the R-help