Hi, using a portion of your data (you could drop S and E) ...   data x; input (s e) (:$20.); if findc(s,'-') then s=catt(s,'-15'); if findc(e,'-') then e=catt(e,'-15'); start = input(s,anydtdte10.); end   = input(e,anydtdte10.); format st: end mmddyy10.; datalines; 6/23/2015 7/20/2015 7/14/2015 9/5/2015 8/4/2015 TBD 8/20/2015 11/2/2015 8/27/2015 TBD 9-Feb 1-Jul 15-Jan 1-Jul 15-Jan 1-Jul 3-Feb 1-Jul 9-Feb 1-Jul 1-Jul TBD 1-Jul TBD 8/27/2015 TBD 10/6/2015 TBD 10/13/2015 11/2/2015 ;   data set X ... Obs        s         e                 start           end   1    6/23/2015     7/20/2015    06/23/2015    07/20/2015   2    7/14/2015     9/5/2015     07/14/2015    09/05/2015   3    8/4/2015      TBD          08/04/2015             .   4    8/20/2015     11/2/2015    08/20/2015    11/02/2015   5    8/27/2015     TBD          08/27/2015             .   6    9-Feb-15      1-Jul-15     02/09/2015    07/01/2015   7    15-Jan-15     1-Jul-15     01/15/2015    07/01/2015   8    15-Jan-15     1-Jul-15     01/15/2015    07/01/2015   9    3-Feb-15      1-Jul-15     02/03/2015    07/01/2015  10    9-Feb-15      1-Jul-15     02/09/2015    07/01/2015  11    1-Jul-15      TBD          07/01/2015             .  12    1-Jul-15      TBD          07/01/2015             .  13    8/27/2015     TBD          08/27/2015             .  14    10/6/2015     TBD          10/06/2015             .  15    10/13/2015    11/2/2015    10/13/2015    11/02/2015     ... View more

Hi.  Another approach is to treat the XLS file as a data set.  Use a LIBNAME statement and use the SHEET name within the XLS file as shiown in the data step   ...   libname x 'z:\test.xls';   data gas_el (keep=dt days eld delc gpd dgas); set x.'sheet1$'n; eld = kwh/days; delc = elec/days; gpd = therms/days; dgas = gas/days; label eld = 'KWHS/DAY' gpd = 'THERMS/DAY' delc = 'ELEC $/DAY' dgas = 'GAS $/DAY' dt   = 'DATE' days = 'DAYS USE' ; format dt monyy. gas elec delc dgas dollar7.2 eld gpd 5.1; run;   libname x clear;   data set GAS_EL ...                  DAYS                   ELEC    THERMS/        GAS Obs    DATE      USE    KWHS/DAY      $/DAY      DAY        $/DAY    1    DEC14     34       13.6        $2.45       8.1       $4.95   2    NOV14     29       11.2        $1.77       5.6       $3.65   3    OCT14     32       12.6        $2.16       2.2       $2.04   4    SEP14     30       13.1        $2.30       0.9       $1.29   5    AUG14     28       15.4        $2.67       0.6       $1.21   6    JUL14     33       14.8        $2.53       0.5       $1.08   7    JUN14     27       15.0        $2.58       0.7       $1.45   8    MAY14     32       13.3        $2.49       1.5       $2.19   9    APR14     32       11.9        $1.76       4.4       $4.45  10    MAR14     29       11.0        $3.11       9.9       $9.71  11    FEB14     29       11.4        $2.67      10.3       $7.57  12    JAN14     29       12.6        $2.87      11.9       $7.81 ... View more

Hi. Here's SAS code that uses tow annotate data sets to show the zip centroids plus your geocoded data.  There are not many comments so ask questions if there's code you do not understand.    I have one comment ... rather than use the map data sets in the MAPS library, start uing the maps in the MAPSGFK library.  As you can see, if you use MAPSGFK.US_COUNT IES, you no longer have to go through the older method of combining annnotate data with the map/project the combined data/separate the files.  Also, you do not have to convert your annotate data lat/long in degrees to radians (HOORAY).   Here's the code and the maps are attached. * read Harford County data ... create annotate data set ; libname hc 'z:\temp\geocoded_data.xls'; data hc; length color style $20; retain xsys ysys '2' hsys '3' function 'label' color 'green' size 2        style 'wingdings' text '6c'x when 'a' state 24; do until (last);    set hc.'geocoded_data$'n (keep=x y _score_) end=last;    output; end; * add location of Harford CC (from Google Maps); y = 39.5583912; x =-76.2842196; style = 'wingdings 3'; text  = '70'x; size  = 6; output; run;   libname hc clear;   * project data set hc; proc gproject data=hc out=anno_hc dupok   parmin=mapsgfk.projparm parmentry=us_counties;   id state; run;   * find lat/long of zips in Harford County ... create annotate data set; data centroids; retain xsys ysys '2' hsys '3' function 'label' color 'yellow' size 3        style 'wingdings' text 'ab'x when 'a'; set sashelp.zipcode (keep=state county x y); where state eq 24 and county eq 25; run;   * project data set centroids; proc gproject data=centroids out=anno_zip dupok   parmin=mapsgfk.projparm parmentry=us_counties;   id state; run;   * graphics options ... output to a PNG file; goptions reset=all gunit=pct ftext='calibri' dev=png gsfname=gout xpixels=1500 ypixels=1000;   filename gout 'z:\harford.png';   * draw a map ... two annotate data sets used; proc gmap data=mapsgfk.us_counties map=mapsgfk.us_counties anno=anno_zip; where state eq 24 and county eq 25; id county; choro county / nolegend statistic=first anno=anno_hc; run; quit;   * add a LEGEND to describe the points; legend1 origin=(3,15)pct value=none shape=bar(.00001,.00001)pct mode=share label=(h=2.5        j=l f='wingdings 3' '70'x f='calibri' ' HARFORD CC'        j=l f='wingdings'   'ab'x f='calibri' ' ZIP CENTROIDS'        j=l f='wingdings'   '6c'x f='calibri' ' GEOCODED (GREEN SCORE 50+, RED SCORE <50)');     filename gout 'z:\harford_leg.png';   * change annotate colors based on gecoding score; data anno_hc; set anno_hc; select;    when (_score_ lt 50) color='red';    when (_score_ ge 50) color='green';    otherwise; end; run;   * draw a map ... two annotate data sets used; proc gmap data=mapsgfk.us_counties map=mapsgfk.us_counties anno=anno_zip; where state eq 24 and county eq 25; id county; choro county / statistic=first legend=legend1 anno=anno_hc; run; quit; ... View more