Ksharp Tracker

Re: proc pls : saving the model

Ksharp — Sat, 16 May 2026 07:12:04 GMT

https://communities.sas.com/t5/Statistical-Procedures/score-for-PLS-procedure-issue/m-p/676606

The only thing you need to do is using SOLUTION option to save parameters estimator and score test data as PROC GLM .

Re: A SAS Macro to Use ODS to Direct Graphics Output to a Specified Folder and Avoid the WARNING: GP

Ksharp — Sat, 16 May 2026 07:06:49 GMT

Did you used:

ods listing gpath='c:\temp';

https://blogs.sas.com/content/sgf/2015/11/20/how-to-increase-the-resolution-of-your-sas-graphics-output/

https://blogs.sas.com/content/graphicallyspeaking/2017/06/08/getting-started-output-formats/

And SAS_EXECFILEPATH macro variable is not working in batch mode (I tested), you need to set parameter

-set SAS_EXECFILEPATH "c:\temp\"

in batch process file.

Re: Is It Acceptable to Use Different Box-Cox Transformations for Each Variable?

Ksharp — Wed, 13 May 2026 06:52:53 GMT

Sorry. Your question is too tough for me.
If your requirement is regarding FORECASTING , you could post your question at Forecasting Forum:

https://communities.sas.com/t5/SAS-Forecasting-and-Econometrics/bd-p/forecasting_econometrics

Re: Is It Acceptable to Use Different Box-Cox Transformations for Each Variable?

Ksharp — Wed, 13 May 2026 06:37:06 GMT

My two cents:
I think it is OK.
Box-Cox Transformations is just to transform non-normal distribution into normal distribution. It has nothing to do with other variables. But the explanation of parameter estimator of that variable is another story !

Re: Separating the same dates in a graph

Ksharp — Tue, 12 May 2026 07:15:19 GMT

Why are you using VLINE ? why not use SCATTER?

VLINE would get MEAN of obs which have the same date.

And you'd post a desired graph to explain what you are looking for ?

libname x  v9 'C:\Users\xiakeshan\Documents\Downloads';

data reset;
 set x.reset;
run;

proc sgplot data=reset ;
scatter x=Date y=Result/jitter group=Structure ;
xaxis interval=Auto type=time label='Sampled Date' LABELATTRS=(weight=bold size=9)
VALUEATTRS=(family=arial size=9 style=normal weight=normal );
yaxis values=(0 to 100 by 5) LABELATTRS=(weight=bold size=9)
VALUEATTRS=(family=arial size=9 style=normal weight=normal );;
refline 29 / axis=y label=('Expected Result') lineattrs=(color=red pattern=dash)
labelattrs=(size=12) labelloc=outside;
title1 j=c " design Structure" bold;
run;

Re: Transposing a long dataset to wide with multiple variable types

Ksharp — Wed, 13 May 2026 06:23:20 GMT

Check the MERGE skill proposed by me, Arthur.T and Matt :

http://support.sas.com/resources/papers/proceedings15/2785-2015.pdf

And could you post your sample data and the desired output to explain this question better?

data have;
input id char $ num date :date9. time :time. ;
format date date9. time time.;
cards;
1 ss 23 1jan2010 09:20:21
1 ss 24 2jan2010 09:20:22
1 ss 25 3jan2010 09:20:23
2 ss 26 4jan2010 09:20:24
2 ss 27 5jan2010 09:20:25
;



data temp;
 set have;
 by id;
 if first.id then n=0;
 n+1;
run;
proc sql noprint;
select distinct(catt('temp(where=(n=',n,') rename=(
 char=char_',n,' num=num_',n,' date=date_',n,' time=time_',n,'))'))
 into :merge separated by ' '
   from temp;
quit;
data want;
merge &merge ;
by id;
drop n;
run;

And a more simple and convenient way is using proc summary+idgroup, but it has limit that each group have 100 obs at most.

data have;
input id char $ num date :date9. time :time. ;
format date date9. time time.;
cards;
1 ss 23 1jan2010 09:20:21
1 ss 24 2jan2010 09:20:22
1 ss 25 3jan2010 09:20:23
2 ss 26 4jan2010 09:20:24
2 ss 27 5jan2010 09:20:25
;
proc sql noprint;
select max(n) into :n from
(select count(*) as n  from have group by id);
quit;

proc summary data=have nway;
class id;
output out=want idgroup( out[&n] (char num date time)=);
run;

Re: CV for ridge regression

Ksharp — Mon, 11 May 2026 10:01:38 GMT

Rick,
OP used the method you are mentioned, check the first post.
OP want " Ridge shrinkage criterion" , but PROC GLMSELECT didn't display it .

Re: CV for ridge regression

Ksharp — Mon, 11 May 2026 09:22:18 GMT

Sorry. I can't help you. It seems a tough task.
Maybe @Rick_SAS @StatDave could give you a hint.

Re: CV for ridge regression

Ksharp — Mon, 11 May 2026 08:33:38 GMT

OK. Check the following code was what you are looking for ?

data acetyl;
input x1-x4 @@;
x1x2 = x1 * x2;
x1x1 = x1 * x1;
datalines;
1300  7.5 0.012 49   1300  9   0.012  50.2 1300 11 0.0115 50.5
1300 13.5 0.013 48.5 1300 17   0.0135 47.5 1300 23 0.012  44.5
1200  5.3 0.04  28   1200  7.5 0.038  31.5 1200 11 0.032  34.5
1200 13.5 0.026 35   1200 17   0.034  38   1200 23 0.041  38.5
1100  5.3 0.084 15   1100  7.5 0.098  17   1100 11 0.092  20.5
1100 17   0.086 29.5
;
 


/****** K-Fold CV ****/
%macro k_fold_cv(k=);
filename score temp;

ods select none;

proc surveyselect data=acetyl group=&k out=have seed=123;
run;

%do i=1 %to &k ;
data training;
 set have(where=(groupid ne &i)) ;
run;
data test;
 set have(where=(groupid eq &i));
run;

proc reg data=training outest=b ridge=0.02 noprint;
   model x4=x1 x2 x3 x1x2 x1x1;
   code file=score residual ;
quit;

/*Score test data*/
data score;
set test;
%include score;
run;
/*Calculate PRESS*/
proc sql;
create table press as
select uss(R_x4)/count(*) as press from score;
quit;

data score&i;
 merge b(where=(_TYPE_='RIDGE') keep=_TYPE_ _RIDGE_) press;
 retain id &i ;
run;
%end;
data k_fold_cv;
 set score1-score&k;
run;

ods select all;
%mend;

%k_fold_cv(k=3)

Re: CV for ridge regression

Ksharp — Mon, 11 May 2026 06:36:45 GMT

" I cannot perform automatic Ridge cross-validation in SAS,"
I think so. But maybe some expert of sas could know something new .

"Is there an existing macro that performs K-folds Ridge cross-validation?"
You could replace PROC LOGISTIC in my code with PROC REG + RIDGE= option and get what you need .

Re: CV for ridge regression

Ksharp — Sun, 10 May 2026 02:16:21 GMT

You could use RIDGE= option of PROC REG to perform ridge regression, check Rick's blog:
https://blogs.sas.com/content/iml/2013/03/20/compute-ridge-regression.html

About CV you could write your own code , check my code:
https://communities.sas.com/t5/SAS-Data-Science/10-fold-corss-validation/m-p/366193
https://communities.sas.com/t5/SAS-IML-Software-and-Matrix/How-to-do-k-fold-CV-with-replacements-replication/m-p/303938

Re: How do organize he data to use ANOVA

Ksharp — Sat, 09 May 2026 06:07:28 GMT

Post some your sample data to better explain your question.
ANOVA is only for continuous/consequence variable, for ordinary check WILCOXON test which is non-parameter version of ANOVA (proc npar1way + wilcoxon option).

Re: Cross-validation for 10 groups

Ksharp — Mon, 04 May 2026 01:40:06 GMT

I think "cv=split(10)" is good enough.

and I quoted:

CV=TESTSET(SAS-data-set)
specifies a test set of observations to be used for validation (formally, this is called “test set validation”
rather than “cross validation”).

or using "CV=ONE" or "CV=BLOCK(10)"

Re: Issue with proc format

Ksharp — Sat, 02 May 2026 03:02:34 GMT

Make two version of code, one for varieties (I49L and Hass) , another for varieties (others), and make a macro to wrap them up .

Re: SAS fcopy failing to copy .xls file (creates corrupted 1KB file)

Ksharp — Sat, 02 May 2026 02:39:11 GMT

You also could try this code:

%let in1=H:\folder1\file1.xls;
%let ou1=H:\folder2\file1.xls;

data _null_;
rc=rename("&in1","&ou1",'file');
run;

Re: How do I delete observations with dates occurring within another date for the same ID?

Ksharp — Sat, 02 May 2026 01:57:50 GMT

You want to keep the obs which have the same AdmitDate?

data have;
infile cards truncover;
input ID	(AdmitDate	DischargeDate) (: mmddyy10.);
format AdmitDate	DischargeDate mmddyy10.;
cards;
1	1/1/2025	2/4/2025
1	1/1/2025	4/4/2025
1	1/3/2025	1/3/2025
1	1/17/2025	1/17/2025
2	3/4/2025	3/4/2025
3	6/7/2025	6/9/2025
3	6/8/2025	6/8/2025
;
/*keep dates that are the same as the admit date*/
proc sql;
create table dup as
select * from have 
 group by id,AdmitDate
  having count(*)>1;
quit;


data temp;
 set have;
 do date=AdmitDate	to DischargeDate;
  output;
 end;
 keep id date;
 format date mmddyy10.;
run;
proc sort data=temp nodupkey;
by id date;
run;
data temp2;
 set temp;
 by id;
 if first.id or dif(date) ne 1 then group+1;
run;
proc summary data=temp2;
by id group;
var date;
output out=temp3(drop=_: group) min=AdmitDate max=DischargeDate ;
run;


/*combine them together*/
proc sql;
create table want as
select * from dup
union
select * from temp3;
quit;

Re: How do I delete observations with dates occurring within another date for the same ID?

Ksharp — Sat, 02 May 2026 01:39:01 GMT

Tom,
As I said , the most important thing I care about is readable for code, not length, not performance . not others ....
Unless the running time of code is WAY too much I can't accept, in that case I would consider about enhancing the performance of code .
Your code is not easy to read than mine, don't you think so ?

Re: How do I delete observations with dates occurring within another date for the same ID?

Ksharp — Fri, 01 May 2026 13:21:16 GMT

Tom,
I just want to keep my code simple and readable.
I really don't care about the code is long or short .
I think the short code would bring side effect . The readable code is the most important thing.

Re: How do I delete observations with dates occurring within another date for the same ID?

Ksharp — Fri, 01 May 2026 02:05:13 GMT

data have;
infile cards truncover;
input ID	(AdmitDate	DischargeDate) (: mmddyy10.);
format AdmitDate	DischargeDate mmddyy10.;
cards;
1	1/1/2025	2/4/2025
1	1/3/2025	1/3/2025
1	1/17/2025	1/17/2025
2	3/4/2025	3/4/2025
3	6/7/2025	6/9/2025
3	6/8/2025	6/8/2025
;

data temp;
 set have;
 do date=AdmitDate	to DischargeDate;
  output;
 end;
 keep id date;
 format date mmddyy10.;
run;
proc sort data=temp nodupkey;
by id date;
run;
data temp2;
 set temp;
 by id;
 if first.id or dif(date) ne 1 then group+1;
run;
proc summary data=temp2;
by id group;
var date;
output out=want(drop=_: group) min=AdmitDate max=DischargeDate ;
run;

Re: How to plot "A" outside of the plotting area in base graphics?

Ksharp — Tue, 28 Apr 2026 06:53:40 GMT

/*****************************/
%let label1=A(*ESC*){sup 'mw -1'}  ;
data _anno;
length label $ 200;
 drawspace="graphpercent"; function="text"; textweight="normal"; textsize=12;textcolor="black"; width=200;
 x1=5; y1=95;label="&label1."; output;
 run;
proc sgplot data=sashelp.class sganno=_anno;   
scatter x=weight y=height;   
xaxis label=' ';
yaxis label=' ';
run;