About Ksharp

Ksharp · ‎08-14-2025

Try to use FINDW() instead of making a macro variable. data test; input id var1-var3; format var3 f4.2; cards; 1 1 2 3 ; %macro MakeSummary(lib=WORK, dset=, varlist=, outfmt=SummaryFmt); proc sql noprint; create table _varfmt as select name, format as fmtname from dictionary.columns where upcase(libname) eq "%upcase(&lib)" and upcase(memname) eq "%upcase(&dset)" and findw("%upcase(&varlist.)", upcase(strip(name))) ; quit; %mend; %MakeSummary(lib=WORK, dset=test, varlist=Var1 Var2 Var3, outfmt=SummaryFmt);

Ksharp · ‎08-13-2025

Check this for interact effect explanation: https://support.sas.com/kb/38/384.html Check this for relative risk: https://support.sas.com/kb/23/003.html

Ksharp · ‎08-13-2025

%let label1=XAXIS Label(*ESC*){sub 'L'} ; %let label2=YAXIS Label(*ESC*){sup 'L'} ; data _anno; length label $ 200; drawspace="layoutpercent"; function="text"; textweight="normal"; textsize=12;textcolor="black"; width=200; x1=50; y1=2.5;label="&label1."; output; x1=2.5; y1=50;rotate=90;label="&label2."; output; run; proc sgplot data=sashelp.class sganno=_anno; scatter x=weight y=height; xaxis label=' '; yaxis label=' '; run;

Ksharp · ‎08-12-2025

I should clarify that DARRTMAP dataset is applied to " group=year" , is unable to apply to "_region " ,so you can not change style of "_region ". Anyway a workaround way is using YAXISTABLE statement to mock YAXIS values: /* Here is my data */ data have; input code region $ 3-16 growth_rate year; datalines; 0 North America -5.92 2020 0 North America -2.18 2021 0 North America -7.47 2022 0 North America -2.17 2023 0 North America -9.64 2024 1 South America 0.36 2020 1 South America 2.69 2021 1 South America -8.13 2022 1 South America 1.22 2023 1 South America -1.34 2024 2 Asia -1.53 2020 2 Asia 2.18 2021 2 Asia 5.64 2022 2 Asia 6.80 2023 2 Asia 8.70 2024 3 Europe -9.68 2020 3 Europe 3.91 2021 3 Europe 5.48 2022 3 Europe 3.52 2023 3 Europe 3.69 2024 4 Africa 3.29 2020 4 Africa 7.97 2021 4 Africa 2.48 2022 4 Africa 7.36 2023 4 Africa 1.31 2024 5 Oceania 8.00 2020 5 Oceania -6.05 2021 5 Oceania 6.91 2022 5 Oceania -1.09 2023 5 Oceania -7.46 2024 ; run; /* First, I want to determine the order in which the values will be displayed */ proc sort data=have(where=(year=2024)) out=temp; by descending growth_rate ; run; data fmt; set temp(rename=(region=start)); retain fmtname 'fmt' type 'c'; length label $ 80; label=repeat(' ',_n_)||strip(start); keep fmtname type start label; run; proc format cntlin=fmt; run; /* Here, I determine the colors of the hbars*/ data have2; set have; _region=put(region,$fmt32.); run; /* Now, I create my range attribute map to customize the bar colors */ data dattrmap; infile cards truncover; input id $ value $ fillcolor $ textfamily &$40. textcolor $ textstyle $ textweight $; datalines; year1 2020 cxf7ebde year1 2021 cxefdbc6 year1 2022 cxe1ca9e year1 2023 cxd6ae6b year1 2024 cxb4771f year2 2020 cxdeebf7 year2 2021 cxc6dbef year2 2022 cx9ecae1 year2 2023 cx6baed6 year2 2024 cx1f77b4 code 0 . Times New Roman Uni red Italic Bold code 1 . Arial blue Normal Normal code 2 . Arial blue Normal Normal code 3 . Arial blue Normal Normal code 4 . Arial blue Normal Normal code 5 . Arial blue Normal Normal ; data have3; set have2; if strip(_region)="North America" then do;region1=_region;year1=year;end; else do;region2=_region;year2=year;end; run; proc sort data=have3; by _region year; run; /*Next, I try to apply different fonts to format the different levels of the '_region' variable and the values that appear next to each bar by including the option 'datalabelfitpolicy=none'*/ proc sgplot data=have3 dattrmap=dattrmap noautolegend; title "Growth rate"; hbarparm category=region1 response=growth_rate /group=year1 attrid=year1 groupdisplay=cluster clusterwidth=0.9 datalabel DATALABELFITPOLICY=NONE datalabelattrs=(size=8.5pt family='Times New Roman Uni' color=red style=Italic weight=bold) barwidth=1 nooutline transparency=0 dataskin=none; hbarparm category=region2 response=growth_rate /group=year2 attrid=year2 groupdisplay=cluster clusterwidth=0.9 datalabel DATALABELFITPOLICY=NONE datalabelattrs=(size=8.5pt family='Arial' color=blue style=Normal weight=Normal) barwidth=1 nooutline transparency=0 dataskin=none; xaxis display=(nolabel) min=-12 max=12; yaxis display=(nolabel novalues) colorbands=even ; yaxistable _region/nolabel y=_region position=left valueattrs=(size=10pt) textgroup=code textgroupid=code; run;

Ksharp · ‎08-12-2025

Check this: https://support.sas.com/kb/25/006.html https://support.sas.com/kb/69/654.html

Ksharp · ‎08-12-2025

Assuming I understood your question. Try this one: data have; input case_id control_id; datalines; 1 10 1 125 1 1850 1 928 3 125 3 8 3 1276 3 1011 3 10 4 10 4 1011 ; %let n=3; *each case_id has 3 control_id at most; data have2; set have; call streaminit(1234); random=rand('uniform'); run; proc sort data=have2; by case_id random; run; data want; if _n_=1 then do; if 0 then set have; declare hash h(); h.definekey('control_id'); h.definedone(); end; set have2; by case_id; if first.case_id then count=0; if h.check() ne 0 and count<&n. then do; count+1; h.add(); output; end; run;

Ksharp · ‎08-11-2025

Sure. Of course. You can scan all these LOG file and output the row started with ERROR or WARN. And attach the log file name together. %let path= c:\temp\log ; * the path stored LOG files; filename x "&path"; data log_file; did=dopen("x"); do i=1 to dnum(did); file=dread(did,i); if lowcase(scan(file,-1,'.'))='log' then do; log_file=cats("&path",'\',file); output; end; end; keep log_file; run; data want; set log_file; length fname filename $ 200; infile dummy filevar=log_file filename=fname end=last length=len; filename=fname; do while(not last); input log $varying200. len; if left(log) in: ('ERROR' 'WARN') then output; end; run;

Ksharp · ‎08-11-2025

You could save the LOG of sas programm into a file ,and check if there are any row started with ERROR or WARNING, if there were ,put an ERROR message into a macro variable. Here is an example: filename x temp; proc printto log=x new; run; data x; set sashelp.class; array x{*} _numeric_; y=x{10}; run; proc printto; run; %let rc=the code is good ; data _null_; infile x; input; if _infile_ in: ('ERROR' 'WARN') then call symputx('rc','the code is failed'); run; %put &=rc ;

Ksharp · ‎08-11-2025

I think ballardw was right. There could be some error data at OBS=177. For example: obs x1 x2 177 1 0 when you have formula "target=x1/x2;" ,that could lead to an error.

Ksharp · ‎08-11-2025

You could try to use CALL EXECUTE instead of "%include" . data MATCHING_TABLE_FUNC ; input Target $ Formula & $20.; cards; v1 x1 + x2 /x3 + x4 v2 x1 + x2 + x4 v3 x1 + x2 + x3 v4 x1 + x2 ; data TRNP_SERIES_LIST ; input date :date9. x1 x2 x3 x4; format date date9.; cards; 01jan2010 1 2 3 4 01jan2010 4 3 2 1 ; data MATCHING_TABLE_FUNC2; set MATCHING_TABLE_FUNC; function=catx('=',target,formula)||';'; keep function; run; data _null_; set MATCHING_TABLE_FUNC2 end=last; if _n_=1 then call execute('Data Formula_ ; set TRNP_SERIES_LIST ;'); call execute(strip(function)); if last then call execute('run;'); run;

Ksharp · ‎08-11-2025

Can you post dataset "TRNSP_SERIES_LIST " ,especially variable "Target" and "Formula" .

Ksharp · ‎08-09-2025

Sure. Of course. Just change the following line: %let path= c:\temp\a ; into %let path= %sysfunc(pathname(work)) ; BTW, you also need to change the path of RTF file, otherwise the RTF file would be stored at WORK library.

Ksharp · ‎08-08-2025

data a; input category $ group $ v lower upper p $; if group='Ole' then v1=v; if group='Ola' then v2=v; max=max(v,lower,upper,0); cards; Trig Ole -49.3 -52.3 -46.2 p<0.01 Trig Ola -64.3 -68.3 -58.2 p<0.01 APOC Ole -69.3 -75.3 -60.2 p<0.01 APOC Ola -79.3 -85.3 -68.2 p<0.01 VLDL Ole 49.3 42.3 55.2 p<0.01 VLDL Ola 40.3 35.3 48.2 p<0.01 ; proc sgplot data=a; styleattrs datacolors=(yellow gold) ; vbarparm category=category response=v1 /fillattrs=(color=white) groupdisplay=cluster group=group nooutline datalabel=v1 discreteoffset=-0.1; vbarparm category=category response=v2 /fillattrs=(color=white) groupdisplay=cluster group=group nooutline datalabel=v2 discreteoffset=0.1; scatter x=category y=max/group=group groupdisplay=cluster markerattrs=(size=0) datalabel=p datalabelpos=top datalabelattrs=(color=black) labelstrip; vbarparm category=category response=v / groupdisplay=cluster group=group nooutline limitlower=lower limitupper=upper LIMITATTRS=(color=black) name='x' ; keylegend 'x'/location=inside position=top; xaxis display=(nolabel); run; /*************************************/ data b; infile cards truncover; input group day v lower upper text $20.; cards; 1 0 200 180 210 1 1 150 120 160 1 2 120 110 128 1 3 110 102 118 1 6 109 101 112 1 10 105 100 110 2 0 210 190 220 2 1 160 130 170 2 2 110 120 138 2 3 120 112 128 2 6 119 111 122 2 10 115 110 120 . 8 124 . . TreatA|(n=26) . 8 102 . . TreatB|(n=26) ; proc sgplot data=b; series x=day y=v/group=group nomissinggroup; scatter x=day y=v/group=group nomissinggroup yerrorupper=upper yerrorlower=lower markerattrs=(symbol=circlefilled); xaxis type=linear values=(0 to 10 by 1) display=(nolabel noticks) valuesdisplay=('Baseline' 'Day28' 'Day56' 'Day84' '' '' 'Day252' '' '' '' 'Day308'); text x=day y=v text=text/strip contributeoffsets=none splitchar='|' splitpolicy=splitalways; run;

Ksharp · ‎08-07-2025

And Here is PROC REPORT solution. data report; input id label & $20. a & $20. b; cards; 1 Residual Statistics Observations 184 1 Residual Statistics Minimum 184 1 Residual Statistics Mean 184 1 Residual Statistics Maximum 184 1 Residual Statistics Std Dev 184 2 Fit Statistics Objective 184 2 Fit Statistics AIC 184 2 Fit Statistics AICC 184 2 Fit Statistics BIC 184 ; %let path= c:\temp\a ; *the path stored the picture of report and desired RTF file; ods _all_ close; options printerpath=png nonumber nodate papersize=(4.5cm 7cm) ; title; ods printer file="&path.\report.png" style=journal dpi=300; proc report data=report nowd noheader style(report)={frame=box}; columns id label a b; define id/order noprint; define label/order noprint; compute before label; if label ne lag(label) then n+1; if n=1 then len=0; else len=10; x=' '; line x $varying10. len; line label $40.; endcomp; run; ods printer close; options nodate nonumber; ods listing gpath="&path." style=htmlblue image_dpi=300; ods graphics /width=3in height=3in reset=index noborder imagename='FAS' outputfmt=png ; proc sgplot data=sashelp.class noautolegend ; reg x=height y=weight/cli clm; run; ods graphics /width=3in height=3in reset=index noborder imagename='PPS' outputfmt=png; proc sgplot data=sashelp.class; scatter x=age y=weight; run; ods graphics /width=3in height=3in reset=index noborder imagename='SS' outputfmt=png ; proc sgplot data=sashelp.class noautolegend ; ellipse x=height y=weight/group=sex; scatter x=height y=weight/group=sex; run; ods _all_ close; options papersize=A4 leftmargin=0in rightmargin=0in leftmargin=0in rightmargin=0in; ods rtf file="&path.\want.rtf" dpi=300 style=htmlblue; data x; x=' ';y=' ';output; x=' ';y=' ';output; run; title; proc report data=x nowd noheader style={outputwidth=100% }; column x y; define x/display; define y/display; compute y; n+1; if n=1 then do; call define('x','style','style={ preimage="&path\FAS.png" bordertopcolor=white borderbottomcolor=white borderrightcolor=white borderleftcolor=white}'); call define('y','style','style={ preimage="&path\PPS.png" bordertopcolor=white borderbottomcolor=white borderrightcolor=white borderleftcolor=white}'); end; if n=2 then do; call define('x','style','style={ preimage="&path\SS.png" bordertopcolor=white borderbottomcolor=white borderrightcolor=white borderleftcolor=white}'); call define('y','style','style={ preimage="&path\report.png" bordertopcolor=white borderbottomcolor=white borderrightcolor=white borderleftcolor=white}'); end; endcomp; run; ods rtf close;

Ksharp · ‎08-07-2025

You want this ? Here is GTL solution. Another is using PROC REPORT check the next post. data report; input id label & $20. a & $20. b; cards; 1 Residual Statistics Observations 184 1 Residual Statistics Minimum 184 1 Residual Statistics Mean 184 1 Residual Statistics Maximum 184 1 Residual Statistics Std Dev 184 2 Fit Statistics Objective 184 2 Fit Statistics AIC 184 2 Fit Statistics AICC 184 2 Fit Statistics BIC 184 ; %let path= c:\temp; *the path stored the picture of report and desired RTF file; ods _all_ close; options printerpath=png nonumber nodate papersize=(4.5cm 7cm) ; title; ods printer file="&path.\report.png" style=journal dpi=300; proc report data=report nowd noheader style(report)={frame=box}; columns id label a b; define id/order noprint; define label/order noprint; compute before label; if label ne lag(label) then n+1; if n=1 then len=0; else len=10; x=' '; line x $varying10. len; line label $40.; endcomp; run; ods printer close; %sganno data sganno; %SGIMAGE( IMAGE="&path.\report.png", ANCHOR="CENTER" , BORDER="FALSE", IMAGESCALE="FIT" , HEIGHT=100, HEIGHTUNIT="PERCENT", DRAWSPACE="LAYOUTPERCENT", X1=50, Y1=50,ID="BAR" ) run; proc template; define statgraph y2axis; begingraph; layout lattice /rows=2 columns=2; * columngutter=10 rowdatarange=union row2datarange=union; layout overlay; histogram height / scale=percent yaxis=y; densityplot height / normal(); endlayout; layout overlay; scatterplot x=weight y=height; endlayout; layout overlay; loessplot x=weight y=height; endlayout; layout overlay / walldisplay=none yaxisopts=( display=none ) xaxisopts=( display=none ) ; scatterplot x=weight y=height/markerattrs=(color=white); annotate / id="BAR"; endlayout; endlayout; endgraph; end; run; ods _all_ close; options papersize=A4 leftmargin=0in rightmargin=0in leftmargin=0in rightmargin=0in; ods rtf file="&path.\want.rtf" dpi=300 style=htmlblue; proc sgrender data=sashelp.class template=y2axis sganno=sganno; run; ods rtf close;

Online Status	Offline
Date Last Visited	17 hours ago