BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
ksb39
Calcite | Level 5

I am trying to test for differences in the ages at when two different behaviors emerge in nestlings birds. Age distributions of the two behaviors are reasonably normally distributed.  We have the ages of both behaviors recorded for each of 53 individuals, so I could use a paired t-test, but need to control for non-independence of sibling groups. The 53 nestlings come from 21 different nests, thus some nestlings were not independent (siblings had the same parents, habitat, etc.,).  Here is my code:

Proc GLIMMIX data=...;

CLASS nest ID behavior;

MODEL age = behavior / dist=poisson link=log solution;

RANDOM  nest;

RANDOM residual / type=vc subject=ID;

run;quit;

Does this seem appropriate?  Thanks!

1 ACCEPTED SOLUTION

Accepted Solutions
SteveDenham
Jade | Level 19

This should model what you are looking for:

Proc GLIMMIX data=...;

CLASS nest ID behavior;

MODEL age = behavior / dist=poisson link=log solution;

RANDOM  nest;

RANDOM behavior / residual type=vc subject=ID;

run;quit;

I hope this helps.

Steve Denham

View solution in original post

1 REPLY 1
SteveDenham
Jade | Level 19

This should model what you are looking for:

Proc GLIMMIX data=...;

CLASS nest ID behavior;

MODEL age = behavior / dist=poisson link=log solution;

RANDOM  nest;

RANDOM behavior / residual type=vc subject=ID;

run;quit;

I hope this helps.

Steve Denham

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

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