options nonotes nosource nosource2 errors=0; data temp; infile '****\testing.txt' ; /*path for data*/ input X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 Z Y index; run; data saveres; Descr='-2 Log Likelihood'; value=0; nullstatus=0; rep=0; run; libname mylib 'C:***'; /*path for library*/ %include '***\glimmix.sas'; /*path for glimmix.sas*/ %MACRO model1(number); %let count=0; %do %while(&count<=&number); data use; set temp; where index=&count; run; %glimmix(data=use,procopt=method=ml, options=NOPRINT, stmts=%str( model Y=Z; random X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 / type=toep(1) s;), error=binomial, link=logit, out=fitted); data tosave; set _fitstats; nullstatus=0; rep=&count; where descr="-2 Log Likelihood"; run; data saveres; set saveres tosave; run; run; %glimmix(data=use,procopt=method=ml, options=NOPRINT, stmts=%str( model Y= ;), error=binomial, link=logit, out=fitted); data tosave; set _fitstats; nullstatus=1; rep=&count; where descr="-2 Log Likelihood"; run; data saveres; set saveres tosave; run; %let count=&count+1; %end; data mylib.save; set saveres; run; %MEND model1; %model1 (number=1000); /*aligned with nper.lrt specified in SASmatrix in R*/ quit; options notes source source2 errors=20; proc export data=mylib.save outfile="****\testing.xls" dbms=excel replace; run;