PROC IMPORT OUT= WORK.YOURNAME
            DATAFILE= "C:\FOLDERS\FILENAME.XLS" 
            DBMS=excel2000 REPLACE;
     *GETDELETED=NO;
RUN;
proc fsprint data = YOURNAME;
run;

proc sort data = YOURNAME;
by GROUPINGVARIABLE;
run;

TITLE1 'GLIMMIX model';
proc glimmix data = YOURNAME empirical;
        class GROUPINGVARIABLE INDEPENDENT1 INDEPENDENT2 INDEPENDENT3;
		*ABOVE ONLY INCLUDES CATEGORICAL VARIABLES, BUT MUST INCLUDE GROUPING VARIABLE;
		model RESPONSE = INDEPENDENT1 INDEPENDENT2 INDEPENDENT3 INDEPENDENT4/solution dist=binomial link=logit CL;
        random intercept /subject=GROUPINGVARIABLE TYPE=vc;
		nloptions tech=nrridg;
		*uses newton-raphson with ridging optimization technique, previous line may not be necessary for many data sets;
		Title 'Glimmix model';
output out=Glimmixconditional pred = p Pearson=PEARSRESID UCL=UPPER LCL=LOWER;
Run;
Proc export data=Glimmixconditional
	OUTFILE="C:\FOLDERS\Glimmixconditional.xls"
	DBMS=EXCEL2000 REPLACE;
run;