BookmarkSubscribeRSS Feed
lidove
Calcite | Level 5

Dear  All,

I am a beginner of SAS. I used glimmix procedure to model how a binary response variable was affected by several interested factors. The output gave inconsistent results of the significance test of one factor, "gender", from the "Solutions for Fixed Effects" and "Type III Tests of Fixed Effects" (please see the attached results). Could anybody tell me:  1) what's the reason for my results? Is it due to the nature of my dataset which is unbalanced?  2) And under this situation, how do I intepret my results?  which result should I belive, "Solutions for Fixed Effects" or "Type III Tests of Fixed Effects"?

Thanks for the help.

Solutions for Fixed Effects

Effect        Nestlingsex       Estimate    Standard Error        DF          t            P

Intercept                             2.7637         0.3415                 446.8    8.09     <0.0001

……

Gender            female       -0.8945         0.3448                  9237   -2.59     0.0095

Gender             male              0

……

Type III Tests of Fixed Effects

Effect         Num DF   Den DF   F Value       Pr > F

……

Gender          1              9237      0.04           0.8491

……

3 REPLIES 3
SteveDenham
Jade | Level 19

That seems odd.  Can you share your PROC GLIMMIX code, and the entirety of your Type3 output?  I have a hunch you have some interactions with gender that might be significant, but I'm really just guessing.

Steve Denham

lidove
Calcite | Level 5

Thank you so much, Steve. The work is to examine what kinds of factors would affect whether a nestling bird would be fed during each feeding event.  The SAS code and type 3 output are as follows:

proc glimmixdata = sasuser.EACHALLOCATION  ;

class  adult  nestling gender  nest position  begging CB  No_feeding_event;

title GLIMMIX Each Feeding preferenc;

model F_N(event='1') = adult  gender  position begging  Tsize  CB adult*gender  adult*position  adult*begging adult*Tsize  adult*CB    gender*position        gender*begging    gender*CB  gender*Tsize position*begging   position*CB  position*Tsize   begging*CB begging*Tsize   Tsize*CB NO_begging_nestlings        /solution  ddfm=kr link = logit   dist = binomial ;

random intercept  no_feeding_event / subject=nest;

run;


Type III Tests of Fixed Effects
EffectNum DFDen DFF ValuePr > F
Adult269251.690.1845
Gender192370.040.8491
Position19237844.43<.0001
Begging19237338.44<.0001
Tsize192372.060.1516
CB132.610.220.6394
Adult*Gender292373.660.0259
Adult*Position292370.170.8460
Adult*Begging292374.290.0138
Tsize*Adult292374.480.0114
Adult*CB178070.180.6747
Gender*Position192372.610.1060
Gender*Begging192372.120.1454
Gender*CB154631.470.2253
Tsize*Gender192370.010.9216
Position*Begging192379.410.0022
Position*CB192371.810.1785
Tsize*Position192371.430.2313
Begging*CB192372.190.1390
Tsize*Begging192370.320.5746
Tsize*CB192370.230.6333
No_begging_nestlings1218545.46<.0001

SteveDenham
Jade | Level 19

It doesn't depend on the significance of the interaction--the mere inclusion of the interaction in the model will lead to differences between the solution values and the F test values.  I found this to be true even with balanced data.  Thus, I would trust the F test values, for a couple of reasons.  It considers the effect after "factoring out" all the other effects--it is meant to be an orthogonal test to all the other tests.  The theory part of my brain seems to be non-functional this morning, but I know this difference should be apparent to me in the face of including interactions.  The agreement will probably get worse with imbalance, especially with a Kenward-Rogers adjustment.

Steve Denham

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
  • 1368 views
  • 0 likes
  • 2 in conversation