SAS BANDPLOT
data GTL_GS_bandGroup; format Date Date9.; do i=0 to 334 by 1; date='01jan2009'd+i; if mod (i, 30) =0 then freq=1; else freq=0; Drug='Drug A'; Val = 16+ 3*sin(i/90+0.5) + 1*sin(3*i/90+0.7); upper=val*1.1; lower=val*0.95; output; Drug='Drug B'; Val = 11+ 3*sin(i/90+0.5) + 1*cos(3*i/90+0.7); upper=val*1.1; lower=val*0.95; output; Drug='Drug C'; Val = 12+ 3*cos(i/90+0.5) + 1*sin(3*i/90+0.7); upper=val*1.1; lower=val*0.95; output; end; run; ODS HTML ; proc template; define statgraph _bandplot; begingraph / /*DATACOLORS=(CXFF0000 green blue) DATACONTRASTCOLORS=(red)*/ subpixel = on; layout overlay / xaxisopts = ( display = (ticks tickvalues)) yaxisopts = ( label = "Response"); bandplot x = date limitupper = upper limitlower = lower/ group = drug display = (fill outline) outlineattrs = ( thickness = 2 pattern = solid ) fillattrs = (transparency = 0.5) name = "a" extend = true; discretelegend 'a' / location = inside across = 1 halign = right valign = top; endlayout; endgraph; end; run; proc sgrender data = GTL_GS_bandGroup template = _bandplot; run; ods _ALL_ close;
本文来自博客园,作者:Iving,转载请注明原文链接:https://www.cnblogs.com/SAS-T/p/15367211.html