BookmarkSubscribeRSS Feed
Ziquilin
Calcite | Level 5

Dear SAS users,

We have some doubts with a longitudinal analysis using GLIMMIX and hope someone could advise us. We have daily movement data (Euclidean distances) of 45 individuals of a large bird of prey monitored during seven years which, in total, summed 18574 observations. Birds are evenly distributed in three different regions and no bird moved between regions. We are interested to know if daily movements differed among regions after controlling for bird age. The model we constructed was:


PROC glimmix DATA=COMP.PREDISP method=laplace;

CLASS birdidentity region year;

model distdiar= region age region*age /dist=XXX s;

random int / sub=birdidentity;

random int /sub=year;

LSMEANS REGION/DIFF CL;

run;


Model results (Type III tests of fixed effects) show 15298 denominator degrees of freedom. The model is considering both within and between individual variation, so we are wondering if this procedure “artificially” inflates (in our particular case) degrees of freedom to the point that results could be wrong. Any suggestion?


Thank you very much in advance


David

3 REPLIES 3
1zmm
Quartz | Level 8

PROC GLIMMIX defaults to the containment method to determine denominator degrees of freedom for models with only random effects like yours.  You should consider first sorting your data by your random effects--BIRDIDENTITY and YEAR--so that all observations for a particular bird across the observation period are ordered together.  Then consider using another method for determining denominator degrees of freedom (for example, specifying the MODEL statement option, DDFM=KENWARDROGER or DDFM=SATTERTHWAITE).

SteveDenham
Jade | Level 19

Watch out though for these two options for DDFM= when the method is LAPLACE or QUAD--they are not compatible (see Table 41.4 in the GLIMMIX documentation).  Try DDFM=BW (between/within).

Steve Denham

Ziquilin
Calcite | Level 5

Thanks a lot. Both solutions seem to work, although the BW option with LAPLACE method needs a nested design of random effects -birdidentity(region)- to correct the d.f. As we are not especially interested in obtaining very accurate estimates or AIC values the DDFM=KR (or SAT) with residual ML or PL is also a good option.

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 3 replies
  • 2151 views
  • 8 likes
  • 3 in conversation