topic How do I set up PROC GLIMMIX for simple biology problem? in SAS Procedures

How do I set up PROC GLIMMIX for simple biology problem?

r1231 — Fri, 25 Nov 2016 08:59:02 GMT

In a paper on reintroduced cranes, I have data like this: ID, Sex, Year, Dispersal* (e.g., 1-01, M, 2002, 1)

I want to show that females disperse more than males.

I used chi-square. Reviewer indicated this violated assumption of independence because of multiple years of data for same birds and that I should use GLIMMIX.

What are lines of code needed to set up model?

SAS version 9.x (latest)

*Dispersal: 1 = in reintroduction area; 0 = out (dispersed).

Thanks, Richard

Re: How do I set up PROC GLIMMIX for simple biology problem?

Ksharp — Fri, 25 Nov 2016 11:09:42 GMT

Or proc mixed. make ID as Subject variable, Year as a random effect. You should call @Steve or @lvm

Re: How do I set up PROC GLIMMIX for simple biology problem?

r1231 — Sun, 27 Nov 2016 12:50:56 GMT

id = identification code for each bird.

out = dispersal value (=1 if out of reintro. area); each bird has data in 1 or more years (observations).

This code:

proc glimmix;

title 'difference in dispersal by sex including all birds';

class sex id;

model out = sex / link=logit dist=binomial;

random int / subject=id;

run;

Produces these results: Type III Tests of Fixe Effects:
Effect = Sex, Num DF = 1, Den DF = 568, F Value = 5.26, Pr > F = 0.0221

Does this look correct? Thank you.

Re: How do I set up PROC GLIMMIX for simple biology problem?

Ksharp — Sun, 27 Nov 2016 13:04:30 GMT

Sorry. I am not expert about proc glimmix.
You should post it at Statistical Forum. 
@Steve @lvm ...   are there .