MWSUG Conference Locations Used SAS/GRAPH in 2011 to knock out an image for a Zazzle t-shirt (production run of 1!) commemorating the 36th anniversary of SAS Global Forum.   So, with MWSUG 2019 kicking off next Sunday in Chicago, here's a quick SAS ODS Graphics take on a chart to mark that conference's 30th anniversary.   Probably could've used the newer SGMAP procedure, but after seeing tips in write-ups by Sanjay Matange, Rick Wicklin, Robert Allison, and Louise Hadden, decided to opt for PROC SGPLOT and polygon/scatter plots!   Zazzle T-Shirt Mockup SGPANEL Detail SAS CODE * Fun with SAS ODS Grapics: Maps of MWSUG Conference Locations (1990-2019); data mwsugCities(keep=year city st); * Get locations of past conferences from MWSUG website; filename mwsug url "http://mwsug.org/pastconf.html"; * Web page with data on past conferences; infile mwsug truncover; input @"<tr><td colspan=2><b>" year 4. +3 cityst $char255.; * Scrape year, city, state; if year^=.; city=tranwrd(scan(cityst,1,","),"St. Louis","Saint Louis"); * Get city, recode St. Louis; st=scan(scan(cityst,1,'<'),-1,",< "); * Get state codes or names; do s=1 to 56; * Translate state names into codes; if st=fipnamel(s) then st=fipstate(s); end; proc sql; * Add conference location for 2019; insert into mwsugCities values(2019,"Chicago","IL"); * Map cities to colors (Tableau 10 Light RGB codes from gist.github.com/chebee7i/4041f848ee2710d500e6); create table Hosted as select city, count(*) as hosted from mwsugCities group by 1 order by 2 desc, 1; create table CityAttrMap as /* Create ODS Graphics attribute map of cities to colors for charts */ select "city" as id, city as value, /* 13 colors for 13 unique cities */ scan("CXAEC7E8 CXFFBB78 CX98DF8A CXFF9896 CXC5B0D5 CXC49C94 CXF7B6D2 CXC7C7C7 CXDBDB8D CX9EDAE5 TAN GOLD TURQUOISE",monotonic()," ") as fillcolor from hosted; * Add lat/long to conference locations; create table mwsugCitiesGeo(rename=(x=xC y=yC)) as select m.year, u.* from mwsugCities m, sasfiles.uscity u where m.city=u.city and m.st=u.statecode order by m.year; libname mymaps '/folders/myfolders/mymaps/'; * Get outlines of MWSUG states from SAS US map; data mwsugStates; set mymaps.us(where=(statecode in ('IL' 'IN' 'IA' 'KS' 'MI' 'MN' 'MO' 'NE' 'ND' 'OH' 'SD' 'WI'))); seq=_n_+1; * Maintain original order; length pid $8; keep state segment x y pid seq; pid=put(state, 3.0)||put(segment, 3.0); * Make unique polygon ID for each state+segment combo; proc sql; * Duplicate state outline data for each year and add year/city; create table mwsugStates2 as select * from mwsugStates, mwsugCitiesGeo(keep=year City) order by seq; data mwsug; * Merge state outlines and city points together; set mwsugStates2 mwsugCitiesGeo; proc sql noprint; * Calc aspect ratio of points in map; select (max(y)-min(y))/(max(x)-min(x)) into :aspect from mwsug; * Calc from: blogs.sas.com/content/sastraining/2017/03/21/create-pretty-map-proc-sgplot/; ods listing image_dpi=250 gpath='/folders/myfolders/MWSUG2019/images'; * Image setup parameters; ods graphics on / reset antialias width=14in height=17in imagename="MWSUG" noborder ANTIALIASMAX=8400 LABELMAX=8400; proc sgpanel data=mwsug aspect=&aspect subpixel noautolegend dattrmap=CityAttrMap; * Draw panel of maps!; title height=54pt "MWSUG@30"; * 30 years of MWSUG conferences!; panelby year / rows=6 columns=5 onepanel uniscale=all noheaderborder headerbackcolor=white novarname headerattrs=(weight=bold size=20pt color=black) noborder spacing=5; rowaxis display=none offsetmin=.0 offsetmax=.0; colaxis display=none offsetmin=.0 offsetmax=.0; polygon x=x y=y id=pid / fill outline group=city attrid=city lineattrs=(color=white); * Ploygon plot for maps, scatter plot for city points/labels; scatter x=xc y=yc / markerattrs=(symbol=circlefilled color=black size=14pt) datalabel=city datalabelattrs=(weight=bold size=16pt color=black); ... View more

  Over on Twitter, Robert Allison reminds us today's a special day with his chart of 13 Years of Friday the 13th! So, to commemorate the occasion, here's a quick SAS ODS Graphics take on the box office hauls of the films in the Friday the 13th movie franchise, complete with a "slashed" (aka broken) y-axis.   SAS CODE *==> Fun w/SAS ODS Graphics: Friday the 13th Movie Franchise Box Office (excludes 2009 "reboot") Movie poster images from Google search, Box office figures from Wikipedia; %let datasymbols=; * Map movies to their corresponding thumbnail movie poster images; %macro gensymbols; %do p=1 %to 11; symbolimage name=symP%sysfunc(putn(&p,z2)) /* Offset cropped images so tops approximate height */ image="/folders/myfolders/Friday13th/%sysfunc(putn(&p,z2)).jpg" / voffset=-.5; %let datasymbols=&datasymbols symP%sysfunc(putn(&p,z2)); %end; %mend; data movies; * Create SAS dataset from inline data; infile datalines dlm="|" dsd; input title : $50. released : yymmdd10. boxoffice; format released yymmdd10. boxoffice dollar6.2; year=year(released); x+1; * Used to link box movie info to movie poster images; datalines; Friday the 13th|1980-05-09|59.75 Friday the 13th Part 2|1981-04-30|21.72 Friday the 13th Part III|1982-08-13|36.69 Friday the 13th: The Final Chapter|1984-04-13|32.98 Friday the 13th: A New Beginning|1985-03-22|21.93 Friday the 13th Part VI: Jason Lives|1986-08-01|19.47 Friday the 13th Part VII: The New Blood|1988-05-13|19.17 Friday the 13th Part VIII: Jason Takes Manhattan|1989-07-28|14.34 Jason Goes to Hell: The Final Friday|1993-08-13|15.94 Jason X|2002-04-26|16.96 Freddy vs. Jason|2003-08-15|114.9 ; * Create composite vbar/scatter plot (use movie poster images as markers); ods listing gpath='/folders/myfolders'; ods graphics on / reset antialias width=10in height=8in attrpriority=none imagename="FridayThe13th"; proc sgplot data=Movies noborder noautolegend; %gensymbols; styleattrs datasymbols=(&datasymbols); vbarbasic year / barwidth=.5 response=boxoffice nooutline fillattrs=(color=red) fill datalabel; scatter x=year y=boxoffice / group=x markerattrs=(size=75pt); inset "Box Office (Millions)" / title=" Friday the 13th Movie Franchise " textattrs=(size=20pt) titleattrs=(size=26pt) position=top valuealign=center; xaxis display=(nolabel noticks noline) type=discrete offsetmin=.05 offsetmax=.05; * "Slash" the y-axis (broken-axs); yaxis display=(nolabel noticks) offsetmin=.001 offsetmax=.001 ranges=(0-62 100-117) grid tickvalueformat=dollar4.;   ... View more

Earthrise (Bubble/Ellipse Plot) Earthrise (Apollo 11) Figured I should do something to mark the 50th anniversary of the Moon Landing, even if it's only a silly Earthrise bubble/ellipse plot!    * Fun w/SAS ODS Graphics: Apollo 11 Earthrise Bubble/Ellipse Plot See the real Apollo 11 deal at history.nasa.gov/ap11ann/kippsphotos/6550.jpg; ods graphics / height=5in width=5in antialias; * Original NASA image was 600 x 600; data moon; * Generate a point for the earth; x=300; y=185; r=50; * x, y, radius for bubble plot of "earth"; proc template; * "Earth Rise" chart (bubble + eclipse plot); define statgraph ellipseparm; begingraph / opaque=true border=false backgroundcolor=black pad=0 subpixel=on; * Black background; layout overlayequated / equatetype=square opaque=true border=true backgroundcolor=black wallcolor=black WALLDISPLAY=(FILL) commonaxisopts=(viewmin=0 viewmax=600) /* Limit X-Y axes range to 0=600 */ xaxisopts=(display=none offsetmin=0 offsetmax=0) yaxisopts=(display=none offsetmin=0 offsetmax=0 reverse=true); * Got coordinates from MS-Paintbrush, so use reverse y-axis; bubbleplot x=x y=y size=r / fillattrs=(color=dodgerblue) display=(fill) dataskin=sheen relativescale=false; * Dodger Blue Earth; ellipseparm xorigin=300 yorigin=6250 semimajor=6000 semiminor=6000 slope=0 / display=(fill) fillattrs=(color=cx696969); * Dim Gray Moon (semimajor=semiminor=circle); endlayout; endGraph; end; run; proc sgrender data=moon template=ellipseparm; ... View more

  In advance of Independence Day, here's a quick SAS ODS Graphics bar-and-scatter-plots take on Three Flags by Jasper Johns. Enjoy the holiday, all!   CODE * Fun with SAS ODS Graphics: "Three Flags" for July 4th A SAS ODS Graphics take on "Three Flags" by Jasper Johns (whitney.org/collection/works/1060); data stripes; * Generate points for flag stripes; do y1=1 to 14; x=mod(y1,2); * Alternate red and white stripes (bar length of 1 or 0); output; end; data canton; * Generate points for canton (blue rectangle in upper left corner); do y2=1 to 14; x=.4625*(y2<=7); * Blue over first 7 stripes covering 46.25% of flag width; output; end; data stars; * Generate points for 48 stars (arranged in 6 rows of 8 stars - painted in 1958); do r=1 to 6; y3=.79+(r-1)*1.08; * Distribute vertically over first 7 stripes; do c=1 to 8; x=.0355+0.96*.4625/8*(c-1); * Distribute horizontally over width of canton; output; end; end; data flag; * Merge flag objects together; set stripes canton stars; ods listing gpath='/folders/myfolders'; * Save generated flag image; ods graphics on / reset antialias width=7.38in height=5in imagename="oneflag" border=off; * Width is 1.476 x height; proc sgplot data=flag noautolegend pad=0 nowall noborder subpixel; * Generate a flag - two bars charts + scatter plot; styleattrs backcolor=CXF0F0D8; * Background (Beige); symbolchar name=uniStar char='2605'x; * Unicode value for 5-pointed star; hbarbasic y1 / response=x barwidth=1 nooutline fillattrs=(color=CXB22234 transparency=.2) nozerobars discreteoffset=.5; * Stripes (Red); hbarbasic y2 / response=x barwidth=1 nooutline fillattrs=(color=CX3C3B6E) nozerobars discreteoffset=.5; * Canton (Blue); scatter x=x y=y3 / y2axis markerattrs=(symbol=uniStar color=CXF0F0D8 size=41pt); * Stars (Beige); xaxis display=none offsetmin=0 offsetmax=0; yaxis display=none offsetmin=0 offsetmax=0 type=discrete; y2axis display=none type=linear min=0 max=13 offsetmin=0 offsetmax=0; data threeflags; * Stitch trio of flag images together, each about 24 % smaller than previous one; x=.5; y=.5; proc template; * Use GTL to overlay the three flag images on top of a "dummy" scatter plot; define statgraph image; begingraph / subpixel=on border=false; layout overlay; scatterplot x=x y=y / markerattrs=(size=0); * "Dummy" scatter plot; drawimage "/folders/myfolders/oneflag.png" / border=true borderattrs=(color=black thickness=4pt) height=100 heightunit=percent scale=fitheight drawspace=graphpercent; * Back flag is 100% of original image size; drawimage "/folders/myfolders/oneflag.png" / border=true borderattrs=(color=black thickness=.05pt) height=76 heightunit=percent scale=fitheight drawspace=graphpercent; * Middle flag is about 25% of original image size; drawimage "/folders/myfolders/oneflag.png" / border=true borderattrs=(color=black thickness=.05pt) height=58 heightunit=percent scale=fitheight drawspace=graphpercent; * Front flag is about 25% of previous image size; endlayout; endgraph; end; ods graphics on / reset antialias width=7.38in height=5in imagename="threeflags" border=off; * Width is 1.476 x height; proc sgrender data=threeflags template=image;   Original Masterpiece    Bar and Scatter Plot Knockoff ... View more

  Happy Easter, SAS Peeps!   BUILT-FROM-ELLIPSES EASTER BUNNY - BEFORE FILL     SAS CODE   * Fun w/SAS ODS Graphics: Built-From-Ellipses Easter Bunnies; %macro bunny(); ods _all_ close; * Animated GIF setup (pause for 2 secs on 1st frame); options papersize=('5 in', '5 in') printerpath=gif animation=start nodate nonumber animduration=2 animloop=YES NOANIMOVERLAY; ods printer file='/folders/myfolders/peeps.gif' ; * Animated GIF image filename; ods graphics / height=5in width=5in imagefmt=gif antialias; %do f=1 %to 21; * Generate 21 "frames" (Easter Bunnies) for animated GIF; %let title=%substr(%str( HAPPY EASTER, PEEPS!),1,&f); * Title (one letter at a time); %if &f=1 %then %let title=HAPPY EASTER, PEEPS!;; * Full title on first frame; %let letter=%substr(%str( HAPPY EASTER, PEEPS!),&f,1); * New letter in title; /* Lifewire: The Colors of Easter https://www.lifewire.com/the-colors-of-easter-4126619 Light Plum, Sky Blue, Turquoise, Lemon, Spring Green, Mint, Pale Pink, Light Salmon, Violet, Lilac, Robin's Egg Blue, Lavender */ %let color=CX%scan(%str(dda0dd,87ceeb,40e0D0,fff700,00ff7f,bdfcc9,ffc0cb,ffa07a,ee82ee,c8a2c8,00cccc,b573dc,dda0dd,87ceeb,40e0D0,fff700,00ff7f,bdfcc9,ffc0cb,ffa07a,00ff7f),&f,","); data peeps; * Generate point for text plot to show each new letter in title; x=220; y=640; label="&letter"; proc template; * Draw built=from-ellipses Easter Bunny (using points obtained via MS-Paintbrush grid); define statgraph ellipseparm; begingraph / opaque=true border=false backgroundcolor=WHITE; * Smoke white background; entrytitle "&title"; layout overlayequated / opaque=true border=false backgroundcolor=WHITE wallcolor=WHITE WALLDISPLAY=(FILL) xaxisopts=(display=none offsetmin=0.05 offsetmax=0.05) yaxisopts=(display=none offsetmin=0.05 offsetmax=0.05 reverse=true); * y-values lower at top in Paintbrush, so reverse; ellipseparm xorigin=220 yorigin=640 semimajor=220 semiminor=220 slope=0 / fillattrs=(color=&color.) display=(fill); * Body (semimajor=semiminor=circle); ellipseparm xorigin=220 yorigin=397 semimajor=185 semiminor=120 slope=0 / fillattrs=(color=&color.) display=(fill); * Head; ellipseparm xorigin=130 yorigin=148 semimajor=85 semiminor=210 slope=-.4 / fillattrs=(color=&color.) display=(fill); * Left Ear (use -slope to tilt left); ellipseparm xorigin=310 yorigin=148 semimajor=85 semiminor=210 slope=.4 / fillattrs=(color=&color.) display=(fill); * Right Ear (use +slope to tilt right); ellipseparm xorigin=220 yorigin=430 semimajor=15 semiminor=15 slope=0 / fillattrs=(color=black) display=(fill); * Nose; ellipseparm xorigin=155 yorigin=392 semimajor=15 semiminor=15 slope=0 / fillattrs=(color=black) display=(fill); * Left Eye; ellipseparm xorigin=285 yorigin=392 semimajor=15 semiminor=15 slope=0 / fillattrs=(color=black) display=(fill); * Right Eye; textplot x=X y=y text=label / position=center strip=true textattrs=(color=white size=60pt weight=bold); * Display new letter of message; endlayout; endGraph; end; run; proc sgrender data=peeps template=ellipseparm; run; options animduration=.6; %end; options printerpath=gif animation=stop; * Animated GIF wrap-up; ods printer close; %mend; %bunny; * Draw Bunnies!; GETTING ELLIPSE POINT COORDINATES WITH PAINT To get the rough coordinates of the points needed for the ellipses, I just cut-and-pasted a picture of a 'Peeps' cereal box into Microsoft Paintbrush and hovered over the pertinent points to display the coordinates in the status bar.     ... View more

    To mark the close of the old year and beginning of the new, here's a short program for a Pac-Man inspired New Year's Eve countdown timer. Who says you shouldn't use Pie Charts?  Happy New Year, all!   Btw, for a compilation of past Fun-With-SAS-ODS-Graphics posts, check out Oh, There’s No Place Like SAS ODS Graphics for the Holidays! Silly, yes, but also illustrates some techniques you may find useful for creating reproducible and scalable graphics in your Day Job!   * Fun With SAS ODS Graphics: Pac-Man inspired New Years Eve 2019 Pie Chart countdown timer; data slices; * Create SAS data for pie charts; do r=1 to 10; * 10-9-8-7-6-5-4-3-2-1 countdown (seconds); do p=0 to 20 by 2.5, 20 to 0 by -2.5; * 18 charts for each second (mouth closed=0% to wide-open=20%, and back again!); n+1; * Assign unique # to each chart (2 slices/chart); slice=1; pct=p; output; * 1st slice for "mouth"; slice=2; pct=100-p; output; * 2nd slice for remainder; end; end; proc template; * Template for New Years Eve countdown timer chart (pie chart + text); define statgraph pacman; mvar Yr remainSec unicodeSec xLoc Msg; * Parameters: year, remaining seconds, Unicode secs, x-location (%), Happy New Year message; begingraph / subpixel=on border=false datacolors=(black CXFFFF00) backgroundcolor=black pad=0in; * Like game, make "Pac-Man" yellow, with black "mouth"; layout region / border=false backgroundcolor=black outerpad=1in; * Provide some padding for text; piechart category=slice response=pct / datalabelcontent=none centerfirstslice=true display=(fill); drawtext textattrs=(size=36pt color=CX00ff00 weight=bold) Yr / y=100 anchor=top widthunit=percent width=100 yspace=graphpercent xspace=graphpercent; * Green year; drawtext textattrs=(size=16pt color=CX00ff00 weight=bold) remainSec / y=5 anchor=bottom widthunit=percent width=100 yspace=graphpercent xspace=graphpercent; * Green seconds remaining; drawtext textattrs=(size=24pt color=CXFFFF00 weight=bold family="Arial Unicode MS") {unicode unicodeSec} / y=50.5 anchor=right x=xLoc widthunit=percent width=100 yspace=graphpercent xspace=graphpercent; * Yellow "dots" with seconds remaining; drawtext textattrs=(size=24pt color=CX1919A6 weight=bold) Msg / justify=center y=50 anchor=center x=50 widthunit=percent width=50 yspace=graphpercent xspace=graphpercent; * Blue-ish Happy New Year message; endlayout; endgraph; end; run; ods _all_ close; * Animated GIF setup (pause for 1.5 secs on 1st frame); options papersize=('3 in', '3 in') printerpath=gif animation=start nodate nonumber animduration=1.5 animloop=NO NOANIMOVERLAY; ods printer file='/folders/myfolders/PacMan.gif'; * Animated GIF image filename; ods graphics / height=3in width=3in imagefmt=gif antialias; %macro genframes; /* Generate 180-frame coundown (10 seconds * 18 frames/second) */ %do s=10 %to 1 %by -1; /* Generate set of 18 frames for each second in 10-second countdown */ %let UnicodeSec="%scan(2460 2461 2462 2463 2464 2465 2466 2467 2468 2469, &s, ' ')"x; /* Unicode values for animated dots filled with remaining secs (10, 9, ..., 1) */ %let RemainSec="Time Remaining: 00:%sysfunc(putN(&s,z2.))"; /* Display remaining secs at bottom of image */ %do c=1 %to 18; /* Generate frames for one second interval */ %let n=%eval((&s-1)*18+&c); /* Determine pie chart # (used to select data generated earlier) */ %let xLoc=%scan(100 98 96 94 92 90 88 86 84 82 80 78 76 74 72 70 68 66, &c, ' '); /* X-coordinates (%'s) for animated dot that show seconds remaining */ proc sgrender data=slices(where=(n=&n)) template=pacman; /* Create pie charts! */ run; options animduration=.056; /* Try to make total duration of 10 seconds-ish (180 images x .056 sec/image) */ %end; %end; %mend; %let Yr="2018"; * Display old year (2018); %let Msg=""; * No Happy New Year message until end; %genframes; * Issue macro call to generate 180 frames!; options animduration=1.5; * Display "GAME OVER" frame at bottome of screen for 1.5 seconds; %let RemainSec="GAME OVER"; proc sgrender data=slices(where=(n=180)) template=pacman; * Create pie chart|; run; options animduration=3; * For 3 seconds, show final frame...; %let Yr="2019"; * Display new year (2019) at top of screen; %let Msg="HAPPY NEW YEAR!"; * Display "HAPPY NEW YEAR" over pie chart in middle of screen; %let RemainSec="PRESS START"; * Display "PRESS START" at bottom of screen; proc sgrender data=slices(where=(n=180)) template=pacman; * Create pie chart!; run; options printerpath=gif animation=stop; * Animated GIF wrap-up; ods printer close;   ... View more

  Not sure whether to be proud or ashamed of this  but here are Microsoft OneDrive and Google Drive links to view the PowerPoint for my SCSUG 2018 presentation, "Oh, There’s No Place Like #SAS ODS Graphics for the Holidays!" (caution: 33+ MB, so download with care!). The PowerPoint includes static images + animated GIFs, as well as links to the SAS code that produced them. The 8MB accompanying paper includes Twitter screenshots of images, as well as links to the SAS code. Btw, Microsoft's PowerPoint Online displays the slides better than Google's viewer does, but you'll have to download the big .pptx file and open it in "real" PowerPoint client software to get the animated GIFs to animate (static versions are displayed online). Happy Holidays!      ... View more

ANIMATED GIF OUTPUT (Click on image or magnifying glass to view larger image at full speed!)   How Santa solves his massive traveling salesman problem remains a mystery (SAS Viya on quantum computers?). So, in the spirit of Santa's DasherBoard (SAS VA-based, from @TravisMurphy), Optimal Santa Route (SAS VA-based, from, @FalkoSchulz), and Santa's Information Dashboard (SAS/GRAPH-based, from @ RobertAllison_S AS), here's my decidedly less-ambitious SAS ODS Graphics-based 'Traveling Santa Tour' through the capital cities of the contiguous 48 states using an optimized 'flight plan' lifted from a PROC OPTNET example. Happy Holidays, all!   CODE * Fun With SAS ODS Graphics: Animated "Traveling Santa Tour" Thru U.S. Capital Cities; data FlightPlan; * Create SAS dataset of cities in Santas "flight plan"; infile datalines dlm=',' dsd; * Cities are output of PROC OPTNET example in SAS documentation; input capital : $30.@@; * See "Example 2.7 Traveling Salesman Tour through US Capital Cities"; if capital^=''; * support.sas.com/documentation/cdl/en/ornoaug/66084/HTML/default/viewer.htm#ornoaug_optnet_examples07.htm; seq+1; * Add sequence number for tour path; datalines; Montgomery, Montgomery, Tallahassee, Columbia, Raleigh, Richmond, Washington, Annapolis, Dover, Trenton, Hartford, Providence, Boston, Concord, Augusta, Montpelier, Albany, Harrisburg, Charleston, Columbus, Lansing, Madison, Saint Paul, Bismarck, Pierre, Cheyenne, Denver, Salt Lake City, Helena, Boise City, Olympia, Salem, Sacramento, Carson, Phoenix, Santa Fe, Oklahoma City, Austin, Baton Rouge, Jackson, Little Rock, Jefferson City, Topeka, Lincoln, Des Moines, Springfield, Indianapolis, Frankfort, Nashville, Atlanta, Montgomery ; proc sql; * Get Geo info for cities from SAS USCITY map dataset (SASFILES points to SAS Map datasets); create table capitalsGeo(keep=seq capitalCity stateCode x y lat long rename=(x=xC y=yC)) as select f.seq, upcase(case when f.capital='Carson' then 'Carson City' else f.capital end) as capitalCity, /* Fix: add 'City' to 'Carson' */ u.StateCode, u.x, u.y, u.lat, u.long from flightplan f join sasfiles.uscity u on u.capital='Y' and f.capital=u.city order by seq; * Create dataset with info needed to chart Santas route; data CapitalsGeoLines(keep=seq xC yC x1 y1 x2 y2 CapitalCity: StateCode: miles); set CapitalsGeo end=eof; * Get data for "FROM" city from prior observation; x1=lag(xC); y1=lag(yC); lat1=lag(lat); long1=lag(long); CapitalCity1=lag(CapitalCity); StateCode1=lag(StateCode); x2=xC; y2=yC; * X2/Y2=points for "TO" City; if ^eof then CapitalCityLbl=CapitalCity;* Do not repeat starting city name at end of tour; miles=geodist(lat1, long1, lat, long, 'M'); if _n_>1; * No line for first point; seq=seq-1; * Re-adjust seq # to account for skipped point; data usa(keep=state segment x y pid); * Get outlines of 48 contiguous U.S. states from SAS US dataset; set sasfiles.us(where=(statecode not in('AK','HI', 'PR'))); * Remove AK HI PR; length pid $ 8; * State outlines are one or more polygons; pid=put(state, 3.0)||put(segment, 3.0); * Generate unique polygon identifier; * For more info, see 'Micro Maps' at blogs.sas.com/content/graphicallyspeaking/2015/04/19/micro-maps/; data capitalGeoUSA; * Merge city points and state outline data; set usa capitalsGeoLines; proc template; * Template for Santa Tour Map (chart is composite of multiple plots); define statgraph Maps; mvar FromCity ToCity Distance TotDistance; * Macro variables for inset box; begingraph / subpixel=on border=false; * Microsoft Windows 10 Santa emoji image (no art budget!); symbolimage name=santa image='/folders/myfolders/Windows10Santa.png'; layout overlayEquated / walldisplay=none border=false xaxisopts=(offsetmin=0.0 offsetmax=0.0 display=none) yaxisopts=(offsetmin=0.0 offsetmax=0.0 display=none); * 1. Polygon plot of 48 contiguous U.S. states; polygonPlot x=x y=y id=pid / display=(fill outline) fillattrs=(color=CXc4cace) outlineattrs=(color=white); * "Metallic silver" fill; * 2. Vector plot of lines between cities visited; vectorplot x=x2 y=y2 xorigin=x1 yorigin=y1 / lineattrs=(color=black); * 3. Scatter plot of cities visited - circles; scatterplot x=xC y=yC / markerattrs=(symbol=circle color=black) datalabel=capitalCityLbl datalabelattrs=(color=black); * 4. Scatter plot of last city visited - Win10 Santa emoji marker; scatterplot x=xS y=yS / markerattrs=(symbol=santa color=black size=30); * 5. Chart title; entry halign=center " A Traveling Santa Tour of U.S. State Capitals " / valign=top textattrs=(size=15 color=CX296e01 weight=bold); * "Metallic green"; * 6. Inset box of from/to cities and distances in lower left corner; layout gridded / rows=4 order=columnmajor valign=bottom halign=left border=true; entry halign=left "From: " / textattrs=(size=11 color=CXa62c2b weight=bold); * "Metallic red"; entry halign=left "To: " / textattrs=(size=11 color=CXa62c2b weight=bold); entry halign=left "Distance: " / textattrs=(size=11 color=CXa62c2b weight=bold); entry halign=left "Total: " / textattrs=(size=11 color=CXa62c2b weight=bold); entry halign=left FromCity / textattrs=(size=11 color=CXa62c2b weight=bold); entry halign=left ToCity / textattrs=(size=11 color=CXa62c2b weight=bold); entry halign=left Distance / textattrs=(size=11 color=CXa62c2b weight=bold); entry halign=left TotDistance / textattrs=(size=11 color=CXa62c2b weight=bold); endlayout; * 6. Larger Win10 Santa emoji image in lower right corner; drawimage "/folders/myfolders/Windows10Santa.png" / anchor=bottomright x=100 y=0 drawspace=wallpercent height=18 heightunit=percent scale=fitheight; endlayout; endgraph; end; run; ods _all_ close; * Create animated GIF for 'Santa Tour'; options papersize=('11 in', '6.6 in') printerpath=gif animation=start nodate nonumber animduration=2 animloop=NO NOANIMOVERLAY; ods printer file='/folders/myfolders/Santa.gif'; ods graphics / height=6.8in width=11in imagefmt=gif antialias labelmax=2000; * Specify labelmax to minimize city name "collisons"; %macro drawlines; proc sql noprint; * Get number of frames to draw (# of from/to pairs); select count(*) into :numlines from CapitalsGeoLines; quit; %do l=1 %to &numlines; * Draw "frames" containing image with cities visited thus far; %if &l=1 | &l=&numlines %then options animduration=2; %else options animduration=.75;; * Make first/last frames longer; data capitalGeoUSAwork; * Dataset for frame 1 has 1 city, 2 has 2 cities, ...; set capitalGeoUSA(where=(seq<=&l)); totmiles+miles; * Tally cumulative miles; if seq=&l then do; * Create points, macro variables for last line drawn; xS=xC; yS=yC; * X/Y coordinates for Santa Claus location; call symput('FromCity',trim(CapitalCity1)||', '||StateCode1); * From city/state; call symput('ToCity',trim(CapitalCity)||', '||StateCode); * To city/state; call symput('Distance',compress(put(miles,comma6.),' ')||' mi'); * Miles (current leg); call symput('TotDistance',compress(put(totmiles,comma6.),' ')||' mi'); * Cumulative miles; end; * Render frame; proc sgrender data=capitalGeoUSAwork template=Maps; run; %end; %mend; %drawlines; * Draw the frames; options printerpath=gif animation=stop; * All done - close animated GIF file; ods printer close;   MICROSOFT WINDOWS 10 SANTA EMOJI (Windows10Santa.png)   ... View more

COMPLETE IMAGE   DETAIL After seeing Sanjay Matange's blog post on the use of a text plot to create a clock face, thought I'd take a shot at creating time-and-temperature clock faces from temperature forecast XML data obtained from weather.gov (cold week ahead in Chicago!). OK, it's no threat to the cool Apple Watch Weather app, but it does have a literal "heat map" going for it.    Btw, presenting temperature data "radially" was the subject of a recent Microsoft Research study. And way back in 2012, Sanjay Matange and Robert Alison respectively presented their SAS ODS Graphics and SAS/GRAPH takes on Ken Kleinman's interesting R-based approach.   CODE *==> Fun w/SAS ODS Graphics: Time and Temperature "Watch Faces"; *==> Grab 7 days of forecasted hourly temperatures from weather.gov XML (Chicago); filename TimeTemp url 'https://forecast.weather.gov/MapClick.php?lat=41.8843&lon=-87.6325&FcstType=digitalDWML'; filename map '/folders/myfolders/TimeTemp.map'; * XML map to grab times/temps; libname TimeTemp xmlv2 xmlmap=map; data timetemp; * Create dataset with 5 days of hourly temperatures; set timetemp.times; * Merge time and temperature data; set timetemp.temps(where=(type='hourly')); * Just hourly temps; hr=hour(time); * Hour (0-23); if hour(time)< 12 then AmPm='AM'; else AmPm='PM'; * AM or PM?; hr12=mod(hr,12); * Change 0-23 to 0-11; day=put(datepart(time),downame3.); * Grab day of week; if '12nov2018'd<=datepart(time)<='16nov2018'd; * Select 5 complete days; *==> Calc x/y points for "watch faces"; data watch; array hrs[12] _TEMPORARY_ (3 2 1 0 11 10 9 8 7 6 5 4); * Map 30-degree wedges to hour of day; do w=0 to 11; * Create 12 polygons for wedges; hr12=hrs[w+1]; * Find hour for 30 degree wedge; seq=0; * Polygon ID; x=0; y=0; output; * Polygons start at 0,0; do a=2*constant("pi")/12*(w-.5) to 2*constant("pi")/12*(w+.5) by 0.01; seq+1; * Polygon sequence number; x=cos(a); * Assume circle is radius=1; y=sin(a); output; end; xd=x; yd=y; * 0,0 to xd,yd are "divider line" points for wedges; x=.; y=.; output; xd=.; yd=.; a=2*constant("pi")/12*w; * xt, yt are points for temperature values; xt=1.125*cos(a); * Place temperatures outside of unit circle (radius=1); yt=1.125*sin(a); output; xt=.; yt=.; end; proc sql; * Merge temperature and "watch face" data; create table tempswatch as select w.*, put(t.temp,3.) as txt, t.day, t.temp, t.time, t.ampm from timetemp t join watch w on t.hr12=w.hr12 order by datepart(t.time), w, seq; ods _all_ close; ods listing gpath='/folders/myfolders' image_dpi=300; ods graphics / reset imagename='TimeTemp' height=15in width=6in noborder antialias imagefmt=png; * Create chart - AM/PM panels for each day; proc sgpanel data=tempswatch noautolegend aspect=1 subpixel; panelby ampm day / onepanel uniscale=all columns=2 layout=lattice noheader noborder sort=data; * Display day of week, AM/PM indicator in upper left corner; inset Day ampm / nolabel position=topleft textattrs=(color=black weight=bold); * Polygon plot hourly wedges - fill colors vary with temperature; polygon x=x y=y id=w / fill nooutline colorresponse=temp colormodel=(lightblue darkblue); * Vector plot divider lines; vector x=xd y=yd / lineattrs=(color=white thickness=1.5 pattern=solid) noarrowheads; * Text plot temperature values; text x=xt y=yt text=txt / position=center textattrs=(color=black weight=bold) strip; colaxis display=none min=1.2 max=1.2; rowaxis display=none min=1.2 max=1.2;   XML MAP <?xml version="1.0" encoding="UTF-8"?> <SXLEMAP version="1.2" name="SXLEMap"> <!-- Create TIME/TEMP SAS Datasets From Weather.gov XML --> <TABLE name="TIMES"> <TABLE-PATH>//data/time-layout/start-valid-time</TABLE-PATH> <COLUMN name="TIME"> <PATH>//data/time-layout/start-valid-time</PATH> <TYPE>numeric</TYPE> <DATATYPE>date</DATATYPE> <FORMAT width="19">E8601DT</FORMAT> <INFORMAT width="19">B8601DT</INFORMAT> </COLUMN> </TABLE> <TABLE name="TEMPS"> <TABLE-PATH>//data/parameters/temperature/value</TABLE-PATH> <COLUMN name="TYPE" retain="YES"> <PATH>//data/parameters/temperature/@type</PATH> <TYPE>character</TYPE> <DATATYPE>string</DATATYPE> <LENGTH>20</LENGTH> </COLUMN> <COLUMN name="TEMP"> <PATH>//data/parameters/temperature/value</PATH> <TYPE>numeric</TYPE> <DATATYPE>integer</DATATYPE> <LENGTH>8</LENGTH> </COLUMN> </TABLE> </SXLEMAP> ... View more

      Fretting over what to wear when you go back to work on Tuesday? The Emily Post Institute has advice on Wearing White After Labor Day. Happy Labor Day!   (Update: Added PERCENT keyword to PIECHART statement)   * Fun With SAS ODS Graphics: Can You Wear White After Labor Day? Reference: Wearing White After Labor Day, emilypost.com/advice/wearing-white-after-labor-day; proc format; * Special formats; picture mmmdd other='%b%0d' (datatype=date); * Dates=mmmdd; picture ddd other='000 Days'; * Days=nnn Days; data When2WearWhite(keep=DateRange Days); * Create data to chart; LaborDay=holiday('LABOR',year(today())); * When is Labor Day this year?; MemorialDay=holiday('MEMORIAL',year(today())); * When is Memorial Day this year?; length DateRange $ 30.; * Date range label for chart; DateRange=put(MemorialDay,mmmdd.)||'-'||put(LaborDay,mmmdd.); DaysInYear=intck('days', /* How many days in this year? */ intnx('years',today(),-1,'e'),intnx('years',today(),0,'e')); call symput('degrees', /* Compute degrees relative to piechart origin for start of Memorial-Labor Day pie slice */ put(360-(substr(put(MemorialDay,julian5.),3,3)/DaysInYear*360-90),3.)); Days=LaborDay-MemorialDay+1; * Days between Memorial Day & Labor Day (Definitely OK to wear white!); output; Days=DaysInYear-Days; * Rest of year (Maybe OK to wear white?); DateRange='Jan01-'||put(MemorialDay-1,mmmdd.)||','||put(LaborDay+1,mmmdd.)||'-Dec31 ?'; output; ods graphics on / reset width=8in height=5in; * Create a pie chart with two slices (Memorial-Labor Day, rest of year); proc template; define statgraph laborday; begingraph / datacolors=(white lightgray); * Memorial Day to Labor Day=White (OK), rest of year=LighGray (OK??); entrytitle 'CAN YOU WEAR WHITE AFTER LABOR DAY?'; layout region; /* Pie wedges proportional to number of days in date range */ piechart category=DateRange response=Days / /* Position Memorial-Labor Day wedge as if year is a clock w/Jan. 1 at 12:00 */ start=&degrees categorydirection=clockwise dataskin=crisp datalabelcontent=(category response percent) datalabelattrs=(color=black) datalabellocation=callout; endlayout; endgraph; end; run; proc sgrender data=When2WearWhite template=laborday; format Days ddd.; run;   ... View more