BookmarkSubscribeRSS Feed
altijani
Quartz | Level 8

I need to create a heat map for the US and its territories. I am using this code, which is giving me what I need for the "main land", but not the territories. I am using my own data for the values that create the heat map's color for each state.

 

Any help?

 

 

 

 

proc sort data=mapsgfk.us out=USBorders2 equals;

by statecode;

run;

proc sort data=sashelp.us_data(keep=density_2010 statecode) out=USDensity;

by statecode;

run;

data USDensity2;

merge USBorders2 USDensity;

by statecode;

if statecode ne '' then density_2010 = 100 ;

if statecode = 'IL' then density_2010 = - 20; 

if statecode = 'IN' then density_2010 = 88;

if statecode = 'IA' then density_2010 = 68;

if statecode = 'KS' then density_2010 = - 20; 

if statecode = 'KY' then density_2010 = 87;

if statecode = 'MI' then density_2010 = 41;

if statecode = 'MN' then density_2010 = - 20; 

if statecode = 'MO' then density_2010 = - 20; 

if statecode = 'NE' then density_2010 = - 20; 

if statecode = 'OH' then density_2010 = 46;

if statecode = 'WI' then density_2010 = - 20; 

if statecode = 'CT' then density_2010 = 90;

if statecode = 'DE' then density_2010 = 68;

if statecode = 'DC' then density_2010 = - 20; 

if statecode = 'ME' then density_2010 = - 20; 

if statecode = 'MD' then density_2010 = - 20; 

if statecode = 'MA' then density_2010 = 59; 

if statecode = 'NH' then density_2010 = - 20; 

if statecode = 'NJ' then density_2010 = 53; 

if statecode = 'NY' then density_2010 = 89;

if statecode = 'PA' then density_2010 = 62;

if statecode = 'RI' then density_2010 = - 20; 

if statecode = 'VT' then density_2010 = - 20; 

if statecode = 'AL' then density_2010 = 68;

if statecode = 'FL' then density_2010 = 49;

if statecode = 'GA' then density_2010 = - 20; 

if statecode = 'NC' then density_2010 = - 20; 

if statecode = 'SC' then density_2010 = - 20; 

if statecode = 'TN' then density_2010 = 87;
if statecode = 'VA' then density_2010 = - 20; 

if statecode = 'WV' then density_2010 = - 20; 

if statecode = 'AR' then density_2010 = 52; 

if statecode = 'CO' then density_2010 = 48;

if statecode = 'LA' then density_2010 = 27;

if statecode = 'MS' then density_2010 = 48;

if statecode = 'NM' then density_2010 = 77;

if statecode = 'OK' then density_2010 = 59; 

if statecode = 'TX' then density_2010 = 38;

if statecode = 'ND' then density_2010 = - 20; 

if statecode = 'AK' then density_2010 = 72;

if statecode = 'AS' then density_2010 = - 20; 

if statecode = 'AZ' then density_2010 = 38;

if statecode = 'CA' then density_2010 = - 20; 

if statecode = 'HI' then density_2010 = - 20; 

if statecode = 'ID' then density_2010 = 49; 

if statecode = 'MT' then density_2010 = 30;

if statecode = 'NV' then density_2010 = 43;

if statecode = 'OR' then density_2010 = 41;

if statecode = 'SD' then density_2010 = - 20; 

if statecode = 'UT' then density_2010 = 48; 

if statecode = 'WA' then density_2010 = 50; 

if statecode = 'WY' then density_2010 = 61;

if statecode = 'VI' then density_2010 = 49;

if statecode = 'PR' then density_2010 = 83;

run;

 

 

 

 

%annomac /* Make the SAS/GRAPH annotate macros available */

%centroid(mapsgfk.us,centers2,id) /* Polygon center is used to place each label */

data map2(drop=statecode segment);

set USDensity2 centers2(rename=(x=xCen y=yCen) in=a);

if a then Label=fipstate(input(substr(id,4,2), 2.)); /* State postal code */

id = catx('-', id, segment); /* Combine ID and Segment to make unique ID */

if label = 'ID' then ycen + -.025; /* Adjust a few label coordinates */

if label = 'MI' then ycen + -.025;

if label = 'HI' then ycen + -.01;

if label = 'NH' then ycen + -.012;

if label = 'VT' then ycen + .01;

if label = 'MD' then do; ycen + -.02; xcen + .037; end;

if label = 'AK' then ycen + .01;

if label = 'DE' then do; ycen + -.005; xcen + .025; end;

if label = 'CT' then do; ycen + -.009; xcen + .032; end;

if label = 'NJ' then do; xcen + .015; end;

if label = 'RI' then do; xcen + .015; end;

if label = 'DC' then delete;

run;

 

options orientation=landscape;

ods graphics on / height=7.6in width=12.8in;

ods pdf file='nv_highlow.pdf' notoc dpi=300;

 

proc sgplot data=map2 noautolegend;

title height=5 'Heat Map, 2019';

footnote height=3 "Notes, Notes, Notes" c=gray;

footnote height=3 "State Name" c=gray;

polygon x=x y=y id=id / outline lineattrs=(color=blue)

 

colorresponse=density_2010 dataSkin=matte



fill colormodel=(white red white green) name='map2';

scatter x=xcen y=ycen / markerchar=label markercharattrs=(size=10);

gradlegend 'map2';

label density_2010 = ' High Risk     Low Risk';

xaxis offsetmin=0.01 offsetmax=0.01 display=none;

yaxis offsetmin=0.01 offsetmax=0.01 display=none;

run;

ods graphics on / reset=all;
3 REPLIES 3
GraphGuy
Meteorite | Level 14

Which other 'territories' are you wanting, and where would you like for them to show up in the map?

 

You'll likely have to subset them from another map, and then re-scale them, and apply x/y offsets, to get them to show up 'around' the main US 48 states (like Hawaii and Alaska often are).

altijani
Quartz | Level 8
Just PR and VI for now.






GraphGuy
Meteorite | Level 14

mapsgfk.us_all contains Puerto Rico and Virgin Islands ...

 

proc gmap data=mapsgfk.us_all map=mapsgfk.us_all;
id id;
choro statecode;
run;

 

us_all_map.png

sas-innovate-2024.png

📢

ANNOUNCEMENT

The early bird rate has been extended! Register by March 18 for just $695 - $100 off the standard rate.

 

Check out the agenda and get ready for a jam-packed event featuring workshops, super demos, breakout sessions, roundtables, inspiring keynotes and incredible networking events. 

 

Register now!

How to Concatenate Values

Learn how use the CAT functions in SAS to join values from multiple variables into a single value.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 3 replies
  • 877 views
  • 1 like
  • 2 in conversation