%macro clem(var1=,var2=);

proc glimmix data=main.Final_year4;

class status;

model &var1=status

/dist=normal /*binomial*/ link=identity/*logit*/ s ddfm=kr;

lsmeans status / ilink diff;

ods output

/*ModelInfo=x1var1 ParameterEstimates=x2var1 CovParms=x3var1*/

Diffs=DIF_RESULT1 LSMeans=LS1;

run;

data &var1.;

set LS1 DIF_RESULT1;

outcome_name="&var1";

v="vs.";

if Mu ne . then analysis_type="Descriptive Stats";

if Mu = . then do;

analysis_type="Anova";

if Probt < 0.05 then sig="* ";

if Probt < 0.01 then sig="** ";

if Probt < 0.001 then sig="***";

if Probt =. then sig="";

end;

drop effect;

run;

data &var2;

retain outcome_name analysis_type status v _status;

set &var1.;

run;

%mend clem;

%clem(var1=year4_outcome1,var2=a);

%clem(var1=year4_outcome2,var2=b);

%clem(var1=year4_outcome3,var2=c);

%clem(var1=year4_outcome4,var2=d);

data em;

run;

data allresults;

set a em b em c em d ;

run;