If I understand your goal, these are the steps I think you need to do: import zipcode shapefiles into SAS using PROC MAPIMPORT After import, you will convert  these lat / long coordinates to x/y using atan using code.  IMPORTANT – for GMAP, the longitude coordinates for the western hemisphere need to be positive.  So you may need to use abs() to convert long from negative to positive.      2. merge zipcode shapefile with dataset that identifies zipcodes with regions If there are zipcodes that are not associated with special regions, you will still need to assign them unique region numbers.  Otherwise the unassigned zipcodes will be reduced into one mega-region in the PROC GREDUCE.      3. do PROC GREDUCE based on region      4. do PROG GPROJECT      5. create choropleth map with PROC GMAP Hope this helps. ... View more

Eric, Go to Lex Jansen's wonderful website - www.lexjansen.com  and search for the SUGI paper "Constructing ROC Curves with the SAS System" by Andrew Stead.  The paper contains a nice macro for calculating and plotting ROC curves. Renee ... View more

Try  this: drop the first symbol statment and change SYMBOL2  to SYMBOL1 so there is only one symbol statement then in the PLOT statement drop the overlay and just plot the lower and upper  - PLOT LOW_UP2 * VISN; this should plot the lower and upper bound and the HILOTJ option will calculate the mean of the lower and upper and join it for each timepoint Another option is to reshape the data so each time point has 3 observations - low, mean, high.  Then plot that with one symbol statemement using HILOCTJ. Read the documentation on HILO to understand all of the options. Hope all that is clear.  For some reason I can't paste SAS code here, otherwise I would paste a snippit for clarity. Renee ... View more