BookmarkSubscribeRSS Feed
Poultry
Calcite | Level 5

Hello all,

I've a data set including two fixed factors (house, strain) with 3 levels of each. The response variable is a lameness score of birds with scores ranging from 0 to 4 (0, 1, 2, 3, 4). I wanted to know if there is a effect of house and/or strain for each score separately. I considered using glimmix with following code

proc glimmix data=work;

class house strain;

model score=house strain/dist=multinomial link=glogit ;

estimate 'A vs B' house 1 -1/ bycat;

estimate 'A vs C' house 1 0 -1/ bycat;

estimate 'B Vs C' house 0  1 -1/ bycat;

estimate 'D vs E' strain 1 -1/ bycat;

estimate 'D vs F' strain 1 0 -1/ bycat;

estimate 'E vs F' strain 0  1 -1/ bycat;

run;

However Sas keeps giving me wrong degrees of freedom for main effects and i think the whole analysis was wrong. It would be of great help if somebody could tell me what i'm doing wrong or what would be the appropriate code...

thaks

1 REPLY 1
SteveDenham
Jade | Level 19

Well, assuming no interaction, something odd is happening.  By specifying link=glogit, you should be getting nominal data, and the estimates should be correct--if everything else is working.  So some design questions: I assume 3 houses and 3 strains.  How many observations are available?  What degrees of freedom are coming up in the output?  Maybe with those we can identify the source of the oddness.

Steve Denham

SAS Innovate 2025: Call for Content

Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!

Submit your idea!

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
  • 1 reply
  • 1478 views
  • 0 likes
  • 2 in conversation