Help using Base SAS procedures

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

Reply
Occasional Contributor
Posts: 8

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

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

Super User
Posts: 10,041

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

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

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

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.

Super User
Posts: 10,041

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

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


Ask a Question
Discussion stats
  • 3 replies
  • 184 views
  • 0 likes
  • 2 in conversation