About Ksharp

Ksharp · ‎08-20-2025

data paths; length path $200; input path $; datalines; /dwh_actuariat/sasprocess/prodmens/HYFIprep/output /dwh_actuariat/sasprocess/prodmens/HYFIp/output ; run; %macro process; proc sql noprint; select path into :path_list separated by ' ' from paths; quit; %put &=path_list; %let i = 1; %do %while (%qscan(&path_list, &i. ,%str( )) ne ); %let current_path = %scan(&path_list, &i. , %str( )); %put &=current_path; %let i = %eval(&i. + 1); %end; %mend process; %process;

Ksharp · ‎08-20-2025

I did not see an option you could do that either. The workaround way I could image is using PROC SGPLOT to make a picture for LEGEND and insert this picture into PROC SGRENDER by rotating it 90 degree. But that would cost you a lot of time and you need to make legend be consistent with original graph. Here is an example: %let path=%sysfunc(pathname(work)); *the path stored the picture ; *Make a picture for Legend; data x; input label $ x y s_male_x s_male_y male $ male_x male_y s_female_x s_female_y female $ female_x female_y; cards; SEX 0 0 1 0 Male 1.5 0 2.5 0 Female 3.2 0 ; run; title; ods listing gpath="&path." style=htmlblue image_dpi=300; ods graphics/width=250px height=40px reset=index imagename='legend' outputfmt=png noborder; proc sgplot data=x noautolegend pad=0; *noborder; scatter x=x y=y/markerchar=label markercharattrs=(size=20); scatter x=s_male_x y=s_male_y/markerattrs=graphdata1(symbol=squarefilled size=20); scatter x=male_x y=male_y/markerchar=male markercharattrs=(size=20); scatter x=s_female_x y=s_female_y/markerattrs=graphdata2(symbol=squarefilled size=20); scatter x=female_x y=female_y/markerchar=female markercharattrs=(size=20); xaxis display=none; yaxis display=none; run; %sganno data sganno; %SGIMAGE( IMAGE="&path.\legend1.png", ANCHOR="CENTER" , IMAGESCALE="FITHEIGHT", ROTATE=90, HEIGHT=10, HEIGHTUNIT= "PERCENT", LAYER="FRONT", DRAWSPACE="WALLPERCENT", X1=5, Y1=65,ID="BAR" ) run; proc template; define statgraph y2axis; begingraph; layout overlay ; barchart category=age response=weight/group=sex groupdisplay=cluster ; annotate / id="BAR"; endlayout; endgraph; end; run; ods graphics/reset=all; proc sgrender data=sashelp.class template=y2axis sganno=sganno; run;

Ksharp · ‎08-19-2025

You could faster your code via BY statement. For example: DATA mcont1; SET pairedcontinuous; WHERE Variables='cont1'; RUN; proc sort data=mcont1; by Variables; run; /*ODS OUTPUT ..... <---this could save the output as a dataset*/ PROC TTEST data= mcont1; by Variables; PAIRED _0*_1 ; RUN; DATA mcat1; SET pairedcategorical; WHERE Variables='cat1'; RUN; proc sort data=mcat1; by Variables; run; /*ODS OUTPUT ..... <---this could save the output as a dataset*/ PROC FREQ data=mcat1 ; by Variables; TABLE _0*_1; EXACT mcnem; RUN;

Ksharp · ‎08-19-2025

For this scenario , you could match this condition by operator BETWEEN AND. proc sql; create table CLIN_AE (SUBJID char(11), AEDECOD char(100), AESTDAT num format=date9.); insert into clin_ae (subjid, aedecod, aestdat) values ('101-001-001','Rash','30MAY2025'd) values ('101-001-001','Rash','18JUN2025'd); create table SAFE_AE (USUBJID char(11), DECOD char(100), STDTC num format=date9.); insert into safe_ae (usubjid, decod, stdtc) values ('101-001-001','Rash','30MAY2025'd) values ('101-001-001','Rash','19JUN2025'd); create table compare as select a.*, b.* from clin_ae as a full join safe_ae as b on a.subjid=b.usubjid & a.aedecod=b.decod & b.stdtc between a.aestdat and a.aestdat+1; quit; data compare; set compare; if aestdat^=stdtc then DISCREP='Start Dates not Matching'; run;

Ksharp · ‎08-18-2025

And PROC COPY could get you faster: proc copy in=work out=Server1; select DiagnosisTable; run;

Ksharp · ‎08-18-2025

/* 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; 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 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'*/ /*However, I did not find a way to exclude to exclude missing values in the keylegend statement*/ /* As I want to just plot the legend for "North America", I have tried to do this by using NAME option only for "year1" of the attribute map and the option EXCLUDE=(" ") in keylegend. However, I still see a gray square with "." (dot) as a legend for the squares representing the hbarparam of "year1". I also checked the SAS manual and there is no "nomissinggroup" option for hbarparm. I also tried using EXCLUDE = (""), EXCLUDE = ("."), EXCLUDE = (.), or EXCLUDE = (missing), but nothing seems to work. /* Here is what I have tried*/ proc sgplot data=have3 dattrmap=dattrmap noautolegend; title "Growth rate"; /*Here I use NAME to create the legend for year1*/ hbarparm category=region1 response=growth_rate /group=year1 attrid=year1 NAME="legend_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; /*Here I create the hbarparam for year2, but without the NAME option as I do not want to display values for year2 in the legend*/ 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; keylegend "legend_year1"/title="" position=bottom noborder EXCLUDE=(".") down=1 titleattrs=(size=8.5pt) valueattrs=(size=8.5pt) autoitemsize; run;

Ksharp · ‎08-18-2025

Another solution is using SAS/IML , if you have it. data have; input Have1-Have3 Expect :$32.; cards; 100 200 250 Have1 200 100 275 Have2 300 900 100 Have3 250 25 120 Have2 50 8 65 Have2 400 700 200 Have3 ; proc iml; use have; read all var _num_ into x[c=vname]; close; want=vname[x[,>:<]]; create want from x want[c=(vname||'want')]; append from x want; close; quit;

Ksharp · ‎08-18-2025

Then use the right most variable P30: compute P30; if dif_p10>0.05 then call define ('P10',"style", "style={background=red}");

Ksharp · ‎08-18-2025

https://communities.sas.com/t5/Statistical-Procedures/Dealing-with-nonproportionality-in-a-Cox-PH-Model/m-p/968428#M48690

Ksharp · ‎08-18-2025

"I didnt see any information to read about this rule " Why not make an example to test this rule ?

Ksharp · ‎08-18-2025

Yes. You are right. As long as the variable in COMPUTE compute dif_p90; is at the right side of the variable you refer to (e.x. p90)

Ksharp · ‎08-17-2025

PROC REPORT scan table from left to right. So when you want to refer to left column base on right column, you need use right column within COMPUTE block. data have; input P_IND_Prob mon P10 P20 P30 P40 P50 P60 P70 P80 P90; cards; 0 202507 0.1 0.15 0.25 0.31 0.49 0.46 0.6 0.75 0.8 0.97 0 202506 0.08 0.13 0.14 0.44 0.52 0.65 0.74 0.82 0.91 0.94 ; run; proc sort data=have; by mon; run; data have2; set have; by mon; dif_p10 =abs(p10-lag(p10)); dif_p20 =abs(p20-lag(p20)); dif_p30 =abs(p30-lag(p30)); dif_p40 =abs(p40-lag(p40)); dif_p50 =abs(p50-lag(p50)); dif_p60 =abs(p60-lag(p60)); dif_p70 =abs(p70-lag(p70)); dif_p80 =abs(p80-lag(p80)); dif_p90 =abs(p90-lag(p90)); run; Proc format ; value highlight_S4_Fmt 0.05<-high = 'red' ; Run; proc report data=have2 missing nowd style(report)={frame=box font_size=8pt bordercolor=black borderwidth=2px} ; column mon P10 P20 P30 P40 P50 P60 P70 P80 P90 dif_p10 dif_p20 dif_p30 dif_p40 dif_p50 dif_p60 dif_p70 dif_p80 dif_p90 ; define mon/DISPLAY ; define P10/DISPLAY f=percent10.4 ; define P20/DISPLAY f=percent10.4; define P30/DISPLAY f=percent10.4; define P40/DISPLAY f=percent10.4; define P50/DISPLAY f=percent10.4; define P60/DISPLAY f=percent10.4; define P70/DISPLAY f=percent10.4; define P80/DISPLAY f=percent10.4; define P90/DISPLAY f=percent10.4; define dif_p10/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p20/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p30/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p40/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p50/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p60/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p70/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p80/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; define dif_p90/DISPLAY f=percent10.4 style=[background=highlight_S4_Fmt.]; compute dif_p90; if dif_p10>0.05 then call define ('P10',"style", "style={background=red}"); if dif_p20>0.05 then call define ('P20',"style", "style={background=red}"); if dif_p30>0.05 then call define ('P30',"style", "style={background=red}"); if dif_p40>0.05 then call define ('P40',"style", "style={background=red}"); if dif_p50>0.05 then call define ('P50',"style", "style={background=red}"); if dif_p60>0.05 then call define ('P60',"style", "style={background=red}"); if dif_p70>0.05 then call define ('P70',"style", "style={background=red}"); if dif_p80>0.05 then call define ('P80',"style", "style={background=red}"); if dif_p90>0.05 then call define ('P90',"style", "style={background=red}"); endcomp; Run;

Ksharp · ‎08-14-2025

Did you try PROC APPEND ? proc append base=spe.limits data=limits5 force; run;

Ksharp · ‎08-14-2025

FreelanceReinh ， Good Point. Learned it happy.

Ksharp · ‎08-14-2025

If you really need to make your code worked, just add %UNQUOTE into it to get rid of delta character generated by %str(). data test; input id var1-var3; format var3 f4.2; cards; 1 1 2 3 ; %macro MakeSummary(lib=WORK, dset=, varlist=, outfmt=SummaryFmt); %LET MultRespVariablesList=%str(%')%sysfunc(tranwrd(%upcase(&varlist.),%str( ),%str(',')))%str(%'); %PUT &=MultRespVariablesList; 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 upcase(name) IN (%unquote(&MultRespVariablesList.)); quit; %mend; %MakeSummary(lib=WORK, dset=test, varlist=Var1 Var2 Var3, outfmt=SummaryFmt);

Online Status	Offline
Date Last Visited	15 hours ago