The short answer is yes. You can read more here: https://www.sas.com/content/dam/SAS/support/en/sas-global-forum-proceedings/2018/1991-2018.pdf and/or by accessing the online SAS documentation for proc causal med.   Good luck! ... View more

Hi Rick, sorry to be bugging you about this. It seems there is no easy way to do a test for differences in predicted rank sum values.   So, what do you think about a weighted jonckheere-terpstra test instead? b is continuous non-normal - but from what I can find on the internet, the current version of proc freq with a jt option can handle ties.   proc freq data=have; table a*b / jt; weight z; run;   Thanks! ... View more

Please also see: https://stats.stackexchange.com/questions/67364/weighted-wilcoxon-ranksum-test#new-answer?newreg=9171852dfa82430baef6b8055609eb70 ... View more

That is: "the wilcoxon rank sum test should compare differences in predicted rank sums of b for the classification variable a based on the weight statement z (probability of a given c)". ... View more

  z is the probability of being in group a (my class variable in proc npar1way which has the values 0 or 1) given c.   proc logistic data=have namelen=32 plots=roc; model a(event='1') = c / stb; output out=want prob=z; run;   My question relates to this paper: https://www.nejm.org/doi/full/10.1056/nejmoa1110717   The authors state that they used wilcoxon rank sum tests to compare differences in two treatment regiments weighted by the probability of receiving one of the treatments (same as I am doing above). I don't understand the next step - how is it possible to weight the wilcoxon rank sum test by the probability of the observed level of the classification variable? ... View more

The problem was that proc pls renames the X variables to col1, col2 in the outmodel. I added rename to the outmodel statement - making it compatible with &xvars - and now the macro runs without errors.   proc pls data=data_a method=pls outmodel=est1(rename=(col1=age col2=male col3=female col4=nohtn col5=htn col6=hscrp)) lv=2; class female baseline_mi dm_all m5 htn inactive; model &yvars = &xvars; output out=outpls p=yhat1 yresidual=yres1 xresidual=xres1-xres4 xscore=xscr yscore=yscr stdy=stdy stdx=stdx h=h press=press t2=t2 xqres=xqres yqres=yqres; run; ... View more

Thanks for spending time on this. I really appreciate it.   Code:    data data_a; set data; if _N_ <= 15; n=_N_; run;   %global xvars yvars predname resname xscrname yscrname num_x num_y lv; %let xvars= age female htn hscrp; %let yvars= bmi; %let ypred=yhat1; %let yres=yres1; %let predname=yhat; %let resname=res; %let xscrname=xscr; %let yscrname=yscr; %let num_y=1; %let num_x=4;   proc pls data=data_a method=pls outmodel=est1 lv=2; class female baseline_mi dm_all m5 htn inactive; model &yvars = &xvars; output out=outpls p=yhat1 yresidual=yres1 xresidual=xres1-xres4 xscore=xscr yscore=yscr stdy=stdy stdx=stdx h=h press=press t2=t2 xqres=xqres yqres=yqres; run;   %let lv=2;   options mprint; %plot_scr(outpls); %plotxscr(outpls,max_lv=2); %get_wts(est1,dsxwts=xwts); /* see log below */ %plot_wt(xwts,max_lv=2); %getxload(est1,dsxload=xloads); %pltxload(xloads,max_lv=2); %get_bpls(est1,dsout=bpls); %get_vip(est1,dsvip=vip_data);   data eval; merge bpls vip_data; run; proc print data=eval; run;   Log for %get_wts with option mprint:   1864 %get_wts(est1,dsxwts=xwts); NOTE: There were 15 observations read from the data set WORK.DSOUT. NOTE: There were 15 observations read from the data set WORK.PLTANNO. NOTE: PROCEDURE GPLOT used (Total process time): real time 17.50 seconds cpu time 1.79 seconds MPRINT(GET_WTS): data xwts; MPRINT(GET_WTS): set est1(keep=_TYPE_ _LV_ age female htn hscrp); ERROR: The variable age in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable female in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable htn in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable hscrp in the DROP, KEEP, or RENAME list has never been referenced. MPRINT(GET_WTS): if _TYPE_='WB' then output; NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.XWTS may be incomplete. When this step was stopped there were 0 observations and 1 variables. WARNING: Data set WORK.XWTS was not replaced because this step was stopped. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.00 seconds MPRINT(GET_WTS): proc transpose data=xwts out=xwts; MPRINT(GET_WTS): run; NOTE: There were 0 observations read from the data set WORK.XWTS. NOTE: The data set WORK.XWTS has 1 observations and 1 variables. NOTE: PROCEDURE TRANSPOSE used (Total process time): real time 0.02 seconds cpu time 0.00 seconds MPRINT(GET_WTS): data xwts; MPRINT(GET_WTS): set xwts; MPRINT(GET_WTS): if _NAME_='_LV_' then delete; MPRINT(GET_WTS): n=_n_-1; MPRINT(GET_WTS): run; NOTE: There were 1 observations read from the data set WORK.XWTS. NOTE: The data set WORK.XWTS has 1 observations and 2 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.03 seconds MPRINT(GET_WTS): data xwts; MPRINT(GET_WTS): set xwts; MPRINT(GET_WTS): rename col1=w1; MPRINT(GET_WTS): run; WARNING: The variable col1 in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 1 observations read from the data set WORK.XWTS. NOTE: The data set WORK.XWTS has 1 observations and 2 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.00 seconds MPRINT(GET_WTS): data xwts; MPRINT(GET_WTS): set xwts; MPRINT(GET_WTS): rename col2=w2; MPRINT(GET_WTS): run; WARNING: The variable col2 in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 1 observations read from the data set WORK.XWTS. NOTE: The data set WORK.XWTS has 1 observations and 2 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.00 seconds   --------------------------------------   Looking at work.est1 (see attached) produced by proc pls, there are indeed no &xvars but instead col1-col6  -  it seems to me like something is not adding up.     ... View more

Ok, now I understand - Thanks.   Have you tried the pls macro? I copied the code from: https://support.sas.com/rnd/app/stat/papers/plsex.pdf   It seems to be working fine until I submit %get_wts(est1,dsxwts=xwts), which produces an error message:   NOTE: There were 15 observations read from the data set WORK.DSOUT. NOTE: There were 15 observations read from the data set WORK.PLTANNO. NOTE: PROCEDURE GPLOT used (Total process time): real time 14.57 seconds cpu time 1.04 seconds   ERROR: The variable age in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable dbp in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable female in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable htn in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable bmi in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable chol in the DROP, KEEP, or RENAME list has never been referenced. ERROR: The variable hscrp in the DROP, KEEP, or RENAME list has never been referenced. NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.XWTS may be incomplete. When this step was stopped there were 0 observations and 1 variables. WARNING: Data set WORK.XWTS was not replaced because this step was stopped. NOTE: DATA statement used (Total process time): real time 0.03 seconds cpu time 0.01 seconds   The macros %get_bpls and %get_vip also produce error messages:   2041 %get_bpls(dsoutmod,dsout=bpls); NOTE: There were 16 observations read from the data set WORK.DSOUTMOD. NOTE: The data set WORK.EST_WB has 4 observations and 24 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.00 seconds   NOTE: There were 16 observations read from the data set WORK.DSOUTMOD. NOTE: The data set WORK.EST_PQ has 4 observations and 24 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.03 seconds NOTE: IML Ready ERROR: AGE is not in the scope of variables for the data set.     */-----------------------------------------------------------------------*/   As far as I can see in the macro, there is no renaming of x-variables so I don't understand the problem. Could it be related to the number of characters in the names of the x variables? Here is the code I submitted:   data data_a; set data; if _N_ <= 15; n=_N_; run; %global xvars yvars predname resname xscrname yscrname num_x num_y lv; %global xvars yvars; %let xvars= age dbp female htn bmi chol hscrp; %let yvars= sbp; %let ypred=yhat1; %let yres=yres1; %let predname=yhat; %let resname=res; %let xscrname=xscr; %let yscrname=yscr; %let num_y=1; %let num_x=15;   proc pls data=data_a method=pls outmodel=dsoutmod lv=4; class female htn; model &yvars = &xvars; output out=outpls p=yhat1 yresidual=yres1 xresidual=xres1-xres15 xscore=xscr yscore=yscr stdy=stdy stdx=stdx h=h press=press t2=t2 xqres=xqres yqres=yqres; run;   %let lv=4; %plot_scr(outpls); %plotxscr(outpls,max_lv=4); %get_wts(est1,dsxwts=xwts); %plot_wt(xwts,max_lv=4); %getxload(est1,dsxload=xloads); %pltxload(xloads,max_lv=4); %get_bpls(dsoutmod,dsout=bpls); %get_vip(dsoutmod,dsvip=vip_data);   data eval; merge bpls vip_data; run; proc print data=eval; run;   ... View more