turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Proc Glimmix : Calculate intraclass correlation ...

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-26-2016 03:56 PM

Hello,

I would like to compute intraclass correlation using Proc Glimmix for a generalized linear mixed modelling (Poisson) with 1 fixed effect (intervention) and 1 random effect (cluster).

I model proportions with the Poisson model, defining log(n) in the offset term :

**proc** **glimmix** data=test method=quad;

class CLUSTER INTERVENTION;

model NEVENTS = INTERVENTION / link=log dist=poisson s offset=logn ;

random intercept / subject=CLUSTER G solution

I understand that the residual variance is the 'Estimate' listed in the "Covariance Parameter Estimates" table (or in the G matrix with only 1 row) but I have maybe some misunderstanding.

I can't find the common variance in order to estimate the ICC for this model: does it correspond to the 'Estimate' per cluster in the "Solution for random effects" ?

The ICC could then be calculated with ( common variance / common variance + residual variance ) but I'm not sure this is appropriate here.

I thought maybe others with more knowledge in GLMM and GLIMMIX may be able to enlighten me.

Regards,

Jude

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-26-2016 04:26 PM

You are specifying a Poisson distribution. There is **no** residual variance term with a Poisson, by definition. The variance estimate you are seeing is the cluster variance. If you fitted the model without a random effect term, you would not get any variance terms in the output.

The variance of the Poisson conditional distribution that you are using is given by the mean. There would be a different mean for each level of intervention.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-26-2016 04:51 PM

Thanks a lot LVM for your answer. I've got it.

To be sure to understand well, does it mean that I can't compute ICC with this distribution? Or is there a way to quantify the degree of similarity in the responses of individuals from the same cluster for the outcome (nevents) ?

Thanks again,

Jude

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

07-25-2017 07:23 PM

Totally worthless response from SAS. How does Stata generate ICCs then?