24898 - Create a map of selected counties within a state with PROC GMAP

Sample 24898: Create a map of selected counties within a state with PROC GMAP

This sample program uses PROC GMAP to create a map that contains the map boundaries for only the counties in the response data set. The STATE and COUNTY variables in the response data set contain the Federal Information Processing Standards (FIPS) codes for the state and counties being mapped.

The graphics output on the Results tab was produced using SAS^® 9.2. Submitting the sample code with releases of SAS prior to SAS 9.2 might produce different results.

 /* Set the graphics environment */
goptions reset=all cback=white border htitle=12pt htext=10pt;  

 /* Create a sample data set */
data products;
   input State County Count;
   datalines;
 17  19   3
 17  29   1
 17  31 124
 17  37   3
 17  43  22
 17  77   3
 17  89   2
 17  97   5
 17 113   6
 17 115   4
 17 119   1
 17 143   3
 17 161   4
 17 167   8
 17 201   4
 18   3   4
 18   5   3
 18  35   1
 18  39   4
 18  63   3
 18  67   3
 18  85   2
 18  89   5
 18  97  32
 18 127   1
 18 141   4
 18 157   2
 18 163   2
 18 167   1
 26  21   1
 26  49   2
 26  65   1
 26  75   5
 26  77   8
 26  81   4
 26  83   1
 26  99  11
 26 101   1
 26 103   1
 26 111   4
 26 125  12
 26 139   4
 26 145   6
 26 155   8
 26 161   5
 26 163  26
 39   9   6
 39  17   6
 39  35  34
 39  45   3
 39  49  28
 39  57   2
 39  61  21
 39  63   3
 39  83   1
 39  89   2
 39  93   1
 39  95   9
 39  99   4
 39 109   2
 39 113  14
 39 149   2
 39 151   3
 39 153   9
 39 159   1
 39 173   1
 55   9   2
 55  25  12
 55  63   1
 55  79  25
 55 101   4
 55 117   2
 55 127   2
 55 139   2
;
run;
 
 /* Create a data set that contains the map boundaries */
 /* for only the counties in the response data set     */
 /* also create a data set for only the states which   */
 /* are contained in the response data                 */
data counties(drop=count test) selstate(keep=state);
   retain test 0 ;
   merge products(in=inprod) maps.counties(in=inmap);
   by state county ;
   if inprod and inmap;

   /* Determine which states are included in the response data */
   if state ^= test then do;
      output selstate;
      test=state;
   end;
   output counties;
 run;

 /* Get the unprojected coordinates for the states */
data states;
   merge maps.states selstate(in=inprod);
   by state;
   if inprod;
run;

 /* Concatenate states with counties map data set */
 /* for projection                                */
data newmap;
   set states counties;
run;

 /* Project the map coordinates and create a new map data set */
proc gproject data=newmap out=newmap;
   id state county;
run;
quit;

 /* Specify distribution format */
proc format;
   value dist 1-5='1-5' 6-10='6-10' 11-high='over 10';
run;

title 'Products Installed in the Great Lakes Region';
title2 h=10pt 'By County';

 /* Generate the map */
proc gmap data=products map=newmap all;
   id state county;
   choro count / discrete;
   format count dist.;
   label count='# of Products';
run;
quit;

These sample files and code examples are provided by SAS Institute Inc. "as is" without warranty of any kind, either express or implied, including but not limited to the implied warranties of merchantability and fitness for a particular purpose. Recipients acknowledge and agree that SAS Institute shall not be liable for any damages whatsoever arising out of their use of this material. In addition, SAS Institute will provide no support for the materials contained herein.

Type:	Sample
Topic:	SAS Reference ==> Procedures ==> GMAP

Date Modified:	2005-08-24 16:06:31
Date Created:	2004-11-11 11:07:57

Product Family	Product	Host	SAS Release
			Starting	Ending
SAS System	SAS/GRAPH	All	n/a	n/a

Support

Sample 24898: Create a map of selected counties within a state with PROC GMAP

Operating System and Release Information