Typically, you would not want to merge your data with the map dataset. I would recommend merging it with mapsgfk.world_attr (which contains 1 observation per country) instead. ... View more

When using the new CT imported from the shapefile, it is already projected, there you will not need to run the "proc gproject" on it (unless you want some other projection, in which case you can use Jeff's suggestion to deproject and get the lat/long values, so you can re-project it using any projection you want).   Here is a simple example of importing the shapefile, and plotting the already-projected map:   proc mapimport datafile="/home/u59782200/data/CT_Planning_Regions/CT_Planning_Regions.shp" out=raw_ct; run; proc gmap data=raw_ct map=raw_ct; id ObjectID; choro planningRE;   Also, I notice that your title mentions maps.counties. Note that the 'maps' library is the old/legacy SAS maps (perhaps from 1996 vintage), and you'll probably want to use the latest 'mapsgfk' library instead (for your maps other than CT). Specifically mapsgfk.us_counties (the dataset label indicates this map is based on 2021 map data in the version of SAS I'm using).   ... View more

The goptions ctitle and ctext "should" be controlling your title and text, but I'm suspecting it is somehow coming out white on a white background instead.   Perhaps try specifying the color on the title statement itself, and see if that makes a difference, such as ...   TITLE1 color=red "test title"; ... View more

data have; length yymm $4; input YYMM Y :percent8.2; format y percent8.2; datalines; 2201 0.18% 2202 0.19% 2203 0.22% 2204 0.04% 2205 0.23% 2206 0.27% 2207 0.26% 2208 0.04% 2209 0.06% 2210 0.09% 2211 0.03% 2212 0.04% 2301 0.03% 2302 5.3% ; data want; set have; length date_string $8; date_string='20'||trim(left(yymm))||'15'; format sasdate yymmn4.; sasdate = input(date_string,yymmdd8.); proc print data=want; format sasdate date9.; axis1 label=none minor=none offset=(0,0) order=(0 to .1 by .01); axis2 label=('YYMM') offset=(2pct,2pct) order=('15jan2022'd to '15feb2023'd by month); symbol1 value=circle height=2pct color=blue; title1 h=5pct "Values, by Month"; proc gplot data=want; format y percentn7.0; plot y*sasdate/ autovref lvref=33 cvref=gray33 vaxis=axis1 haxis=axis2 ; ... View more

data have; input ID Greetings; datalines; 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 2 13 2 14 2 15 3 16 3 17 3 18 3 19 3 20 3 21 3 22 . 23 . 24 . 25 2 ; Proc format; Value Gre 1 = 'Hello' 2 = 'Hi' 3 = 'Dear' ; axis1 label=none minor=none order=(0 to 60 by 20) offset=(0,0); axis2 label=('Type of Greeting') offset=(9,9); pattern1 color=cyan; title1 h=6pct font='albany amt/bold' "Customized 3d Chart"; proc Gchart Data=have; Format Greetings Gre. ; vbar3d Greetings / discrete space=7pct width=5pct raxis=axis1 autoref cref=gray33 lref=33 maxis=axis2 type=percent outside=percent statfmt=comma5.0 ; ... View more

If you care about your colors, and want them to look consistent no matter what type of output you're creating (png, svg, java, activex, javaimg, actximg, pdf, etc) ... I suggest always using the hex rgb color codes. ... View more

Is there a way to get the lines to start at zero? ... View more

As Jeff says, Alaska & Hawaii have been moved/projected/scaled differently from the rest of the US, in this particular map (mapsgfk.us). Therefore, there is no one/simple command that will move/project/scale arbitrary lat/long coordinates and have them line up with the map.   But ... mapsgfk.uscity has applied the same move/project/scale to all the mapsgfk.uscity locations, so they line up with mapsgfk.us ... and therefore if you find the closest city to each of your lat/long locations, you could look up the x/y from mapsgfk.uscity, and plot those locations on the map. Here's a small example showing all the mapsgfk.uscity locations overlaid on mapsgfk.us:   title "mapsgfk.us with mapsgfk.uscity locations overlaid"; proc sgmap mapdata=mapsgfk.us noautolegend plotdata=mapsgfk.uscity; choromap / mapid=statecode; scatter x=x y=y / markerattrs=(size=1px color=red symbol=circlefilled) transparency=0; run;     ... View more

Why would you want any ticks on that axis at all? What purpose do they serve? ... View more

@tc  never disappoints! ... View more

Here's how to do it without combining/projecting/separating the map & data, by using the somewhat new feature of saving the projection parameters for the map and then applying them to the annotate dataset:   data zips; input c zip; cards; 1 65203 2 65211 3 65201 ; proc geocode data=zips out=geocodedzips (rename=(x=long y=lat)); run; data mocounties; set mapsgfk.us_counties (where=(statecode='MO')); run; proc gproject data=mocounties out=mocounties latlong eastlong degrees parmout=projparm; id county; run; proc gproject data=geocodedzips out=geocodedzips latlong eastlong degrees parmin=projparm parmentry=mocounties; id; run; data anno_locations; set geocodedzips; if c=1 then color="red"; if c=2 then color="bib"; if c=3 then color="viyg"; if c=4 then color="bippk"; xsys='2'; ysys='2'; hsys='3'; when='a'; function='pie'; rotate=360; size=.5; style='psolid'; output; length html $300; html='title='||quote(trim(left(c))||'0d'x||trim(left(put(zip,z5.)))); style='pempty'; color='deb'; output; run; pattern1 v=e; proc gmap map=mocounties data=mocounties anno=anno_locations; id county; choro segment / coutline=black levels=1 nolegend coutline=gray99; ... View more

Could you provide the code you used to produce your gchart pie chart? ... View more

Hypothetically, you could draw bars on a 2d map using annotate (where the bar starts at the center of each state, and the height of the bars is scaled based on the data). I'm not sure it would be a great way to present data, and not sure it would be worth the work. ... View more

Examples are a *lot* more useful, when they include the sample data! ... View more

The maps in the blog post show numeric (lat/long) axes - I've never seen Proc SGMAP do that. Are you 100% certain they were created with SGMAP?   ... View more