Programming the statistical procedures from SAS

Modeling paired data with GLMM?

Accepted Solution Solved
Reply
N/A
Posts: 1
Accepted Solution

Modeling paired data with GLMM?

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!


Accepted Solutions
Solution
‎03-06-2012 07:45 AM
Respected Advisor
Posts: 2,655

Modeling paired data with GLMM?

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


All Replies
Solution
‎03-06-2012 07:45 AM
Respected Advisor
Posts: 2,655

Modeling paired data with GLMM?

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

☑ This topic is SOLVED.

Need further help from the community? Please ask a new question.

Discussion stats
  • 1 reply
  • 200 views
  • 0 likes
  • 2 in conversation