Building on "Proc GMap County and State Boundries"

CathyB — Mon, 09 Sep 2019 22:56:48 GMT

Using the excellent instructions recommended in Proc GMap County and State Boundaries" responses, I created a map with heavy state boundaries and lighter county boundaries.

Then, I used the instructions found in this document https://www.lexjansen.com/nesug/nesug03/bt/bt005.pdf to create a colored chrono map with select counties in various colors.

However, the results are a mish mash of dark and lighter lines. The appropriate counties are correctly colored.

proc gmap map=AdaMap data=AdaMap anno=State_outline;
id state county;
choro segment / levels=1
nolegend
coutline=grayaa; *Light grey county lines;
run;

* Create st/co fips var in Map dataset and data dataset for matching;
data AdaData;
set AdaData;
FIPS=cats(W_StFIPS,"_",W_CoFIPS);
put FIPS=;
run;

data AdaMap;
set adamap;
FIPS=cats(State,"_",County);
put FIPS=;
run;

/* Preparing for the chrono map */

proc format;
value Workers_fmt
0-25 = '1 to 25'
26-50 = '26 to 50'
51-100 = '51 to 100'
101-200 = '100 to 200'
201-500 = '201 to 500'
501-1000 = '501-1,000'
1000-500000 = 'Over 1,000'
;
run;

* fill patterns for the map areas;
pattern1 v=s c=pab;
pattern2 v=s c=grpb;
pattern3 v=s c=mogb;
pattern4 v=s c=mopb;
pattern5 v=s c=grb;
pattern6 v=s c=papb;

* create instructions for legend location and appearance;
legend1
origin=(5,10) pct
mode=share
across=1
label=none
value=(j=r)
shape=bar(3,4) pct;

proc gmap
data=AdaData
map=AdaMap (where=(density lt 6))
all;
id FIPS;
choro workers / discrete
legend=legand1;
format workers workers_fmt.;
note ' Ada Map Data set'
j=l ' Ada County Outbound Commuters';
run;
quit;

Here is my code. (Previous code creates the state/county map, and the following creates the colored chrono map.)

Where am I going wrong? Is there an order of procs I'm overlooking that will annotate the state outlines on top of the colored choro map?

......

proc gmap map=AdaMap data=AdaMap anno=State_outline;
id state county;
choro segment / levels=1
nolegend
coutline=grayaa; *Light grey county lines;
run;

* Create st/co fips var in Map dataset and data dataset for matching;
data AdaData;
set AdaData;
FIPS=cats(W_StFIPS,"_",W_CoFIPS);
put FIPS=;
run;

data AdaMap;
set adamap;
FIPS=cats(State,"_",County);
put FIPS=;
run;

/* Preparing for the chrono map */

proc format;
value Workers_fmt
0-25 = '1 to 25'
26-50 = '26 to 50'
51-100 = '51 to 100'
101-200 = '100 to 200'
201-500 = '201 to 500'
501-1000 = '501-1,000'
1000-500000 = 'Over 1,000'
;
run;

* fill patterns for the map areas;
pattern1 v=s c=pab;
pattern2 v=s c=grpb;
pattern3 v=s c=mogb;
pattern4 v=s c=mopb;
pattern5 v=s c=grb;
pattern6 v=s c=papb;

* create instructions for legend location and appearance;
legend1
origin=(5,10) pct
mode=share
across=1
label=none
value=(j=r)
shape=bar(3,4) pct;

My ultimate goal will be to superimpose directional arrows from the source county to the commuted-to counties. If anyone has done something similar, I would like to know how it was done.

Thanks to all.

Cathy

Re: Building on "Proc GMap County and State Boundries"

GraphGuy — Wed, 11 Sep 2019 11:25:36 GMT

This example might help some:

http://robslink.com/SAS/democd102/us_dry_counties.htm

http://robslink.com/SAS/democd102/us_dry_counties_info.htm

topic Re: Building on "Proc GMap County and State Boundries" in Graphics Programming

Building on "Proc GMap County and State Boundries"

Re: Building on "Proc GMap County and State Boundries"