In a paper on reintroduced cranes, I have data like this: ID, Sex, Year, Out* (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 = Out of reintroduction area. 0 = In reintro. area. I tried this: 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? Thanks, Richard
... View more