Hi again, GraphGuy: Sorry to bother you again, but I can't figure out how to have the dots on the map reflect presence/absence information (0=absence, 1=presence). All my attempts have failed thus far. The data and code are below. j -------------------- data coordinates; input lat long present; datalines; 28.9787 -95.1672 0 29.1468 -94.8645 0 29.2313 -89.8375 1 29.4695 -94.0852 0 29.5513 -92.899 1 29.5282 -93.9293 1 29.016 -95.098 1 29.8865 -88.1793 1 29.5088 -93.9617 0 29.4177 -92.587 0 29.1007 -94.02 1 29.5347 -93.5135 0 25.2593 -81.9585 1 28.9765 -91.7018 1 28.1155 -93.6067 1 28.3387 -94.4748 0 28.9765 -91.7018 1 ; run; data us48; set mapsgfk.us_states; if state not in (2 15 72); run; proc gproject data=us48 out=gulf latlon degrees eastlong longmin=-98 longmax=-68 latmin=25 latmax=36.5 parmout=projparm; where density<5; id state; run; /* project the coordinates using the same projection parameters as were used to project the map */ proc gproject data=coordinates out=coordinates latlon degrees eastlong longmin=-98 longmax=-68 latmin=25 latmax=36.5 parmin=projparm parmentry=gulf; id; run; /* convert the projected coordinates into annotate dataset */ data coordinates; set coordinates; xsys='2'; ysys='2'; hsys='3'; when='a'; function='pie'; rotate=360; size=1.0; color='red'; run; title "US Southern Atlantic and Gulf of Mexico"; pattern value=mempty color=blue; proc gmap map=gulf data=gulf all anno=coordinates; id state; choro state / statistic=first nolegend levels=1; run;
... View more