Fold line of length 2*PI*r into a circle Inspired by a @Rick_SAS  blog post, here's a last-minute PI Day chart (hey, it's still PI Day in Hawaii!  ).   *==> Fun w/SAS ODS Graphics: Fold line of length 2*PI*R into a circle; data circle; * Calc x/y points for circle of N points; retain xR yR xRO 0 yRO 1 tR "r=1"; * Used to draw unit circle radius; seq=0; xo=-constant("pi"); yo=-1; x=-xo; y=yo; output; x=.; * Starting line segment of length 2*PI; N=40; * Number of line segments in circle (40+ is good, multiples of 4); d=constant("pi"); * Loop for points on right side of circle; do a=-2*constant("pi")/4 to 2*constant("pi")/4 by 2*constant("pi")/N; x=cos(a); * Assume unit circle with radius=1; y=sin(a); xo=cos(lag(a)); * Prior x/y point; yo=sin(lag(a)); if xo^=. then do; * Output circle line segments; seq+1; d=d-2*constant("pi")/N; * "Unused" line segment (drawn perpindicular to last points on circle); m=-1/(y/x); * Slope of tangent to current point on circle; if m=. then do; * Slope=infinity (x=0)?; xT=x; * Yes, calculate x/y endpoint of "unused" line; yT=y+d; end; else do; * No, calculate x/y endpoint of "unused" line; xT=x+sign(m)*d/sqrt(1+m**2); * Slope determines x direction; yT=y+abs(m)*d/sqrt(1+m**2); * y direction always positive; end; output; x=-x; * Negate x values for left side of circle; xo=-xo; xT=-xT; output; end; end; proc sql; * Determine frames needed for GIF; create table framenums as select distinct seq as frame from circle; create table frames as /* Create frames with needed points */ select frame, seq, x, y, xo, yo, xR, yR, tR, xRO, yRO, case when f.frame=c.seq then xT end as xT, case when f.frame=c.seq then yT end as yT from framenums f, circle c where f.frame>=c.seq and (f.frame=0 or c.seq>0) order by 1, 2; ods _all_ close; * Animated GIF setup; options papersize=('5 in', '5 in') printerpath=gif animation=start nodate nonumber animloop=YES NOANIMOVERLAY animduration=.5; ods printer file='/folders/myfolders/piDayV2.gif'; ods graphics / width=5in height=5in imagefmt=GIF; %macro genCircle(start=, stop=, dur=); options nobyline animduration=&dur; * Create plots; proc sgplot data=frames aspect=1 noautolegend; by frame; where &start<=frame<=&stop; vector x=xR y=yR / xorigin=xRO yorigin=yRO /* Radius */ noarrowheads datalabel=tR datalabelpos=bottom lineattrs=(pattern=shortdash); vector x=x y=y / xorigin=xo yorigin=yo noarrowheads; /* Segments of circles */ vector x=x y=y / xorigin=xT yorigin=yT noarrowheads; /* Remaining "unused" line segment */ xaxis display=(nolabel) min=-3.25 max=3.25 offsetmin=.01 offsetmax=.01 values=(-3.14 -2 -1 0 1 2 3.14) valuesdisplay=("-(*ESC*){Unicode pi}" "-2" "-1" "0" "1" "2" "(*ESC*){Unicode pi}"); yaxis display=(nolabel) min=-3.25 max=3.25 offsetmin=.01 offsetmax=.01 values=(-3.14 -2 -1 0 1 2 3.14) valuesdisplay=(" " "-2" "-1" "0" "1" "2" " "); inset "2(*ESC*){Unicode pi}r" / position=top textattrs=(size=18pt); inset "Circumference" / position=bottom textattrs=(size=18pt); run; %mend; %genCircle(start=0, stop=0, dur=.5); * Start frame (.5 seconds); %genCircle(start=0, stop=20, dur=.15); * Draw circle (.15 seconds/frame); %genCircle(start=20,stop=20, dur=.5); * Final frame (.5 seconds); options printerpath=gif animation=stop; * Animated GIF wrapup; ods printer close;   ... View more

Stepping Feet Illusion With the stepping feet illusion racking up likes on Twitter this week, thought I'd take a crack at recreating it using SAS ODS Graphics!    * Fun w/SAS ODS Graphics: Stepping feet illusion (see en.wikipedia.org/wiki/Stepping_feet_illusion); data barsNmarkers; * Generate points for black lines and yellow & blue markers; do frame=1 to 20 by .1; * Move markers between x=2 and x=18 in increments of .1; y=1; do x=1 to 20; * Draw 20 lines 1 unit high; output; end; x=frame; y=.; * Points for yellow and dark blue markers; yy=.75; yb=.25; output; yy=.; yb=.; end; ods _all_ close; * Animated GIF setup; options papersize=('5 in', '2 in') printerpath=gif animation=start nodate nonumber animloop=YES NOANIMOVERLAY animduration=.075; ods printer file='/folders/myfolders/SteppingFeet.gif'; ods graphics / width=5in height=2in imagefmt=GIF; options nobyline; * Stepping feet illusion for first 13 bars (black bars create contrast); proc sgplot data=barsNmarkers(where=(frame<14)) noautolegend noborder nowall; by frame; symbolchar name=uniRect char='25AC'x; * Unicode for rectangle; dropline x=x y=y / lineattrs=(color=black thickness=11pt); scatter x=x y=yy / markeroutlineattrs=(color=yellow) markerattrs=(color=yellow size=144pt symbol=uniRect); scatter x=x y=yb / markeroutlineattrs=(color=darkblue) markerattrs=(color=darkblue size=144pt symbol=uniRect); xaxis display=none offsetmin=0 offsetmax=0 min=0.5 max=21 type=linear; yaxis display=none offsetmin=0 offsetmax=0 min=0 max=1; * No Stepping feet illusion for last 7 bars (white bars create no contrast); proc sgplot data=barsNmarkers(where=(frame>=14)) noautolegend noborder nowall; by frame; symbolchar name=uniRect char='25AC'x; * Unicode for rectangle; dropline x=x y=y / lineattrs=(color=white thickness=11pt); scatter x=x y=yy / markeroutlineattrs=(color=yellow) markerattrs=(color=yellow size=144pt symbol=uniRect); scatter x=x y=yb / markeroutlineattrs=(color=darkblue) markerattrs=(color=darkblue size=144pt symbol=uniRect); xaxis display=none offsetmin=0 offsetmax=0 min=0.5 max=21 type=linear; yaxis display=none offsetmin=0 offsetmax=0 min=0 max=1; run; options printerpath=gif animation=stop; * Animated GIF wrapup; ods printer close; ... View more

Happy New Year! (note: GIF compressed with compressor.io) To mark the new decade, here's a quick-and-dirty SAS ODS Graphics take on a film leader countdown. Happy New Year, all!   * Fun w/SAS ODS Graphics: "Happy New Year!" Film Countdown Leader (see en.wikipedia.org/wiki/Film_leader); data points(keep=counter angle text x y tx ty vx vy); * Counter=10, 9, ..., 1; pi=constant('pi'); * Angle=degrees of slices (0 to 360 by 15); retain tx ty 0; * x/y=Polygon coordinates, tx/ty=text, vx/vy=vector; do counter=10 to 1 by -1; * Countdown seconds, 10 to 1; text=put(counter,2.); * Format econds remaining; do angle=0 to 360 by 15; * Plot every 15 degrees beginning at top; StartRadians=90/360*2*pi-angle/360*2*pi; * Convert radians to x/y points for plotting; vx=cos(StartRadians)*2; vy=sin(startradians)*2; * Last polygon point x/y coordinates; x=0; y=0; output; * Start drawing shaded polygon at center (0,0); do i=0 to angle by 1; * Add more points to shaded plygon (one point/degree); StartRadians=90/360*2*pi-i/360*2*pi; * Convert radians to x/y points for plotting; x=cos(StartRadians)*2; y=sin(startradians)*2; * Polygon point x/y coordinates; output; end; end; end; run; options nobyline; proc template; * Film countdown leader chart definition; define statgraph filmtemplate; * Composite of scatter, polygon, vector, ellipse, text plots; begingraph / pad=0; layout overlayequated / equatetype=square pad=0 border=false walldisplay=none backgroundcolor=lightgrey opaque=true commonaxisopts=(tickvaluelist=(-1 1) viewmin=-1 viewmax=1) /* Cutoff window at -1 and 1 to trim" polygon area */ xaxisopts=(display=none offsetmin=0.000 offsetmax=0.000 viewmin=-1 viewmax=1) /* Suppress axes/ticks/values */ yaxisopts=(display=none offsetmin=0.000 offsetmax=0.000 viewmin=-1 viewmax=1); scatterplot x=x y=y / markerattrs=(size=0); * Ensure at least some data is plotted; polygonplot x=x y=y id=angle / fillattrs=(color=gray) display=(fill) label=none; * Shaded "radar sweep" area; vectorplot x=vx y=vy xorigin=0 yorigin=0 / arrowheads=false lineattrs=(color=black thickness=4pt); * "Radar sweep" hand; referenceline x=0 / lineattrs=(color=black thickness=4pt); * Line bisecting x axis; referenceline y=0 / lineattrs=(color=black thickness=4pt); * Line bisecting y axis; ellipseparm semimajor=.92 semiminor=.92 xorigin=0 yorigin=0 slope=0 / outlineattrs=(color=white thickness=8pt); * Outer circle; ellipseparm semimajor=.77 semiminor=.77 xorigin=0 yorigin=0 slope=0 / outlineattrs=(color=white thickness=8pt); * Inner circle; textplot x=tx y=ty text=text / textattrs=(size=216pt color=black weight=bold) strip=true contributeoffsets=none; * Seconds remaining; endlayout; endgraph; end; run; * Create animated GIF film countdown leader followed by Happy New Year! message; ods _all_ close; * Animated GIF setup; options papersize=('5 in', '5 in') printerpath=gif animation=start nodate nonumber animloop=YES NOANIMOVERLAY; ods printer file='/folders/myfolders/filmlead.gif'; ods graphics / width=5in height=5in imagefmt=GIF; options animduration=.04; * 25 frames per second; proc sgrender data=points template=filmtemplate; * Generate countdown frames using BY statement; by descending counter angle; run; data year; * Text for Happy New Year message; retain tx ty 0 msg "2020" tx2 0 ty2 .9 msg2 "HAPPY NEW YEAR!"; options animduration=2; * Generate Happy New Year! message frame (2 seconds); proc sgplot data=year aspect=1 noborder opaque nowall noautolegend; styleattrs backcolor=black; text x=tx y=ty text=msg / textattrs=(size=128pt color=white weight=bold) contributeoffsets=none; text x=tx2 y=ty2 text=msg2 / textattrs=(size=30pt color=white weight=bold) contributeoffsets=none position=bottom; xaxis display=none values=(-1 0 1); yaxis display=none values=(-1 0 1); run; options printerpath=gif animation=stop; * Animated GIF wrapup; ods printer close; ... View more

SAS ODS Graphics Knockoff DIY Xmas Tree Card Here's a quick X-ACTO knife-free SAS ODS Graphics take on a neat Xmas craft project (Steph Curry ornaments sold separately). Happy Holidays!   * Fun w/SAS ODS Graphics: "3-Pointer" Xmas Tree Inspired by craft project at https://howaboutorange.blogspot.com/2012/12/diy-geometric-christmas-tree-card.html; data tree; w=.9; * Length of sides of equilateral triangles (leave 10% white space); h=.9*sqrt(3)/2; * Height of equilateral triangles (length of sides * sqrt(3)/2); do r=0 to -5 by -1; * 6 rows; do c=r/2 to r/2+abs(r); * 1-6 triangles per row; up=1; * 3 polygon points for triangles pointing up (dark green); id+1; * Polygon identifier; x=c; y=r+.5*h; output; * Point at triangle top; x=c-w/2; y=r-h/2; output; * Point at base (left); x=c+w/2; y=r-h/2; output; * Point at base (right); up=0; * 3 polygon points for triangles pointing down (light green); id+1; * Polygon identifier; x=c; y=r-1.5*h; output; * Point at trianngle top; x=c-w/2; y=r-h/2; output; * Point at base (left); x=c+w/2; y=r-h/2; output; * Point at base (right); end; end; * Plot the "tree" of triangles (polygon plot); ods graphics / reset antialias width=5in height=5in; proc sgplot data=tree noautolegend aspect=1 subpixel noborder; polygon x=x y=y id=id / colormodel=(cx00FF3E cx00B32C) colorresponse=up fill; xaxis display=none offsetmin=.05 offsetmax=.05; yaxis display=none offsetmin=.05 offsetmax=.05; ... View more

  With the help of code snipped from Stephen Mistler's 2014 Producing Normal Density Plots With Shading, here's a quick SAS ODS Graphics takeoff on Andrew Heiss's 2018 R/ggplot Paranormal Distribution. Happy Halloween, all! Btw, it's probably worth mentioning that if you need a density plot - but don't need to make it look like a crudely-drawn Halloween ghost costume - SAS ODS Graphics' DENSITY statement makes things a breeze!   SAS CODE * Fun w/SAS ODS Graphics: "Para-Normal" Density Plot See Stephen Mistler's Producing Normal Density Plots With Shading at blogs.sas.com/content/sastraining/2014/06/10/producing-normal-density-plots-with-shading and Andrew Heiss's R/ggplot Paranormal Distribution at https://twitter.com/andrewheiss/status/1057682358135021568; data pdf; * Points for density plot; retain lower 0; * Used to create area plot; do X=-5 to 5 by .001; * Generate points to plot; PDF=pdf("Normal", x, 0, 1.5); * SAS probability distribution function (normal); output; end; data EyesMouth; * Points for eyes/mouth of "ghost"; circle=unicode('\u25CF'); * Unicode filled circle for eyes/mouth; input x2 y2 size; * Coordinates/sizes for eyes/mouth (eyes smaller than mouth); cards; -.55 .2 1 Left Eye .55 .2 1 Right Eye 0 .15 3 Mouth ; data Ghost; * Merge data together to create our "ghost"; set pdf EyesMouth; * Plot the "ghost"; ods graphics / ANTIALIASMAX=10100; * Band plot (density curve) + text plot (eyes/mouth); proc sgplot data=ghost noautolegend nowall noborder; title height=18pt "PARANORMAL DENSITY PLOT"; styleattrs backcolor=orange; * Pumpkin-ish color; band x=x lower=lower upper=pdf / fillattrs=(color=white); * Ghost body is band plot (shaded normal density plot); text x=x2 y=y2 text=circle / textattrs=(family="Arial Unicode MS") /* Plot ghost eyes/mouth (circles) */ sizemin=32pt sizemax=96pt sizeresponse=size; xaxis display=(nolabel noline noticks) grid gridattrs=(color=white) valueattrs=(color=white); * Use white grid lines, suppress axis line; yaxis display=(nolabel noline noticks) grid gridattrs=(color=white) valueattrs=(color=white) min=0 offsetmin=.01 offsetmax=.01; inset "HAPPY HALLOWEEN!" / position=bottom textattrs=(size=14pt weight=bold); * Halloween greeting at bottom; run;     ... View more

  In addition to hosting SAS conference papers galore, lexjansen.com displays the sizes of the files that contain them. So, here's a quick box plot of the sizes of the MWSUG 2019 papers by section, including my own Oh, There's No Place Like SAS ODS Graphics for the Holidays!, which is revealed to be quite an outlier at 8,887 Kb (it's chock-full of Twitter images). The accompanying (25.6 Mb) PowerPoint e-Poster - with 27 holiday-themed SAS ODS Graphics charts and the code that created them - is available for viewing or downloading on Microsoft OneDrive. The e-Poster's designed for a 41" x 23" screen, so it'll show better on a bigger, higher-resolution monitor or TV!   SAS CODE * Fun with SAS ODS Grapics: MWSUG 2019 Conference Proceedings - Paper Sizes by Section; data mwsugPapers(keep=size section paper); * Scrape paper sections, names, sizes from lexjansen.com; filename mwsug url "https://lexjansen.com/cgi-bin/xsl_transform.php?x=mwsug2019"; infile mwsug recfm=v truncover column=c; input @; _infile_=tranwrd(_infile_,"c2a0"x," "); * Translate non-breaking spaces in web page to Ascii space; length size 3. section paper $ 255.; retain section paper; s=index(_infile_,'<span class="stream">'); * Grab paper section; if s then section=scan(substr(_infile_,s),2,'<>'); input @"<b>Size</b>:" size; * Grab paper size; if size then output; if index(_infile_,'class="paper') then /* Grab paper name */ paper=scan(_infile_,-3,'<>'); proc format; * Display sizes as megabytes; picture mb low-high='00mB' (multiplier=.001); proc summary data=mwsugPapers nway; * Calc # papers, mean/median paper size by section for table; class section; var size; output out=mwsugStats n=N mean=Mean median=Median; proc sql; * Merge paper detail and summary stats; create table mwsugPapers2 as select t1.*, t2.n, t2.mean, t2.median from mwsugPapers t1 left join mwsugStats t2 on t1.section=t2.section; ods graphics on / antialias height=4in width=8in; proc sgplot data=mwsugPapers2 noautolegend; * Box plots of paper sizes by section, show paper name for far outliers; styleattrs datacolors=(CXAEC7E8 CXFFBB78 CX98DF8A CXFF9896 CXC5B0D5 CXC49C94 CXF7B6D2); * Light colors for boxes; hbox size / category=section datalabel=paper labelfar group=section lineattrs=(color=black) meanattrs=(color=black) medianattrs=(color=black) whiskerattrs=(color=black) outlierattrs=(color=black size=9pt) datalabelattrs=(color=black size=9pt); inset "MWSUG 2019 Paper Sizes" / position=topright textattrs=(size=14pt); xaxis display=(nolabel) offsetmin=0 offsetmax=.02 valuesformat=mb. grid values=(0 to 9500 by 1000) valueshint; yaxis display=(nolabel); yaxistable N Mean Median / stat=mean valueattrs=(color=black size=8pt); * Stats pre-calculated and on every obs, so specify stat=mean; format n mean median comma7.; e-Poster Screenshot ... View more

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