Example 2: Mapping and Annotating Values from the GINSIDE Procedure

goptions reset=global border;

data customer;
  length city $20;
  input lastname$ zip x y city $;

cards;
Smith    27611 1.374164 0.623436 Raleigh  
Jones    27560 1.375948 0.625278 Morrisville 
Doe      27513 1.375279 0.624922 Cary  
Patel    27520 1.369120 0.621970 Clayton  
White    27705 1.377910 0.628629 Durham
Short    27587 1.370373 0.627680 WakeForest
Phillips 27591 1.368124 0.624705 Wendell
Jackson  27597 1.367264 0.625629 Zebulon
;

data states;
  set maps.counties(where=(fipstate(state)="NC" and county=183));
run;

data combined;
  set customer states;
run;

proc gproject data=work.combined out=work.combined dupok;
  id state county;
run;

/*split the data*/
data work.states customer;
  set work.combined;
  if missing(zip) = 1 then output work.states;
  else output customer;
run;

proc ginside map=work.states data=customer out=mapout;
  id state county;
run;
/*see the resulting data*/
proc print;
run;

data points;
   set mapout;
   length function style color $ 8 position $ 1 text $ 20 ;
   retain xsys ysys "2" hsys "3" when "a" text 
"";
   retain rotate 360 style "solid" function "pie" position "5";
   color=
"black"; 
   size=2; 
   if missing(county) then color="red";
   output;
run;
title "Red dots are outside of Wake County";

proc gmap data=states map=states anno=points;
  id county;
  choro county / coutline=black nolegend;
run;
quit;