PROC IMPORT OUT= WORK.NICKY
            DATAFILE= "C:\Nicky\Asst Professorship\conferences\nacw workshop on gee glmm\telemetrydata.XLS" 
            DBMS=excel2000 REPLACE;
     *GETDELETED=NO;
RUN;
proc fsprint data = nicky;
run;

proc sort data = nicky;
by animal;
run;

TITLE1 'GLIMMIX model';
proc glimmix data = Nicky empirical;
        class animal habitat_1 habitat_2 habitat_3;
		model presabsence = logdist_7 logdist_16 habitat_1 habitat_2 habitat_3/solution dist=binomial link=logit CL;
        random intercept /subject=animal 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:\Nicky\Asst Professorship\conferences\nacw workshop on gee glmm\Glimmixconditional.xls"
	DBMS=EXCEL2000 REPLACE;
run;