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

proc sort data = nicky;
by animal;
run;

TITLE1 'GEE model';
proc genmod data = nicky descending;
        class animal habitat_1 habitat_2 habitat_3;
		model presabsence = logdist_7 logdist_16 habitat_1 habitat_2 habitat_3/ dist = binomial corrb;
        repeated subject = animal / corr=CS modelse;
		output out=saskresids predicted=inverselogit reschi=pearsresid stdreschi=stpearsresid STDXBETA=stdxbeta xbeta=logit; 
Run;

Proc export data=saskhomemerge3resids
	OUTFILE="C:\Nicky\Asst Professorship\conferences\nacw workshop on gee glmm\GEE_CS.xls"
	DBMS=EXCEL2000 REPLACE;
run;